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Abstract 

Holographic Thermodynamics and Transport of Flavor Fields 

Andrew Hill O'Bannon 

Chair of the Supervisory Committee: 
Assistant Professor Andreas Karch 
Physics 

We use gauge-gravity duality to study a strongly-coupled non-Abelian gauge theory with 
flavor fields, i.e. fields transforming in the fundamental representation of the gauge group. 
We first study the thermodynamics of the flavor fields. In the grand canonical ensemble at 
zero temperature, we find a second-order transition when the mass of the flavor fields equals 
the chemical potential. We then study the transport properties of the flavor fields at finite 
temperature and density. We introduce external electric and magnetic fields and compute 
the resulting current of flavor charge. From this current wc extract the conductivity, using 
Ohm's law. In addition, wc compute the drag force on the flavor fields at large mass, in the 
presence of a finite baryon density and external electric and magnetic fields. 
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Chapter 1 
INTRODUCTION 

1.1 Motivation: QCD 

All known matter interacts via only four forces: gravity, electromagnetism, the weak nuclear 
force, and the strong nuclear force. Gravity is perhaps the most mysterious. Why is the 
cosmological constant so small? What is dark energy? How do we quantize gravity? The 
Standard Model of particle physics describes the remaining three forces, at least at energies 
below about lO-*^^ electron- Volts. Electromagnetism is extremely well-understood, thanks to 
the spectacular success of Quantum Electrodynamics (QED). The weak force is also well- 
understood, though big questions remain, for example about the exact nature of neutrinos 
and whether the Higgs exists. 

The strong nuclear force also remains mysterious. The accepted theory of the strong 
nuclear force is Quantum Chromodynamics (QCD). QCD is a non-Abelian gauge theory 
with gauge group SU{3). The theory thus has gauge bosons, massless spin-1 fields called 

"gluons," transforming in the adjoint representation of the gauge group. The theory also 
has matter fields, massive spin-1/2 fields called "quarks," transforming in the fundamental 
representation of the gauge group. 

QCD has "asymptotic freedom" [1,2]. An asymptotically free theory is a theory whose 
coupling strength decreases as the renormalization scale increases, and vice-versa. Physi- 
cally, in an asymptotically free theory, particles scattering at higher and higher energy will 
interact more and more weakly, while particles scattering at lower and lower energy will 
interact more and more strongly. In QCD, the renormalization scale at which the coupling 
becomes order one is called Kqcd- This is one of the fundamental energy scales of QCD, 
in addition to the quark masses. When discussing QCD, we will refer to "low-energy" and 
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"high-energy" regimes, where "low" and "high" are with respect to Aqcd, which is roughly 
200 million electron- Volts. 

Standard "textbook" methods for computing observables in quantum field theories rely 
on a valid perturbative expansion: assume the theory is almost free and calculate observables 
in a term-by-term expansion in a small parameter, namely the coupling. These methods 
were invented for QED, whose coupling remains perturbatively small over a very wide range 
of energies (far above the Planck scale, if we naively extrapolate up in energy). 

For QCD, on the other hand, perturbation theory is only reliable in high-energy regimes, 
due to asymptotic freedom. Theorists thus go looking for dials to turn to reach high energy. 
For example, we can introduce a temperature in QCD. At very high temperatures, far above 
^QCD, we expect excitations to scatter with energies on the order of the temperature, and 
hence to be weakly interacting. We can then ask questions and compute the answers reliably 
using perturbation theory. For instance, we can ask the questions of thermodynamics: what 
is the pressure of a thermal bath of quarks and gluons? The pressure of hot QCD has in 
fact been computed perturbatively to very high order in the coupling [3-5]. Similarly, we 
can introduce a finite chemical potential in QCD. For a chemical potential far above Aqcd 
(and far above the temperature), we expect to find a Fermi surface of quarks, and we expect 
excitations about that Fermi surface to scatter with energies on the order of the chemical 
potential. Perturbation theory should thus be reliable again, and indeed, has provided truly 
amazing insights into the nature of extremely dense quark matter [6]. These systems are 
not purely theoretical: the early universe was very hot and neutron stars are very dense. 

Of course, an enormous amount of physics occurs in the low-energy regime where QCD 
is strongly coupled and perturbation theory is no longer reliable. We know from experi- 
ments (for example, deep inelastic scattering) that at high energy the appropriate degrees 
of freedom of QCD are quarks and gluons, but at low energy the particles observed in ac- 
celerators are not quarks and gluons, but hadrons: baryons and mesons. We thus say that 
at low energy quarks are "confined" in hadrons, that is, a single, isolated quark has never 
been observed in vacuum, in dramatic contrast to, say, electrons or muons. 

The fundamental problem of QCD is easy to pose: given the high-energy theory, namely 
the path integral of QCD, derive the low-energy spectrum, i.e. the observed spectrum of 
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hadrons, as well as all of their interactions. To date, all attempts to solve this problem have 
failed. We still do not understand the dynamical origin of confinement in QCD. In one form 
or another, the obstacle to solving the problem is always the fact that the theory is strongly 
interacting at low energy. 

Perturbation theory is not the only tool theorists have, of course. A number of non- 
pcrturbativc techniques have been developed over the past few decades. One of the most 
versatile, and successful, is lattice QCD. In lattice QCD, spacetime is represented by a finite 
number of lattice sites, and quantum fields as degrees of freedom on that lattice^. The QCD 
path integral is then a finite number of finite-dimensional integrals. We can thus compute 
observables by "brute force," using a computer. Indeed, in this fashion we can compute any 
observable we want, at least in principle. 

The many non-perturbative methods developed for QCD can answer an enormous num- 
ber of questions about hadronic matter, but, unfortunately, they have yet to solve the 
fundamental problem of QCD. With lattice QCD, for example, we can compute the low- 
energy spectrum of QCD, but while a numerical simulation may give us an answer, it does 
not explain the dynamical mechanism of confinement. 

Moreover, even for less ambitious questions, many non-perturbative methods have tech- 
nical problems that prohibit their use. For example, the lattice is very good at answering 
questions about thermal physics of equilibrium states, but not very good at answering ques- 
tions about thermal physics out of equilibrium, that is, about time-dependent dynamics. 
Furthermore, lattice QCD also faces technical obstacles even for certain equilibrium states: 
currently, no one knows how to simulate a finite density of quarks in lattice QCD, for in- 
stance. The problem in both cases is essentially the same. Lattice simulations are done in 
Euclidean signature via Monte-Carlo sampling of field configurations, with weighting factor 
e~'^, with S the Euclidean QCD action. Real-time dynamics requires Lorentzian signature, 
hence the weighting factor becomes e*'^ with S the Lorentzian-signature action, and sim- 
ilarly a chemical potential introduces a factor of i into the action, so that in both cases 
the weighting factor oscillates, making Monte-Carlo impossible in practice (though, strictly 



Some standard textbooks on lattice QCD are refs. [7-10]. 
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speaking, not in principle)^. 

The situation is especially dire given the recent experimental results from the Relativistic 
Heavy Ion Collider^ (RHIC). The basic result of RHIC was that, for temperatures around 
250 million electron- Volts, the best description of QCD is ideal hydrodynamics [12-17]. 
In other words, the simulations that best fit (certain appropriate) observables were not 
performed in lattice QCD, but rather in hydrodynamic models, and in fact in hydrodynamic 
models where the shear viscosity was set to zero (hence "ideal"). 

That hydrodynamics was a good description is not difficult to understand. Many people 
were surprised that ideal hydrodynamics worked so well, though. To explain why, we need 
to review some of the physics of hydrodynamics. 

Prom a modern perspective, hydrodynamics is just an effective field theory. In the ocean, 
the mean free path of a single water molecule is on the order of nanometers, while ocean 
waves have wavelengths on the order of meters or longer. If we only want to ask questions 
about the physics of the low-energy degrees of freedom, the water waves, then a description 
using the microscopic degrees of freedom, the water molecules, will probably not be easy to 
use, given that the interaction between water molecules may be very complicated. 

We have a large separation of scales, however, so we can follow the rules of effective 
field theory and write equations of motion^ including every possible term consistent with 
the symmetries. For liquid water the symmetries are rotation and translation invariance. 
The effective theory will be valid up to some UV cutoff, that is, down to some length scale 
somewhere between a meter and a nanometer. Higher-derivative terms will be suppressed 
by inverse powers of the cutoff, so we can keep leading terms up to some order in the cutoff, 
giving us a theory with a manageable (read finite) number of terms. All these terms will 
a priori have undetermined coefficients, however. In hydrodynamics these coefficients are 
called transport coefficients, an example being the shear viscosity, ry. We cannot compute 

^In the actual practice of lattice QCD, the fermions are integrated out and hence represented as a 
determinant in the path integral. This just shifts the problems around, though. Introducing a chemical 
potential, for example, makes the fermion determinant complex, producing the so-called "sign problem," 
which is equivalent to what wc have described. 

good pedagogical review of heavy-ion physics is rcf. [11]. 

*We do not start with an effective Lagrangian because of dissipation. 



5 



anything with our effective theory until we determine these coefficients. In principle, we 
can compute the transport coefficients from a microscopic theory. For most real materials, 
though, such as water, the microscopic theory is very complicated, so no one would even try 
to compute the transport coefficients. Fortunately, in the real world, we can usually just 
measure transport coefficients, and then use the effective theory to make predictions. 

The fact that hydrodynamic modeling of RHIC collisions works at all indicates that 
the mean free path inside the RHIC plasma must be much smaller than the size of the 
plasma [15,17], which is about"'' 10^^"^ meters. A short mean free path in turn suggests that 
the system is strongly-coupled [15,17], so that an excitation will not travel very far before 
experiencing a collision. That the RHIC plasma is strongly-interacting is not surprising, 
given that RHIC produces a plasma at temperatures on the order of Aqcd, where we expect 
the coupling to be order-one anyway. 

The fact that ideal hydrodynamics worked so well was totally unexpected, however. 
Before the experiments at RHIC, the only estimates for the viscosity came from perturbation 
theory [18,19]. The actual visocsity appears to be much smaller than these estimates. The 
big question that arose from RHIC was thus: why does the viscosity of QCD appear to be 
so small at RHIC temperatures? 

Theorists would love to answer this question by computing the viscosity from the under- 
lying microscopic theory, QCD. The problem, once again, is strong coupling. Perturbation 
theory will not be reliable. RHIC collisions involve real-time, out-of-equilibrium physics 
that lattice QCD cannot reliably address. In short, as this point in history, no one knows 
how to compute the viscosity of the RHIC plasma from QCD! 

That could have been the end of the story. At this point, however, like a bolt out of the 
blue, quantum gravity enters the picture, via the anti-de Sitter / Conformal Field Theory 
(AdS/CFT) correspondence [20]. Let us pause to introduce AdS/CFT^, and then discuss 
what it has to do with the viscosity of QCD. 

^At RHIC, gold nuclei are used in collisions. A gold nucleus is about 15 x 10~^^ meters in size, so that 
their collision produces a plasma of the size quoted above. The exact size in a particular collision depends 
on how much the gold nuclei overlap when they collide {i.e. on the impact parameter, or equivalently the 

"ccntrality" of the collision). 

''Two standard reviews of the AdS/CFT correspondence are refs. [21,22]. 
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The original AdS/CFT conjecture was that a conformal field theory, = 4 super- 
symmetric SU{Nc) Yang- Mills theory (SYM), is in fact equivalent to a theory of quantum 
gravity, type IIB string theory, formulated on background geometries that asymptotically 
approach AdS^ x S^, where AdS^ is (4-1-1 )-dimensional anti-de Sitter space and is a 
five-sphere [20]. Anti- de Sitter spaces are defined as spaces of constant negative curvature. 
They are naturally equipped with a radial coordinate, and also with a boundary at some 
asymptotic value of that radial coordinate. The SYM theory in some sense "lives" on this 
(3+l)-dimensional boundary. We will review = 4 SYM theory and AdS spaces below, in 
sections 2.1 and 2.2. 

What do we mean that the two theories are "equivalent?" AdS/CFT is the statement 
that the two theories. A/" = 4 SYM theory and type IIB string theory on AdS^ x S^, are 
identical: they are merely two different descriptions of the same physics. If this is true, 
then all of the physics of one description must map onto all of the physics of the other 
description somehow. We will review the "dictionary" of how to translate between the two 
theories in section 2.4. 

What makes AdS/CFT useful is that it is a duality: when one theory is weakly-coupled, 
the other is strongly-coupled, and vice- versa. AdS/CFT is thus especially useful for answer- 
ing questions about the A/" = 4 SYM theory at strong coupling. To see this explicitly, we 
start with two key entries in the AdS/CFT dictionary [20], 

2 

— ^ = A-rrgsNc, 27r gg = gyn 

Here L is the radius of curvature of both the AdS^ and factors of the background, gg 
is the string coupling, N^. is the number of colors of the AA = 4 SYM theory, and gvM is 
the coupling of the H = 4, SYM theory^, a' is related to the tension, T^, or energy per 
unit length, of strings in string theory by Tg = Xjl-Koi . The string tension is a fundamental 
dimensionful scale in string theory, which, via some dimensional analysis, gives us a length 
scale, the string length Ig = \fo^ . Notice that not all string theory strings have the length 

^We follow the conventions of D-brane physics (which we review in section 2.3.2), for which the definition 
of gYM is actually half the conventional value. The discrepancy comes from the normalization of the 

generators of the gauge group (in our case, SU{Nc)): Tr{TaTt) = dSat- The usual choice in D-brane 
physics is d = 1, whereas the conventional choice in field theory is d = 1/2. 
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is, rather, it is a characteristic scale at which the finite length of strings becomes important. 

Currently, only perturbative string theory is under good theoretical control, so we want 
to take 5s ^ 1. At leading order, we then have classical string theory, in the sense of "tree- 
level in the string genus expansion." If we keep Nc fixed as we take ^ 1, then we see 
that <^ a'^, meaning that the background is very highly curved: the radius of curvature 
is smaller than the string scale! 

Currently, no one knows the full string spectrum for strings propagating in AdS^ x S^, 
because string theory is difficult to quantize in this background, but dimensional analysis 
tells us that the mass-squareds must be proportional to a'~^. When <C a'^, the massive 
string modes will be very light relative to the natural mass scale of the background, L~^. 
Since we do not know the spectrum of these modes, we would like to avoid this limit. 

We are thus motivated to take Nc ^ oo as we take — in such a fashion that ggNc 
remains fixed, and then take an additional limit in which ggNc S> 1, so that S> a'^. 
The massive string states will then be very heavy, and we can justifiably write an effective 
theory of masslcss modes alone. In this classical, low-energy limit, the effective theory of 
type IIB string theory is known: it is type IIB supergravity. 

Notice what this limit translates into in the gauge theory: the 't Hooft limit [23], with 
large 't Hooft coupling. In other words, we define the 't Hooft coupling A = igyj^^N^ 
which is the physical coupling of the gauge theory. We then take two limits. First, we 
take Nc ^ oo with gyj^ — such that A is fixed. Second, we take a strong-coupling limit, 
A — GO. AdS/CFT is telling us that this strongly- coupled gauge theory is in fact equivalent 
to a classical theory of gravity. We thus have a new method for computing observables in 
strongly-coupled gauge theories: classical gravity. 

The first objection, of course, is that = 4 SYM is not QCD! Indeed, the two theories 
are vastly different. A/' = 4 SYM has much more symmetry than QCD, namely supersym- 
metry and conformal symmetry. Supersymmetry is why J\f = 4 SYM has many more fields 
than QCD, namely fermions and scalars in the adjoint representation of the gauge group (in 
a supermultiplet with the gluons). Conformal symmetry is why, where the QCD coupling 
shrinks with increasing energy scale, the coupling of TV = 4 SYM does not change at all 
as the energy scale changes. Furthermore, in TV = 4 SYM, as in any conformal theory, 
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we cannot even define an S-matrix, or even excitations that we can call particles, since we 
cannot define asymptotic states. 

The original AdS/CFT conjecture, however, motivates the more general conjecture of 
"gauge-gravity duality," which is the statement that any theory of quantum gravity (it 
need not be a string theory), if formulated on a spacetime that approaches an AdS space 
asymptotically, is equivalent to a gauge theory formulated on the boundary of that space [24] . 
In other words, gauge-gravity duality proposes that many theories of quantum gravity are 
identical to gauge theories in one lower dimension, again in the sense that the two theories 
will merely be two different descriptions of the same physics. 

Gauge-gravity duality systems are concrete examples of "holography" [25-27], which, 
roughly speaking, is the statement that the physics of any theory of quantum gravity in some 
region of spacetime must be described by degrees of freedom formulated on the boundary 
of that spacetime. In other words, any theory of quantum gravity may be described by 
some degrees of freedom "living" in one lower dimension, hence the name "holography." A 
big hint for holography comes from black hole physics, where the entropy of a black hole is 
proportional to the area of its event horizon, rather than its volume. Notice that holography 
is not constructive, however, asserting only the existence of some lower-dimensional degrees 
of freedom, but saying nothing about their dynamics. Gauge-gravity duality tells us that, 
in some cases, the lower-dimensional degrees of freedom arrange themselves into a gauge 
theory. 

Gauge-gravity duality is an extremely profound statement, raising deep questions about 
both gauge theories and quantum gravity. How do certain gauge theories, which include no 
dynamical gravity whatsoever, "know" about quantum gravity, and, on top of that, emergent 
extra dimensions? Where in the Hilbert space of a gauge theory are the microstates of a 
black hole? If gauge theories are unitary, does that mean black hole evaporation must 
be unitary? (By now even Stephen Hawking concedes that the answer is yes, because of 
AdS/CFT [28].) Indeed, many people believe that, "Gauge-gravity duality is the biggest 
conceptual change in quantum mechanics since the 1930's" [29]. 

Gauge-gravity duality is wonderful and good, but in the interest of full disclosure we 
present a caveat: it has not been proven to be true. In other words, nothing like a rigorous 
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mathematical proof of an equality between theories exists , and hence gauge-gravity duality 
remains a "conjecture." Nevertheless, an enormous amount of evidence has accumulated 
suggesting that gauge-gravity duality is true. For a partial list of such evidence, see ref. [21]. 
The growing mountain of evidence has laid to rest any doubts about the validity of the 
conjecture. 

= 4 SYM docs not describe the real world at all, of course, so we would like some 
way to determine what theories have holographic diials. In particular, wc would like to 
know whether QCD has a holographic dual. If not, then gauge-gravity duality is just a 
theoretical exercise. Most examples of gauge-gravity duality come from D-brane physics in 
string theory, which is one method to generate specific dualities (gauge theory X is dual 
to string theory Y). If someone hands you an arbitrary gauge theory, however, and asks 
whether it has a dual description as a theory of quantum gravity, how do you answer the 
question? What tests could you do? These are open problems; a complete classification of 
which theories have dual descriptions does not yet exist. 

A key clue, though, is the large-A^c limit. For theories with a well-defined large-iVc 
limit, 't Hooft noticed that the Feynman diagrams of the theory, when drawn in a double- 
line notation (basically meaning one line per gauge index), naturally arrange themselves in 
powers of Nc [23]. The power of Nc is, in fact, the genus of the Riemann surface on which 
the diagram can be drawn without any lines crossing: the leading diagrams are order N'^ 
and can be drawn on a sphere, the sub-leading diagrams are order iV° and can be drawn 
on a torus, and so on. In other words, the large-iVc expansion is shockingly similar to the 
genus expansion in perturbative string theory, if we identify 1/Nc with gg. The original 
AdS/CFT conjecture provides an explicit example where the large- AT^ genus expansion is, 
in fact, a perturbative string expansion: using the dictionary, and with L^/a'^ fixed, we 
find Qs oc l/Nc- 't Hooft discovered the \axge-Nc genus expansion over thirty years ago, but 
it remains the biggest piece of evidence suggesting that gauge theories with a well-defined 

*A good question is what would really eonstitute such a proof. Before attempting to prove that two 
theories are equivalent, however, a good idea is to define the two theories involved. We can in principle 
provide a lattice definition for a quantum field theory, but defining a theory of quantum gravity is much 
harder. Given this state of affairs, some people have actually proposed that gauge-gravity duality itself 
could serve as the definition of the theory on the gravity side of the correspondence [20]. 



10 



large-A^c limit may have dual descriptions as theories of quantum gravity. In particular, 
this is one of the main reasons why QCD is believed to have a dual description. 

The holographic dual of QCD will be hard to find, if it exists, however. Perhaps the 
biggest obstacle is that QCD has no large separation of scales. In the original AdS/CFT 
conjecture, the gauge theory was conformal, meaning in particular that the coupling A was 
a free parameter, so we could, quite artificially, dial A to be as big as we like. On the string 
theory side, this made massive string states arbitrarily heavy and the background geometry 
arbitrarily weakly curved, so that classical gravity was then a good approximation to the 
string theory. In QCD, however, the coupling is not a free parameter: it is determined 
dynamically. We cannot artifically produce a large separation of scales. We will be forced 
to deal with a fully quantum theory of gravity, with high spacetime curvatures, many excited 
states (the massive string states), etc. 

What should we do, then? One option is to forge ahead and attempt to construct the 
holographic dual of QCD. We can take a "top-down" approach, deriving from a full quantum 
gravity theory, such as a string theory, a geometry whose dual is QCD, or perhaps whose 
dual approaches QCD in some limit. An example of such a top-down approach appears in 
ref. [30]. Another approach is "bottom- up," that is, demanding that the gauge theory be 
QCD and inferring the holographic dual [31,32]. 

Yet another approach is to abandon QCD and work with a theory whose dual we know 
precisely, namely A/" = 4 SYM. Given how different J\f = 4 SYM and QCD are, though, 
what could this approach possibly teach us about QCD? 

Our philosophy is this: try to make "universal" statements that may be applicable to 
QCD and its holographic dual. These "universal" lessons come in two flavors: qualitative 
and quantitative. 

What are some examples of qualitative universal statements? A state in a gauge theory 
with nonzero entropy will be dual to a geometry with a black hole horizon: the black hole's 
Hawking temperature and entropy are identified with the gauge theory temperature and 
entropy [33,34]. Topology changes in the geometry, for example a Hawking-Page transi- 
tion in which a black hole condenses, generically map to first-order phase transitions in 
the gauge theory, for example a deconfinemcnt transition: this makes sense since on both 
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sides a discontinuous change occurs [34]. As we can see, qualitative universal statements 
usually involve translating gauge theory questions into gravity questions. What does a fi- 
nite quark density look like holographically [35]? What does a quark moving through a 
gauge-theory plasma look like holographically [36, 37]? What does gluon scattering look 
like holographically [38,39]? 

What are examples of quantitative universal statements? Here is where we return to 
the RHIC story. The canonical example of a quantitative universal statement concerns the 
ratio of shear viscosity, rj, to entropy density, s. Hydrodynamics is the correct effective 
description of low-energy exictations of finite-temperature, strongly-coupled = 4 SYM 
for essentially the same reasons as in QCD. In particular, at strong coupling the mean 
free path of excitations will be short. We may not know how to compute rj for QCD at 
RHIC temperatures, but we do know how to compute 77 for A/" = 4 SYM theory using 
AdS/CFT [40-43]. The result is usually expressed as the ratio rj/s = I/Att. This value 
is far smaller than that of any known substance. The big surprise here, though, was that 
the value rj/s = 1/4tt is actually the result for any theory with a gravitational dual [43-45]. 
More precisely, l/47r will be the value of r//s in the limit of any gauge theory for which the 
holographic dual is classical gravity, so that the dual geometry includes a classical black 
hole (finite-coupling corrections, for example, which translate into curvature corrections in 
gravity, may alter the value of ry/s). The smallness and universality of jy/s are suggestive, 
given that QCD has a very small viscosity and, most likely, a holographic dual. Though 
neither a prediction nor an explanation of the RHIC results, the precise number rj/s = I/Att 
serves as an order-of-magnitude starting point for thinking about QCD. Universal quantities 
such as r]/s are hard to find, however. 

In the big picture, then, the goal is not to compute observables in = 4 SYM theory 
and then make a prediction about QCD. Comparing quantitative results between A/" = 4 
SYM theory and QCD is a very difficult task, and probably not fruitful. Instead, we look 
for qualitative universal statements that should be true for the holographic dual of QCD, so 
that if today someone handed us the holographic dual of QCD and asked us to compute the 
phase diagram or the viscosity, we would at least know where to start: a state of nonzero 
entropy is dual to a black hole, deconfinement is dual to black hole condensation, etc. In 
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other words, we develop methods for computing gauge theory observablcs using gravity. 
Along the way, we try to extract quantitative universal statements, which may teach us 
how much of QCD physics is in fact insensitive to the details of its dynamics, that is, 
how much of QCD physics is in fact the same as the physics of any gauge theory with a 
gravitational/holographic dual. 

We summarize our whole philosophy by saying that gauge-gravity duality is a tool for 
studying strongly-coupled gauge theories. We do not ask questions about why the duality 
works. We do not ask questions about quantum gravity and compute the answer using 
the gauge theory. Instead, we begin with a question in a strongly-coupled gauge theory, 
translate to gravity, perform the calculation, and then translate back to gauge theory. In 
this spirit, AdS/CFT is a "black box," and A/" = 4 SYM serves as a "toy model" of QCD. A 
good name for this approach, one that is only partly facetious, is "Applied String Theory." 

In this dissertation we focus on questions in strongly-coupled A/" = 4 SYM theory that 
would be difficult to answer for QCD, using existing methods. First, we will study the 
theory at zero temperature and finite quark density. Second, we will compute a transport 
coefficient, a conductivity. We will not attempt to compare our results to QCD. As we have 
tried to argue, the method is more important than the result! 

Hopefully we have convinced the reader that this line of research is worth pursuing, and 
in the process justified the existence of this dissertation. The other main lesson we hope to 
impart in the sequel is that AdS/CFT is not only useful, it's also fun! 

1.2 Outline 

In chapter 2 we review the original AdS/CFT correspondence. We begin with a review of 
J\f = 4: SYM theory, and then of AdS spaces. We then repeat, in summary, Maldacena's 
original D-brane construction of AdS/CFT [20], and how to compute correlators in the 
gauge theory from the gravity theory [24,46]. We also review the finite-temperature version 
of the correspondence [33,34]. 

In chapter 3 we review how to introduce fundamental-representation "flavor" fields, into 
the AdS/CFT correspondence, using probe D7-branes [47,48]. The on-shell action of these 
D7-branes is infinite: it must be regulated and renormalized. We review the procedure of 
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"holographic renormahzation" of probe D7-branes in AdS^ x S^, which was first described 
in ref. [49] . 

In chapter 4 we turn to the thermodynamics of the flavor fields. We review the known 
phase transitions associated with the flavor fields at finite temperature [50-56], finite den- 
sity [35,57,58], and with finite electric and magnetic fields [59-63]. In all of these cases, the 
calculations on the supergravity side were numerical. The point of this chapter is to present 
an exact supergravity solution representing a state of the gauge theory with zero tempera- 
ture but a finite density of flavor fields, as described in ref. [64]. Using this description, we 
find a second order transition, in the grand canonical ensemble, when the mass of the flavor 
fields equals the chemical potential. 

In chapter 5, we turn to the transport properties of the flavor fields. We present su- 
pergravity solutions representing a state of the gauge theory with finite temperature, finite 
density and finite electric and magnetic fields. The external fields produce a current of 
flavor fields. We compute this current and from it extract a conductivity tensor. Along the 
way, we compute the drag force on the flavor fields. These results were first presented in 
refs. [65,66]. 

We conclude in chapter 6 with open questions and suggestions for future research. 

IMPORTANT NOTE ABOUT UNITS: We only write the AdSrj radius of curva- 
ture, L, explicitly in chapter 2. In the subsequent chapters, we measure lengths in units of 
L, or in other words, we set the radius of AdS^ to be one: L = 1. An important conversion 
from supergravity quantities to field theory quantities is then a'"^ = A. 
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Chapter 2 

THE ADS/CFT CORRESPONDENCE 

In this chapter we review the AdS/CFT correspondence. We first review the ingredients 
that enter into the correspondence: J\f = 4 SYM and AdS spaces, with emphasis on physics 
rather than equations, although we present a few AdS metrics exphcitly. We then recaU 
some very basic facts about strings and D-branes, so that we can review Maldacena's original 
"derivation" of the correspondence, using D3-branes. We review the precise statement of 
the correspondence, as an equivalence of partition functions, and then review the finite- 
temperature version of the correspondence. Most of the material in this chapter can be 
found in the AdS/CFT reviews refs. [21,22] and the string theory textbooks refs. [67-70]. 

2.1 J\f = 4 Super-Yang-Mills Theory 

In this section we will review some important properties of = 4 SYM theory. Wc will only 
consider the theory in fiat (3-|-l)-dimensional space, with mostly-plus Minkowski metric rjni, 
(as opposed to, say, the theory formulated on some curved manifold). 

Supersymmetry is an extension of the Poincare algebra to include anti-commuting sym- 
metry generators. These anti-commuting generators are spinors, called "supercharges." Let 
Af denote the number of left-handed Weyl spinor supercharges. Weyl spinors include two 
components, i.e. two complex numbers. The total number of real supercharges is thus 
2 X 2 X J\f. The super-Poincare algebra is invariant under rotations of the supercharges into 
one another. These rotations form a group U{M)r, called the R-symmetry group (hence 
the subscript R). 

In a supersymmetric quantum theory, the supercharges act on the Hilbert space as 
raising and lowering operators for helicity. If J\f becomes bigger than 4, then the theory 
will have states of helicity greater than 1. Some theories have Af > 4, and include fields of 
higher helicity. An example is TV = 8 supergravity, in which one field is the metric, with 
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helicity 2. 

We are interested in theories with states of heUcity < 1 only. More specifically, we are 
interested in gauge theories without dynamical gravity, hence A/" = 4 is the maximal number 
of supersymmetries we can allow. The maximally supersymmetric gauge theory is A/" = 4 
supersymmetric Yang-Mills theory {M = 4 SYM). The theory of course describes massless 
spin-1 fields, the gauge fields. Supersymmetry then demands that fields of lower spin also 
be included. The = 4 gauge supermultiplet includes the gauge fields, as well as four 
left-handed Weyl fermions and six real scalars (or equivalcntly three complex scalars). By 
supersymmetry, all must transform in the same representation of the gauge group, namely 
the adjoint representation, and all must have the same mass. Gauge invariance forbids a 
mass for the gauge fields, hence the fermion and scalar fields are massless. 

Remarkably, the Lagrangian for A/" = 4 SYM theory is unique! That is, if we demand 
that the theory be renormalizable, then gauge invariance and supersymmetry completely 
determine the form of the action. Contrast this with Af = 1 supersymmetric theories, for 
example, where we have a lot of freedom in choosing the (super)potential. 

We will not reproduce the AA = 4 SYM Lagrangian. It is written explicitly in ref. [22]. It 
includes the usual kinetic term and theta term for the gauge fields, the usual kinetic terms 
for the fermions and scalars (with gauge covariant derivatives), Yukawa couplings between 
the scalars and fermions, and quartic potential terms for the scalars. Supersymmetry of 
course restricts the interaction terms, and A/" = 4 supersymmetry is so restrictive that only 
one coupling, the gauge coupling gvM, appears in the J\f = 4 SYM Lagrangian. In other 
words, qym is not only the strength of the interaction of the fermions and scalars with 
the gauge fields, but also the strength of the Yukawa interactions and the scalars' quartic 
self-interactions. 

Classically, AT = 4 SYM is scale-invariant. The coupling gyM is dimensionless, so to 
specify the classical M = 4, SYM theory we need only specify the gauge group, gyM and 
the theta angle. (Starting now, we will ignore the theta angle.) We do not need to specify 
any dimensionful scale. Of course, the same is true for classical pure Yang-Mills theory, as 
well as classical Yang-Mills plus massless quarks. 

Any textbook on quantum field theory will tell you, however, that to define a quantum 
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field theory we must specify tlie value of tlie coupling at some renormalization scale. The 
dynamics of the theory then determines how the coupling strength changes with energy 
scale, that is, how the coupling "runs." In particular, at some characteristic scale the 
coupling will run to a non-perturbative, i.e. order-one, value. Notice what this means for 
any theory that is classically scale-invariant: quantum effects dynamically generate a scale. 
The classical symmetry, scale invariance, is broken in the quantum theory. Such theories 
thus have "anomalous scale invariance," or simply a "scale anomaly." Pure Yang-Mills and 
QCD have scale anomalies. In QCD, the dynamically generated scale is called Aqcd- 

Amazingly, = 4 SYM theory remains scale-invariant even quantum mechanically! 
That is, once we specify the value of qym at one scale, the coupling does not change: gyM 
will have the same value at any other scale. In other words, M = 4, SYM has no scale 
anomaly, so nothing like Aqcd is generated. 

In fact, scale invariance is part of a larger symmetry: conformal symmetry. A conformal 
transformation acts on the metric as Qf^^ 0.{t, x)^ (7^1, where fl{t, x)"^ is a smooth, positive 
function (which is why we write it as a square). 

A conformal field theory (CFT) is a theory that is invariant under conformal transfor- 
mations. In (H-l) dimensions, theorems exist showing that, if a theory is scale invariant, 
it must in fact be conformally invariant [71,72]. In other words, in (H-l) dimensions, scale 
invariance implies full conformal invariance. In (3-1-1) dimensions, no such theorem exists 
(yet), but in all known examples scale invariance is accompanied by the full conformal group. 
An example is A/" = 4 SYM. 

What is the full symmetry group of A/" = 4 SYM theory, then? In brief terms, we can 
say that the Poincarc group is enhanced to the conformal group, which is 50(4,2). More 
explicitly, the bosonic symmetries are, for spacetime coordinate x'^, 

• Translations: + c for real constant c 

• Lorentz transformations: x^^ — ^ A'^n x'^ for Lorentz transformation matrix A^^, obeying 
Va/sA'^nA^u = r]^y 

• Scale transformations: x^ ^ Xx^ for real, constant, positive A 



17 



• Inversions: x'^ —>■ ^ with |a;p = g^i^x^x' 



The scale transformation and inversions are the "extra" transformations that take us from 
the Poincare group to the full conformal group. The theory of course also has the R- 
symmetry U{A)ji. Locally, the SU{4)ji subgroup is isomorphic to 50(6). The fermions of 
the theory transform as a 4, and the scalars as an antisymmetric 6, of the SU{4)ji. 

J\f = 4 SYM theory also has fermionic symmetries, of course, namely the sixteen real 
supercharges of A/" = 4 supersymmetry. The Poincare group is extended to include these 
supercharges as well as the generators of conformal symmetry, so the complete symmetry 
algebra is called the "superconformal" algebra. If that were all, however, the superconformal 

algebra would not close, so in fact sixteen additional fermionic generators arc required. 
These sixteen additional fermionic generators are called superconformal generators. 

Conformal symmetry severely restricts observables. For example, in a CFT, the two- 
point functions of operators of definite scaling dimension are completely fixed, and the 
three-point functions of such operators are fixed up to some overall constants. An S-matrix 
cannot be defined in a CFT. In scattering, we begin with some infinitely far-separated wave 
packets (particles) that we "shoot in" to some region where the interaction occurs, and 
we then measure some infinitely far-separated wave packets that come out. With scale 
and inversion invariance, however, in which we can exchange the interaction point with 
the point at infinity, we cannot define infinitely far-separated wave packets, and hence we 
cannot define an S-matrix. 

Conformal symmetry does have an upside, though: in any problem, the only scales 
will be those that we introduce ourselves, which can greatly simplify things. For example, 
suppose we study a CFT at finite temperature (and in flat space), and want to find any phase 
transitions. Without doing any work, we know the answer: the theory has no dimensionful 
scale that can set a transition temperature, hence no finite-temperature phase transitions 
can occur] 
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2.2 Anti-de Sitter Space 

In this section we review properties of anti-de Sitter (AdS) spaces. We will not attempt to 
be mathematically rigorous, rather, we hope to provide simple, intuitive pictures sufficient 
for understanding physics. 

de Sitter spaces are defined as spaces of constant positive curvature. Anti-de Sitter 
spaces are defined as spaces of constant negative curvature. In physics language, AdS 
spaces are solutions of Einstein's equations with negative cosmological constant. 

We will be interested in (4-|-l)-dimensional anti-de Sitter space, AdS^. We know we can 
get a space of constant negative curvature from a hyperboloid, so we start in one higher 
dimension and write AdS^ as a hypersurface. That is, we begin with the algebraic equation 
for a hyperboloid in M^, 

4 

x2 + X|-^X| = l2 (2.2.1) 

i=l 

where G R^. Wc can already sec from the signs on the left-hand side that AdS^ has a 
symmetry S'0(4, 2), which will be the isometry group of the AdS^ metric. Here L is the 
radius of curvature of the hyperboloid. 

One way to parameterize a solution of eq. (2.2.1) is to let 




(2.2.2) 



where i = 1, 2, 3, Xi £ M^, and n > 0. These coordinates only cover half of the hyperboloid 
and are thus called "Poincare patch" coordinates. Coordinates that cover the entire hyper- 
boloid do exist, the so-called "global" coordinates, but these will not be important in the 
sequel. In Poincare patch coordinates, the induced metric of the hyperboloid is 

d^A,s = L'^ + ^ Ht' + dx') (2.2.3) 

An explicit calculation of the Ricci scalar gives — so the space has constant negative 
curvature with radius of curvature L, as advertised. 

We can see that the space is essentially (3-|-l)-dimensional flat space plus an extra 
"warped" direction, u. The coordinate u is the radial coordinate of AdS^. For any fixed 
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value of u, we find a (3+l)-dimensional hypersurface whose metric is simply that of M^'^ 
times a conformal factor, u^/L^. A simple intuitive way to think of AdS^ is thus as an 
onion, with each layer of the onion skin being fiat R^'^. As u shrinks, however, the layer of 
the onion skin shrinks as u^/L^. We depict AdS space (crudely) in figure 2.1 (a.). In the 
figure, the vertical direction is u, with u = at the bottom and u = oo at the top. Wc have 
also drawn two spatial Minkowski directions and suppressed all other directions, and drawn 
a "slice" of the space at some finite value of u. 

As u — ^ we find a horizon, since gu — ^ 0. This horizon is called the Poincare horizon. 
The Poincare horizon has zero area because 5^.^. as u ^ 0. Nothing dramatic happens 
in crossing the Poincare horizon: on the other side of the horizon is the other Poincare 
patch, covering the other half of the hyperboloid. In figure 2.1 (a.), the Poincare horizon is 
the point at the bottom. 

The metric of any AdS space will always have a second-order pole for some value of the 
radial coordinate, i.e. the metric will always diverge quadratically. The place where this 
occurs is called the boundary of AdS. In our coordinates the boundary is at u 00. In 
figure 2.1 (a.) the boundary is at the top. 

We can extend the full AdS5 metric to the boundary, that is, we can extract a finite 
boundary metric from the AdS^ metric, as follows. We must first choose a "defining func- 
tion,"' which is a function of u, t and x that must have a second-order zero at the boundary. 
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and otherwise just needs to be smooth and positive. An example is f{u,t,x) = Q{t,x)'^^. 
If we multiply the Minkowski part of the metric by such an f{u, t, x) and then take the limit 
u — oo, we recover a finite boundary metric, = Q,'^{t,x) (— di^ + dx^). A particular 

choice of f{u,t,x) gives a particular boundary metric, but any other choice is equally ac- 
ceptable, so the boundary metric is only defined up to a conformal factor. In other words, 
the bulk metric only determines a conformal class of boundary metrics. 

The boundary (u oo) is in fact an infinite geodesic distance from any point in the bulk 
(any finite u). Physically, AdS spaces are gravitational potential wells. A massive object 
sitting at some value of the radial coordinate will naturally fall away from the boundary, 
which in the Poincare patch means towards the Poincare horizon. 

Another metric that we will need is the AdS-Schwarzschild metric. As the name implies, 
AdS-Schwarzschild is a geometry that has a black hole horizon and that asymptotically 
approaches an AdS geometry. The AdS-Schwarzschild metric is 

where h{u) = 1 — uf^/u'^ and = ttTL^ is the location of the horizon, with T being the 
Hawking temperature of the black hole. AdS-Schwarzschild is again a solution of Einstein's 
equation with negative cosmological constant, and with the same asymptotic boundary 
condition as pure AdS space. We depict AdS-Schwarzschild in figure 2.1 (b.). Notice that 
at u = oo the space looks the same as in figure 2.1 (a.). The space is now cut off at the 
horizon, at u = Uh, which now has finite area, as we have crudely depicted in the figure. 

In what follows, another coordinate system will be very useful: so-called Fefferman- 
Graham coordinates [73]. These are coordinates in which we define a new radial coordinate, 
z, such that g^z is always just 1/z'^, and the boundary is always at z = 0. Explicitly, we 
change coordinates from u to 

-1/2 



z = V2L'^ \^u^ + - ulj (2.2.5) 
The horizon is then at z^ = \/2L^/u/i = s/^/ttT, and the metric becomes 
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The funny factors of \pl guarantee that when T = we have simply z = I? ju. The Poincare 
horizon is then at z — ^ oo. 

Having learned a little about the symmetries of A/" = 4 SYM and about AdSf,^ we can 
already see hints of the AdS/CFT correspondence. The A/" = 4 SYM theory has an 50(4, 2) 
symmetry, which is precisely the isometry group of AdS^. The metric of AdS^ only defines 
a conformal class of boundary metrics. Any physics that occurs on the boundary must be 
invariant under the choice of defining function, that is, must be conformally invariant. The 
AdS/CFT conjecture is, in fact, that a specific conformal field theory, AA = 4 SYM, "lives" 
on the boundary. We can motivate the correspondence much better by studying strings and 
D-branes, though, so now we will quickly review some of their physics. 

2.3 Strings and D-branes 

Most of the material in this section can be found, in far greater detail, in standard string 
theory textbooks [67-70]. 

2.3.1 String Basics 

In its textbook formulation [67-70], string theory supposes that all fundamental particles 
are, in fact, not point particles but extended objects, namely one-dimensional objects called 
strings. Of course, we have yet to see evidence for such strings: to date, as far as we 
know, electrons, muons, neutrinos, etc., really are point particles. To be consistent with 
experiment, then, the length of strings must be very small, far smaller than any distance 
scale ever probed. Indeed, the string length, ig will be on the order of the Planck length^. 
For historical reasons, we usually work with its square, denoted a' = i^. The tension of a 
string, Ts, is set by the fundamental scale, is, as = l/(27ra'). 

We hasten to note that not every string has a length ^5! Rather, the string length Ig 
is a fundamental parameter of string theory. The point is that all string theory strings 
have an energy per unit length (tension) Tg, which encodes some characteristic length scale 

^Actually, in perturbative string theory, the Planck length is much shorter than the string length: 

(■Planck = gV^^s for string coupling < 1. 
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is, giving us an order-of-magnitude estimate of where we could detect that fundamental 
particles are actually extended objects. 

Strings come in two species, closed and open. Closed strings are loops: the string has 
no boundary. Open strings are line segments: the string has a boundary, its two endpoints. 

Strings may also be oriented or unoriented. We will discuss only oriented strings. 

As time passes, a string will sweep out a two-dimensional surface, a Riemann surface, 
called the string worldshcct. As strings have a tension, physically we expect the worldsheet 
to be minimized. The action for a single string, the Nambu-Goto action, is just the area of 
the string worldsheet, which we extremize to find the string equations of motion. 

Strings interact with coupling strength gg. Physically, Qs is the probability for a string 
to split into two strings, or for two strings to join if they cross one another. 

In perturbative string theory, we represent string interactions with diagrams similar to 
Feynman diagrams, except now instead of point particles interacting at point vertices, we 
have strings sweeping out Riemann surfaces as time passes. Such diagrams are arranged 
according to topology. The topology of a Riemann surface is determined completely by its 
genus, so we arrive at the so-called "genus expansion" of string theory. If g is the genus of the 
Riemann surface, b the number of boundaries, and % = 2 — 2g — b the Euler character, string 
diagrams are weighted by a factor gj^. When ^ 1, the genus-zero surface, topologically 
a sphere, dominates the expansion, having % = 2. If we include open strings, the first 
sub-leading correction is the disk, with ^ = 0, 6 = 1 and hence X = 1. The next correction 
would be from closed strings, the torus, with g = 1, b = and x = 0. 

We have depicted string perturbation theory in figure 2.2. We have (crudely) drawn 
the interaction of three closed strings. The quantum mechanical amplitude for such a 
process will have a leading term from a Riemann surface with the topology of a sphere 
and a subleading term from a Riemann surface with the topology of a torus. The "-|- . . ." 
indicates the infinite number of higher-genus contributions. 

At zero string coupling, for strings propagating in fiat spacetime, we can quantize the 
excitations of the string and determine the string spectrum. We find massless modes and 
massive modes. In the closed string sector, one of the massless states has spin 2, so we 
identify this state as the graviton. In the open string sector, among the massless states we 
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Figure 2.2: Cartoon of leading and sub-leading contribution to perturbative string amplitude 
for the interaction of three closed strings. 

find a spin-1 excitation, which wc identify as a gauge field. As for the massive modes, as 
the only dimensionful scale in string theory is ig, we know without doing any work that the 
masses of the excited states must be proportional to ij^. 

Suppose we consider string theory with only bosonic string excitations. Bosonic string 
theory has two fundamental problems. The first is that the spectrum includes a tachyon: 
a state whose mass-squared is negative. The second is that the theory is only consistent 
quantum mechanically in twenty-six spacetime dimensions! Why is this so? String theory 
has a local invariancc, i.e. a gauge invariance, called Weyl invariance. Weyl transformations 
are essentially (but not exactly) conformal transformations on the string worldsheet. Weyl 
invariance is generically anomalous, which is very bad^. Only in special circumstances will 
the Weyl anomaly vanish, for example in twenty-six spacetime dimensions. 

Supersymmetry, that is, adding fermionic excitations on the string, fixes some of the 
problems of the bosonic string. First, the supersymmetric string spectrum now contains no 
tachyon. Second, supersymmetric tring theory is consistently formulated in ten spacetime 
dimensions, a vast improvement over twenty-six (though still far from four). Five consistent 
superstring theories are known. Of importance to us will be the so-called type IIB theory. 
We will not need to know much about type IIB, except for some facts about D-branes which 



^An anomalous global symmetry can have important physical consequences: recall the axial symmetry of 
QCD. An anomalous gauge invariance, however, is a genuine sickness of a quantum theory. 
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we will review in the next subsection. 

Currently, string theory is only a "first quantized" theory: we only know how to quantize 
excitations on a single string. We do not yet have a full understanding of string field theory, 
where individual strings become excitations of some "string field," and where we could 
describe processes in which strings are created and destroyed. 

String theory is a theory of quantum gravity, being quantum mechanical and including 
a graviton in its spectrum, and also describes gauge interactions. String theory is thus a 
candidate for a "theory of everything." The problem string theorists would like to solve, 
then, is how string theory dynamically determines the geometry of spacetime, and all the 
low-energy physics in that spacetime. In other words, by what mechanism does string 
theory dynamically choose a vacuum state? This problem is incredibly difficult because an 
enormous number of possibilities exist. 

For example, does string theory decide that spacetime should be "compactified" to (3+1) 
dimensions, that is, to R^'^ x M, where M is some compact, very tiny six-dimensional 
manifold? If we want M = I supersymmctry in M^'^, for example, we should choose M to 
be a so-called Calabi-Yau manifold, of which about 10^'^'^ are known to exist! The problem is 
actually even worse, since string theory gives us more options. For example, we can choose 
"fluxes" on M, making the number of options even bigger than lO^'^^. 

The range of all possible vacua of string theory is called the string theory "landscape." 
The landscape is the main reason that string theory is not yet a viable theory of everything. 
For us, the key point to take away is that, since we do not yet know how string theory 
dynamically determines a ten-dimensional geometry, the current formulation of string theory 
is not "background-independent." In the current formulation of string theory, we must pick 
some ten-dimensional geometry and write an action for an individual string propagating 
in that geometry. In the back of our minds, we know the background is something like a 
"coherent state of gravitons," although no one knows a precise meaning for those words. 

Fortunately, somewhere in the landscape is AdS^ x S^. The AdS/CFT correspondence 
then gives us a less ambitious approach to the real world: rather than a theory of everything, 
string theory can be a theory of strongly-coupled gauge theories. 

Knowing that string theory is a ten-dimensional supersymmetric theory of gravity, the 
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low-energy (meaning far below the string scale, ij^), classical (meaning leading-order in Qs) 
effective description of string theory must be ten-dimensional supergravity. We understand 
supergravity fairly well, and in particular we actually know how to calculate things in 
supergravity, so it is usually the most useful tool for studying string theory, as either a 
theory of everything or in the AdS/CFT context. 

More specifically, the low-energy effective theory of type IIB string theory is called type 
IIB supergravity (surprise, surprise). We will not need to know too much about type IIB 
supergravity, except that the theory includes the metric, whose kinetic term is the usual 
Einstein-Hilbert term, and some bosonic fields called Ramond-Ramond form fields, whose 
importance becomes clear in the context of D-brane physics. 

2.3.2 D-brane Basics 

D-branes are fundamental objects in string theory. In perturbative string theory, D-branes 
appear as hypersurfaces on which open strings must end: if we are studying open strings, 
some D-branes must be around; if we introduce D-branes, they are necessarily accompanied 
by open string degrees of freedom. 

In the simplest contexts, a Dp-brane is just a (p-l-l)-dimensional generalization of a 
membrane. Like a soap bubble, or the head of a drum, a D-brane has a tension, and (part 
of) its action is simply its (p+l)-dimcnsional volume, or "worldvolumc," which it wants to 
minimize. Our objective in this section is to write this action explicitly. 

In superstring theory, D-branes carry conserved charges, and hence are stable objects. 
The charges are called Ramond-Ramond (RR) charges. In particular, a Dp-brane will act 
as a source for a RR (p-l-l)-form field. In type IIA theory, only Dp-branes with even p 
are stable, while in type IIB theory, only Dp-branes with odd p are stable. In type IIB, 
then, Dl-branes, D3-branes, D5-branes and D7-branes are stable^. A "wrong-dimension" 
Dp-brane can appear in either theory (odd in IIA or even in IIB), but will eventually decay 

D9-brane also exists in type IIB, but by itself is unstable, essentially by Gauss's law: the D9-brane 
will always fill all of space, but all of space cannot carry a net charge. We can cancel D9-brane charge in 
various ways, however, for example by introducing anti-D9-branes, or by introducing a (9-|-l)-dimensional 
orientifold, which can cancel sixteen units of D9-brane charge. In Euclidean signature, type IIB also 
includes D(-l)-branes, also known as D-instantons. 
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into stable objects: lower-dimensional D-branes and strings. 

Dp-branes with p < 3 are electrically charged under the appropriate RR field, while 
Dp-branes with p > 3 are magnetically charged. Dp-branes thus come in electric-magnetic 
dual pairs. To explain this, we consider an example. Dl-branes couple to the RR 2-form 
A^'^^ electrically, meaning with a term J A^"^^ in their action, where the integral is over 
the Dl-brane worldvolume. Let F^^^ = dA^'^^ be the 3-form field strength of A^'^\ and 
F^^) = ★F^-^) the 7-form Hodge dual of F^^\ We define a 6-form potential by F^^^ = dA^^\ 
The D5-branc then has a magnetic coupling, which is a coupling J A^^\ where the integral 
is over the D5-brane worldvolume. The magic number is six: Dp-branes and Dq-branes are 
electric-magnetic duals ii p + q = 6. D3-branes are special: they are self-dual. They couple 
to a 4-form A^^^ whose 5-form field strength is self-dual, F^^^ = dA^^^ = ~kF^^\ 

The Dp-brane action is known to leading order in gg. Stable Dp-branes are supersym- 
metric, so the Dp-brane action includes fermionic and bosonic fields. We will not need the 
fermionic part in what follows, so we will not write it. The bosonic part of the action has 
terms of two kinds, 

Sdp = Sdbi + Swz (2.3.7) 

Here Sdbi is called the Dirac-Born-Infeld part of the action, and S\vz is called the Wess- 
Zumino part of the action. Sdbi describes the Dp-brane's coupling to the closed-string 
fields of the background, namely the metric, dilaton and Neveu-Schwarz-Neveu-Schwarz 
(NS-NS) 2-form. The Wess-Zumino part of the action describes the Dp-brane's coupling to 
the RR forms, including for example the coupling J A^^^\ as well as other terms. In later 
sections, we will only need Sdbi, and only with no dilaton or NS-NS 2-form, so Sdbi is all 
we will write explicitly. 

The full bosonic action is written in the textbooks, refs. [68-70]. Explicit derivations of 
the DBI action from string theory appear in refs. [74-76]. The fermionic terms in Dp-brane 
actions have been studied in refs. [77-79]. 

To write the action, we first need a value for the tension, Tdp- In terms of fundamental 
parameters [80], 

111 . , 
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The Dp-brane action was originally computed from a calculation at leading order in the 
open string sector [74,75], which is the disk, with x = 1, hence the factor gj^ = in T^p- 
The factors of a' follow by dimensional analysis. 

Next, we need to know what fields propagate on the Dp-brane. Dp-branes support two 
kinds of fields: scalars and gauge fields. We will explain the scalars first. 

A Dp-brane will be extended in p-|-l directions, but can move in directions orthogonal to 
its worldvolume. The analogy with a drumhead makes this intuitive: a drumhead extended 
in the x-y plane can fluctuate in the z direction. Wc would thus write a "scalar field" on 
the "worldvolume" of the drumhead, z{t,x,y), to describe these fluctuations. We think of 
z{t, X, y) as saying "give me a point (x, y) on the drumhead and I will tell you where it ends 
up in the transverse direction at time t." 

A Dp-brane is just the same as a drumhead, except it is relativistic and (possibly) higher- 
dimensional. Let denote the Dp-brane's worldvolume coordinates, with /x = 0, 1, . . . ,p. 
We will always have some background geometry, with coordinates and metric guN, 
where M, N = 0, 1, ... 9. We embed the Dp-brane into the background by specifying the 
map {^^), where p of these functions will describe directions along the Dp-brane and 
9 — p of them will describe the fluctuations orthogonal to the Dp-brane. The latter are just 
like z{t, X, y) in the drumhead example. We actually have some "gauge" freedom in choosing 
the X^(^'*). One physically sensible choice is to identify the time direction of the Dp-brane 
with the time direction of the background: X^{^^^) = This is called "static gauge." 
Additionally, in what follows, we will always identify those X^{^^) that have M = 1, . . . ,p 
with the worldvolume coordinates: X^(^^) = ^'^ for M = /j, = 0,1, . . . ,p. The remaining 
9 — p functions will then be scalar fields on the Dp-brane worldvolume. 

For any gauge choice, the induced metric of the Dp-brane is 

C' = 9MN^^ (2.3.9) 

With just the scalar fluctuations, the Dp-brane action is 

Sdp = -Tdj> j cf+i^^-det {gDp), (2.3.10) 

where the integral is just the worldvolume of the Dp-brane, the natural quantity for an 
object with tension to extremize (hopefully minimize). The action is thus the energy of 
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the Dp-brane: the tension is the energy per unit world volume, and the integral is just the 
worldvolume. More specifically, if we think of the action as kinetic minus potential energy, 
then this Dp-brane action is just minus the potential energy. 

We also promised gauge fields. These arise very simply: suppose an open string has both 
ends on the same Dp-brane. The string has a tension, and (in this case) nothing forbids 
the string from minimizing its worldsheet area by collapsing to zero length. What is the 
low-energy effective description of such masslcss degrees of freedom? The answer is: a gauge 
field on the Dp-brane! In particular, with a single Dp-brane, we will have a U (1) gauge field. 
Remarkably, however, if we have Nc coincident Dp-branes, the gauge group is enhanced to 
a non-Abelian group, U{Nc). The full non-Abelian action for such a stack of Dp-branes is 
not completely known. The action for the Abelian case, either a single Dp-brane or just the 
overall C/(l) subgroup of U{Nc), is known, however. The action is of Born-Infeld type: for 
U{1) field strength F^i^, the action is 

Sdp = -Tdp j dF+^^ y - det (g^F + {2'iTa')F^,,) . (2.3.11) 

If we have two D-branes that are separated in some direction, and a string stretching 
in that direction from one D-brane to the other, then on each D-brane the endpoint of the 
string will act as a point charge, i.e. as a source for the worldvolume gauge field. With 
a stack of Nc D-branes, the endpoint of a string will act as a source in the fundamental 
representation of U{Nc) (or anti- fundamental, depending on the orientation of the string). 
The mass of the point charge will simply be the length of the string^ times the string tension. 
The length of the string is just the distance between the D-branes, I, so the mass of the 
point charge will be I /{lira'). 

Suppose we have a stack of Nc D-branes and then separate some number of them, say 
n, in an orthognal direction, so that we end up with one stack oi Nc — n, and another of 
stack of n, D-branes. We have thus introduced new open string degrees of freedom, namely 
the open strings stretched between the two stacks, which will be massive. We will have 

■*We do not mean the "string length" i^, which as we mentioned in section 2.3.1 is a fundamental scale 
of string theory. Indeed, strictly speaking, the string we are describing should have a length much greater 
than is, and hence be very heavy, such that we can treat it as a classical object, i.e. ignore quantum 
fluctuations along the string. 
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two gauge theories, one on each stack of D-branes, with gauge groups U{Nc — n) and U{n), 
as well as massive excitations in the fundamental (or anti-fundamental) representation in 
each theory. Such a procedure is the D-brane description of the Higgs mechanism. Prom 
the point of view of either stack's worldvolume gauge theory, we are studying a non-trivial 
point on the supersymmetric moduli space of the theory, where the gauge group is reduced 
and some degrees of freedom (the "W bosons" ) have become massive. 

From a field theory perspective, the U (1) part of the U {Nc) gauge group always trivially 
decouples. From the D-brane perspective, this U{1) corresponds to the overall center-of- 
mass motion of the entire stack of Nc D-branes [81]. In cases where we are only interested 
in the motion (of strings, for example) relative to the stack, we will ignore the overall U{1), 
and refer to the gauge group of the worldvolume D-brane theory as just SU{Nc). 

The Born-Infeld action was discovered in the 1930's [82]. (Schrodinger himself discovered 
that the Born-Infeld action was invariant under electric-magnetic duality [83].) At the time, 
people were worried that, according to Maxwell electrodynamics, the electric field near a 
point charge grows without bound. They did not know about QED or vacuum polarization, 
which solves this problem in the real world. Instead, they looked for a new theory of 
electrodynamics, one that preserved Lorentz and gauge invariance, that reduced to Maxwell 
theory for small electromagnetic fields, and that had some kind of maximum electric field, 
to prevent infinite fields at a point charge. They were thus led to the (3-1-1 )-dimensional 
Born-Infeld action, which obeys all of these constraints. It describes a nonlinear theory of 
interacting photons, and so is often called nonlinear electrodynamics. The original Born- 
Infeld theory did not include fluctuations of the geometry, however (the worldvolume scalars 
of a Dp-brane): the metric was simply rj^,^. 

Notice a crucial assumption: where Maxwell theory is scale-invariant, Born-Infeld theory 
must have some scale, to set the value of the maximum electric field. The scale introduced 
in Born-Infeld theory always appears in front of the F/^i, factor. In string theory, this scale 
is the inverse string tension, 2Tra'. 

In string theory, the Born-Infeld action is an effective action. Suppose we want to write 
a low-energy action for a Dp-brane. We could imagine integrating out heavy modes: an 
open string ending on the Dp-brane may have a massive excitation on it, which "tugs on" 
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the Dp-branc, exciting the worldvolume fields. At low energy, the appropriate worldvolume 
fields will be the massless excitations, namely the scalars and gauge fields. Following the 
rules of effective field theory, we would write every possible term allowed by the symme- 
tries. The Dp-brane's local symmetries include coordinate transformations on the Dp-brane 
worldvolume and gauge invariance. We must therefore include all powers of F^i^, for exam- 
ple. We should think of the DBI action as the result of summing up all the powers of F^^/, 
with the nice result that the infinite number of terms sums up to form the determinant 
and the square root. The DBI action is thus valid for field strengths F^i, < l/(27ra'): the 
electric field pulls on the ends of an open string, so if the electric field grows larger than 
the string tension, the string will be ripped apart, and our effective description will break 
down. This is the sense in which Dp-branes have a "maximum" electric field. 

In our effective action we must also include powers of derivatives of F/^^,, combined in 
gauge-invariant ways. The derivatives of the field strength must also remain below the string 
scale. If they approach the string scale, then the DBI action receives an infinite number of 
corrections from higher-dimension operators built from derivatives of i^^y. 

At low energy, Born-Infeld theory reduces to Maxwell electrodynamics. Let us illustrate 
this with the Dp-brane action. If F^i, and its derivatives are sufficiently below the string 
scale we can Taylor expand. 



Sdp = -TdpJ (F+H^-detgop 



l-{27ra'f]F^,F^''' + 0{a'^F^) 



(2.3.12) 



The first term is just the worldvolume again, and the leading term in F^i^ is just the 
Lagrangian of Maxwell electrodynamics. At low energy, then, the open string degrees of 
freedom are described just by Maxwell electrodynamics. 

A similar statement applies in the non-Abelian case, with U{1) becoming U{Nc). In 
practical terms, all we do at order F"^ is add a trace over gauge indices. We then recover 
the action of UiNc) Yang-Mills theory^, -r-^ — f trF"^, where we identify the Yang-Mills 

^As mentioned in the Introduction, we are following the conventions of D-brane physics, for which the 
definition of Qym is actually half the conventional value. The discrepancy comes from the normalization of 
the generators of the gauge group: Tr{TaTi,) = dSab- In D-brane physics, d = 1, whereas the conventional 
choice is d = 1/2. 
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coupling in terms of fundamental string parameters, 

9ym = (2W)-2 = (27r)P-2 a'^^-'^' gs (2.3.13) 

Such a simple prescription (just adding a trace) fails at higher orders in F, however, due to 
ambiguities in the ordering of gauge indices. Currently, the non-Abelian DBI action is only 
known with certainty up to ordcr^ [84]. 

Dp-branes are gravitating objects. Wc can see this easily: all Dp-branes have a tension, 
that is, energy per unit volume, and anything with energy must couple to gravity. We can 
also conclude that Dp-branes gravitate from the fact that Dp-branes interact with closed 
strings, for example by emitting a closed string, as in figure 2.3, and the fact that the 
massless mode of the closed string is the graviton. The gravitational physics of Dp-branes 
motivated the original statement of the AdS/CFT correspondence, to which we now turn. 

2.4 Statement of the Correspondence 

In this section, we will review Maldaccna's original D-brane construction and the resulting 
conjecture of the equivalence of two theories, AT = 4 SYM and string theory on AdS^ x 
[20]. Maldacena conjectured that two theories were equivalent, but he did not provide 
a dictionary between them, which we will need to compute anything. A more precise 
statement of the correspondence came shortly after Maldacena's conjecture [24,46]. We will 
review the precise statement of the correspondence in this section. 

The argument for the original AdS/CFT conjecture begins with D3-branes. The logic 
is very simple: D3-branes have two different low-energy descriptions, one in terms of open 
strings and one in terms of closed strings. 

We know the open string description: Nc D3-branes are described, at sufficiently low 
energy, by a U{Nc) Yang-Mills theory formulated on the (3-|-l)-dimensional D3-brane world- 
volume. To identify precisely which Yang-Mills theory, we simply count supercharges. Type 
IIB string theory has 32 supercharges. Introducing the D3-branes breaks half of the super- 
symmetry, hence 16 supercharges remain, the maximal number in (3+1) dimensions. The 

''Schematically, with D representing the gauge-covariant derivative, the derivative corrections are known 

up to order D'^F'^ and D^F'\ 



Figure 2.3: Cartoon of a D-brane emitting (or absorbing) a closed string. 
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theory must thus be the maximahy supersymmetric (3+l)-dimensional Yang-Mills theory: 
M = 4: SYM. Recall that H = 4, SYM has six real scalar fields. These are in fact the six 
worldvolume scalars, describing low-energy fluctuations of the D3-brane in the six directions 
orthogonal to its worldvolume. 

What about the closed-string description of D3-branes? The low-energy description of 
the closed string sector is supergravity, hence D3-branes must have a description in terms 
of the fields of supergravity. For D3-brancs, the relevant fields arc the metric and the RR 
4- form, A^'^'> . We can thus describe D3-branes as a geometry with some 4- form charge. 

So what is the geometry of Nc D3-branes? The answer must preserve the right symme- 
tries. The D3-branes break the Lorentz group of M^'^, 50(9, 1), to Lorentz transformations 
along their worldvolume, S'0(3, 1), and rotations in the remaining six directions, 50(6). 
These constraints are in fact enough to determine the answer [85]: 

dsl^ = Z{rQ)-^l^ {-de + dx^) + Zir^fl'' {drl + ds\^) , (2.4.14) 

with X = (xi, X2, xs), and where we have written the metric for the six transverse directions 
in spherical coordinates, with radial coordinate rg, the distance to the D3-branes. Here ds^g 
is the metric for a unit . To preserve the symmetries, the warp factor Z{r^ can depend 
only on the radial coordinate re, as indicated. Explicitly, the warp factor that solves the 
supergravity equations is 

^(re) = 1 + ^ (2.4.15) 

with the very important definition 

= 47r£?5-/Vc«'^ (2.4.16) 

The explicit solution for F^^^ = d^l^^) is 

= (1 + ^) a dxi A dx2 A dx3 A d{Z{rQ)-^) (2.4.17) 

with ★ the (9-|-l)-dimensional Hodge star. Notice that F^^^ is self-dual, F^^^ = ★F^^), as it 
must be. In the six transverse directions, the D3-branes look like a point charge producing 
A^'^^ flux. We can compute the flux of D3-brane charge by integrating F^^^ over the that 
surrounds the D3-branes: we find /gs i^(^) = Nc, as of course we must. 
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To understand the D3-brane geometry, eq. (2.4.14), we take limits. Suppose we are far 



then becomes that of (9+l)-dimensional Minkowski space. That makes sense: space may 
be warped near the D3-branes, but should approach flat space asymptotically far away. 

What happens as we move close to the D3-branes? That is, what if we take re ^ L? 
Now we drop the "1" in Z{rQ), and the metric becomes 



Comparing to eq. (2.2.3), we immediately identify the metric of AdS^, with radius of 
curvature L, as well as the metric of S^, with the same radius of curvature. The geometry 
near the D3-branes is thus AdS^ x S^. 

Roughly speaking, the geometry is divided into three regions: the asymptotic region, 
which is flat, the so-called "near-horizon" region, which is AdS^ x S^, and the interpolating 
region, called the "throat." The throat acts as a gravitational potential well. Suppose an 
excitation in the near-horizon region travels up the throat and escapes into the asymptotic 
region. An observer in the asymptotic region will measure an energy lower than what the 
excitation initially had, that is, a red-shift will occur. 

The question we now ask is, what low-energy physics will the observer in the asymptot- 
ically flat region see? Here again, low-energy means far below the string scale, ij^. In the 
asymptotic region, the low-energy closed-string physics is just supergravity in fiat space, 
which is a free theory in the infrared. Peering down into the throat, however, the observer 
will see all the modes of the string theory: even very high-energy string modes, if sufficiently 
far down the throat, will be very highly redshifted, hence the observer cannot discard them 
when writing a low-energy effective theory. The non-trivial low-cncrgy dynamics is thus 
that of type IIB string theory at the bottom of the throat. In other words, the low-energy 
closed string description of D3-brane physics is type IIB string theory on AdS^ x . 

We thus have two theories, M = 4, SYM and type IIB string theory on AdS^ x S^, which 
are two different descriptions of the same physics, namely the low-energy physics of a stack 
of Nc D3-branes. Maldacena then conjectured that these two theories, which a priori look 
so different, must in fact be equivalent. 



away from the D3-branes, so re L. We then drop the L'^/r^ term in Z{rQ). The metric 




(2.4.18) 



35 



A quick and easy check of Maldacena's claim (that we actually already did at the end of 
section 2.2) is to compare symmetries, which must be identical if the theories are equivalent. 
The isometry group of AdS5, SO {4, 2), is identical to the (3+l)-dimensional conformal 
group. The isometry of the internal space, the S^, is of course SO{6), which is isomorphic 
to the SU{4:)r part of the R-symmetry of A/" = 4 SYM. 

Notice that the global symmetries of the gauge theory are dual to gauge invariances 
(local transformations) of the gravity theory: the conformal group and R-symmctry group 
were dual to isomctries of the metric. Indeed, this is a generic feature of gauge-gravity 
duality, and arises because gravity admits no global symmetries^. 

Matching symmetries is not a statement about dynamics, however, and thus is far from 
a proof of Maldacena's conjecture. Nevertheless, as explained in the Introduction, a great 
deal of evidence suggests that the correspondence is true. 

An important question, though, is where exactly can we trust our solution for the metric 
and five-form? Remember, these were solutions of classical supcrgravity! We should not 
trust our solution if the curvature reaches the string scale, where our effective description in 
terms of supergravity will break down (the supergravity action will receive stringy correc- 
tions). In other words, we want the radius of curvature to be bigger than the string length, 
L^/ij » 1. We can rewrite eq. (2.4.16) as 

= AngsN,. (2.4.19) 

We understand perturbative string theory best, so we want gs <C 1. If Nc is finite, that 
would mean L'^/£g ^ 1, which is the opposite of what we want. We should thus take 
Nc ^ oo as we take — in such a fashion that AirggNc remains finite, and then chose 
AirggNc to be large, so that L^/ig ^ 1- 

Prom eq. (2.3.13), with p = 3, we can idenfity gyM = ^vr gg. What do the limits above 
translate into in the field theory? The first limit was 5s — > and Nc ^ oo with ggNc fixed. 

'^The statement that gravity admits no global symmetries comes from black hole physics. A particle 
carrying some charge can fall into a black hole. The charge must be destroyed in this process: a black 
hole has no hair, so no observer outside the hole could know that a charge fell into it. Notice we are 
not discussing electric charge, since this is not associated with a global symmetry. Electrically charged 
black hole solutions do exist in gravitational theories, and if such a black hole ate an electrically charged 
particle, then the black hole's charge would change by one unit. Arguments also exist showing that string 
theory admits no global symmetries (see for example section 4.1 of ref. [86]). 
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Defining the 't Hooft coupling A = 2 gyM^c, we can see that this Umit is the 't Hooft limit, 
9ym ~^ ^ ^^'^ Nc ^ oo with A fixed. The next limit is then to take A — > oo. The limits 



where the dual description as a classical theory of gravity is valid become the 't Hooft limit, 
with large 't Hooft coupling. 

As we have repeated a few times now, AdS/CFT gives us a new tool for studying 
strongly-coupled gauge theory: classical gravity. Maldacena conjectured that the two theo- 
ries are equivalent for all values of gs and Nc, however, which is usually called the "strong 
form" of the conjecture. The "weak form" of the conjecture, which is perhaps on more solid 
ground, equates type HB supergravity on AdS^ x with Af = 4 SYM theory in the 't 
Hooft limit, at large 't Hooft coupling. 

How do we actually use gravity to compute observables in the field theory, though? We 
need a more precise statement of the correspondence. Roughly speaking, the idea is that 
two theories being "dual" means their path integrals should be identical. The degrees of 
freedom over which we integrate will be different, but upon performing the path integral the 
result will be the same, hence the theories are physically equivalent. From this perspective, 
a duality transformation is just a change of variables, which of course cannot change the 
result for the integral. 

If we believe the strong form of the conjecture, we should imagine equating the path 
integral of A/" = 4 SYM theory with some kind of path integral for type HB string theory. 
The problem is that no one knows how to define such a path integral for string theory, 
primarily because no one knows what the appropriate degrees of freedom would be. 

If we limit ourselves to the weak form of the conjecture, we can argue that, whatever the 
string path integral should be, supergravity should be a good saddle point approximation 
to it. We can then write the precise statement of the correspondence. Let cf) be some 
supergravity field, which could be a scalar, a vector, the metric, etc. (we will ignore all 
indices). Let be a particular solution of the supergravity equations of motion, with 
leading asymptotic behavior 0q near the boundary. The AdS/CFT correspondence is the 
statement that [24, 46] 




SYM 




(2.4.20) 
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On the left hand side, we have the path integral of the = 4 SYM theory, or more 
precisely, a generating functional: O is some operator, and (f)Q is a source for that operator. 
We compute correlators of O by taking functional derivatives with respect to ^q- On the 
right-hand side is an exponential of the supergravity action, evaluated on the solution (f), 
also known as the "on-shell" action, representing a saddle-point approximation to the full 
string theory path integral. 

To compute a correlator for O in = 4 SYM, then, the AdS/CFT recipe is: 

1. ) Determine which field (f) is dual to O. 

2. ) Solve the supergravity equations for 4>. 

3. ) Plug the solution (p into the supergravity action and exponentiate. 

4. ) Take variational derivatives with respect to the leading asymptotic value (f>Q. 

Notice that we can compute the connected correlators of O by taking the logarithm 
of both sides of eq. (2.4.20) and taking functional derivatives of the on-shell supergravity 
action itself. We will define the notation for connected correlator 

iO) = = ^SsuauA (2.4.21) 

where in the first equality, the denominator is the expectation value of the identity operator 
(also known as the path integral itself), and in the second equality we have invoked the 
correspondence to translate to supergravity quantities. 

The first step in the AdS/CFT recipe can actually be very difficult. Given an operator, 
how do we determine the dual field, and vice versa? No sure-fire recipe exists for how to do 
this, so wc usually have some detective work to do. Most of the time, though, the dimension 
and symmetries of the operator are enough to identify the dual field. 

The precise statement of the correspondence itself, eq. (2.4.20), also has a big problem: 
the quantities on both sides are divergent! On the left-hand side, the generating func- 
tional has ultraviolet divergences. On the right-hand side, the on-shell supergravity action 
generically has divergences due to the infinite volume of AdS space. To make eq. (2.4.20) 
truly precise, we must introduce some consistent regularization and renormalization scheme. 
Which scheme is best often depends on the problem. 



38 



For example, suppose we have a metric that satisfies Einstein's equation and asymp- 
totically approaches AdS. If we evaluate the Einstein-Hilbert action on such a solution, we 
will find a divergence coming from integration over the infinite distance to the boundary. 
A sensible way to regulate the divergence is to introduce a cutoff: using Fefferman-Graham 
coordinates, we would integrate only to z = e rather than z = 0, with plans to send e — 
in the end. We could then renormalize by subtracting the on-shell action for pure AdS, sim- 
ilarly regulated. The two actions will have the same e ^ divergences, which will cancel 
for finite e. The e ^ limit of the difference is therefore finite. 

This method is fairly simple, but requires an explicit "background" solution to cancel 
divergences. Such a background solution may not always be available. A more general 
method does exist, however, that makes no reference to background solutions, called "holo- 
graphic renormalization" [87-92]. In holographic renormalization, we again regulate by 
integrating only to z = e. We then add counterterms on the z = e hypersurface to cancel 
those divergences coming from integration. These counterterms can preserve all symme- 
tries, such as general covariancc and gauge invariancc (when both are restricted to the 
z = e hypersurface), and make no reference to any background solution. 

An extremely important point is that the infrared (long-distance) divergences on the 
supergravity side, from integrating over the infinite distance to the AdS boundary, trans- 
late into the ultraviolet (short-distance) divergences of the field theory. We have found one 
manifestation of a more general property of the correspondence, the so-called "UV-IR re- 
lation" [93,94]. Heuristically, the UV-IR relation states that ultraviolet physics of the field 
theory is dual to physics "near the boundary," while infrared physics in the field theory is 
dual to physics deep inside AdS. In what follows, we will make no distinction between "near 
the boundary" and "the UV," or between "deep inside AdS" and "the IR." 

More precisely, the radial coordinate can in fact be identified as the renormalization 
scale of the field theory [95-97]. A good mnemonic device to recall this is to perform 
a scale transformation, which is an isometry of the AdS metric. We take — Ax^, 
with A some positive real number (not the 't Hooft coupling). For the AdS metric in eq. 
(2.4.18) to remain invariant, we must also take re re, which is the way an energy 

must transform under a scale transformation. A more rigorous derivation also exists, in 
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which the field theory renormahzation-group equations are computed from supergravity in 
AdS [96,97]. The role of the radial coordinate as the renormalization scale is then explicit. 

The UV-IR relation and the role of an extra dimension as an energy scale are believed 
to be general features of gauge-gravity duality. In other words, the way some gauge theories 
"know" about higher-dimensional gravity is through their renormalization group flow. The 
simplest flow is none at all: a conformal theory. Conformal theories are always dual to 
AdS geometries, simply because of symmetry (the conformal group). Many examples of 
non-trivial renormalization flows are also known, usually for theories that are conformal in 
the UV but flow to non-trivial IR fixed points. Holographically, such a flow appears as a 
geometry that is asymptotically AdS but deforms as we move into the bulk (a few examples 
appear in refs. [98,99]). Examples of theories exhibiting confinement in the IR are also 
known (as a start, see refs. [30, 100-102]). Of course, the dual of QCD is not known (yet)! 

We will illustrate the mechanics of the correspondence in greater detail in the next 
chapter, when we discuss probe D-branes. There we will have D-branes living in AdS^ x S^, 
with worldvolume scalars and gauge fields. We will identify the field theory operators dual to 
these fields, perform holographic renormalization for the DBI action, and explicitly compute 
one-point functions for field theory operators from the on-shell DBI action. 

2.5 The Correspondence at Finite Temperature 

The generalization of the AdS/CFT correspondence to finite temperature is fairly intuitive. 
The geometry dual to A/" = 4 SYM at finite temperature is AdS-Schwarzschild (times 
S^) [34]. We identify the Hawking temperature and the entropy of the AdS black hole with 
the temperature and entropy of the field theory [33,34]. Additionally, if we arc interested 
in thermodynamics®, we must also Wick rotate to Euclidean time and compactify the time 
direction into a circle of circumfrence l/T, as is standard for thermal field theory. 

As mentioned above, if the only scale we introduce into a conformal field theory is the 
temperature, no thermal phase transitions can occur. Once we have characterized the state 
of the theory at any one temperature, we are done with thermodynamics: we know the state 

*In the following chapters, we will also be interested in finite-temperature real-time dynamics, which 
requires Lorentzian signature. 
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of the system at any other temperature. With only one dimensionful scale, the dependence 
of thermodynamic functions, like the free energy, entropy, etc., is completely determined by 
dimensional analysis to be some power of the temperature. For example, for A/" = 4 SYM 
at large Nc but zero 't Hooft coupling, the entropy density is s = ^ [33] . 

The factor of is suggestive. Large- A?c pure Yang-Mills theory has a first-order decon- 
finement transition, characterized for example by a sudden jump in the entropy density from 
order to order N^. The heuristic way to understand this is is to think of s as counting 
degrees of freedom. At low temperatures, the appropriate degrees of freedom should be color 
singlets, i.e. glueballs, which provide the order A^^ degrees of freedom, while at sufficiently 
high temperature the appropriate degrees of freedom should be deconfined gluons, which 
provide the order degrees of freedom. If we define "deconfinement" as the leading-order 
behavior of the entropy density at large Nc (and ignore, for example, any statements about 
the potential between two test charges) , then we may declare the free, large-iVc A/" = 4 SYM 
theory to be deconfined at any temperature. While imprecise, this can be a useful way to 
think about the physics of finite-temperature = 4 SYM. 

Returning to eq. (2.4.20), and upon Wick rotation to Euclidean time, we can re-label the 
left-hand side as Z[<Pq], and call this the thermodynamic partition function in the presence 
of the source ^q. In thermodynamics, we define the free energy F by ^[^q] = e~^/-^. We 
can thus identify 

F = T SsugraI^ ^ ^o] (2-5.22) 

Up to a factor of T, the Euclidean on-shell supergravity action is the free energy of the field 
theory at large- A^c and large 't Hooft coupling. 

What about finite-density thermodynamics? In other words, what happens if we intro- 
duce a chemical potential in the field theory? First of all, we know that with two dimen- 
sionful scales, interesting phase transitions become possible. Second, we have a choice of 
ensembles, canonical or grand canonical. Ultimately, the physics does not depend on our 
choice of ensemble, but which ensemble does the on-shell supergravity action describe? For 
charged AdS black holes, the answer is known: the on-shell action is the grand canonical 
potential, also known as the Gibbs free energy [103,104]. We will denote this as Q. The 
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potential in the canonical ensemble is the Hclmholtz potential, which we have denoted as 
F. The two are related, as usual, by a Legendre transform. 

Inserting the AdS-Schwarzschild metric (times S^) into the Einstein-Hilbert action, and 
cancelling divergences by subtracting the action evaluated on AdS^ x S^, we can compute 
the free energy for A/" = 4 SYM at laxge-Nc and large 't Hooft coupling. Prom the free 
energy, we can compute the entropy density s. Surprisingly, the supergravity result for s is 
precisely 3/4 the free- field value [33]! Indeed, some have speculated that the entropy density 
smoothly interpolates from the free-field value to 3/4 the free-field value as the coupling 
increases [105,106]. 

We can actually extract an important lesson for QCD from the supergravity result for 
A/" = 4 SYM, of the "quantitative universal" variety. In QCD, asymptotic freedom suggests 
that at asymptotically high temperatures, the value of s should be that of an ideal gas of 
quarks and gluons. Lattice studies then showed that the value of s at lower temperatures, 
including temperatures explored at RHIC, was about 80% the ideal gas value [107]. Many 
people therefore hoped that perturbative calculations performed for asymptotically high 
temperatures could be extrapolated to RHIC temperatures. The lesson from J\f = 4 SYM is 
this: a value of s close to the free-field value, even 75% of it, is definitely not a sign of weak 
coupling. Indeed, the results from RHIC, for example the large "elliptic flow," indicate that 
the QCD plasma created there is very strongly coupled [15,17]. 
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Chapter 3 

ADDING FLAVOR TO ADS/CFT 

All of the fields of the A/" = 4 SYM theory transform in the adjoint representation of the 
gauge group. In QCD, the gluons are in the adjoint representation, but the quarks are in 
the fundamental representation. To make A/" = 4 SYM look a little bit more like QCD, then, 
we want to introduce additional fields that transform in the fundamental representation, 
i.e. flavor fields. In this section we will describe how this is done in AdS/CFT. Put another 
way, we will address a question of the "qualitative universal" variety: "What do flavor fields 
look like holographically?" 

IMPORTANT NOTE ABOUT UNITS: Starting now, we measure ah lengths in 
units of the AdS^ radius, L. In other words, we use units in which L = 1. An important 
conversion from supergravity quantities to field theory quantities is then a'~^ = A. 

3.1 Adding D7-Branes 

We start at the start, with a stack of Nc D3-branes in fiat space. As described in section 
2.3.2, we know that the open string description leads to a gauge theory on the D3-brane 
worldvolume, at sufficiently low energies. If D3-branes are the only D-branes present, then 
both ends of an open string must end on the D3-branes. One endpoint acts as a point 
charge in the fundamental representation while the other endpoint acts as a point charge in 
the anti-fundamental representation. The fundamental and anti-fundamental together form 
the adjoint representation. In the situation we are describing, nothing forbids the string 
from collapsing to zero length, forming a massless excitation in the adjoint representation 
of SU{Nc). We want fields in the fundamental representation, so we therefore need open 
strings with only one endpoint on the D3-branes. The other end must go somewhere, so we 
must introduce some other D-branes. 

What D-branes should we introduce, and where should we put them? As mentioned in 
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section 2.3.2, one way to produce massive excitations in the fundamental representation of 
the gauge group is to separate some number of D3-branes from the initial stack, so that we 
then have one stack of, say, n D3-branes and another oiNc—n D3-branes. We will then have 
open strings stretched from the n D3-branes to the Nc — n D3-branes, and hence massive 
fundamental-representation excitations. The physics of the worldvolume theory is a little 
different from what we want, however. This construction is interpreted in the field theory 
as moving to a non-trivial point on the moduli space of = 4 SYM, where the gauge group 
is broken from SU{Nc) to SU{Nc — n). The massive fundamental-representation fields are 
W-bosons, i.e. massive vector fields. These do have fermionic superpartners, which we 
could call quarks, but we would like to introduce flavor fields that are fermions or scalars, 
not vectors, and we would like to do so without having to move onto the moduli space of 
the theory and reducing the gauge group. How can we do this? 

Let us jump to the answer [48], and then explain why it's a good choice. We introduce 
a number Nj of D7-branes, oriented relative to the D3-branes as described by the following 
array, with the coordinates of M^'^: 

Xq Xi X2 X^ X4 X5 Xq X-j Xg Xg 

D3 X X X X (3.1.1) 

DTxxxxxxxx 

The D7-branes overlap the D3-branes in all of the D3-brane's directions, but also extend 
into four more directions orthogonal to the D3-branes. Notice that two directions, X^ and 
Xg, are orthogonal to both stacks of D-branes. 

Two things motivate this choice. The first thing is supcrsymmctry. We know that 
adding D-branes will break at least half the supersymmetry, and possibly more. We want 
to preserve as much supersymmetry as possible, to guarantee stability and retain as much 
control over the system as possible. General analyses of intersecting D-branes reveal that 
precisely half the supersymmetries will be preserved if we have four or eight directions in 
which one D-brane is extended but the other is not [68-70] . In terms of an array like the one 
above, we need four or eight directions in which only one row has an "x ." We thus introduce 
the D7-branes as above, so that in the directions X4 to X7, one D-brane is extended while 
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the other is not. 

The second thing that motivates our choice is the issue of dimensionaUty. If our only 
demand was to preserve half the supersymmetries, we could have introduced D5-branes, 
arranged as 

Xq Xi X2 X3 X4 X5 Xq X-j Xg 

D3 X X X X (3.1.2) 

D5 X X X XXX 
Such an array again has four directions with only one D-branc, X3, X4, X5, and Xq. Imagine 
an open string stretched between the D3-brancs and the D5-branes, however. The endpoint 
on the D3-branes cannot propagate throughout the entire D3-brane worldvolume, but only 
along the (2+l)-dimensional intersection with the D5-branes, i.e. in the Xq, Xi, and X2 
directions. We would like our flavor fields to propagate in the same 3+1 dimensions as the 
adjoint fields, just like the quarks in QCD. That, together with supersymmetry, uniquely 
determines the choice of D7-branes above. 

By introducing the D7-branes, we have introduced open strings whose endpoints end 
on the D3-branes, and at low energies on the D3-branc worldvolume will behave as point 
charges in the fundamental representation of SU{Nc). We can separate the two stacks 
of D-branes, so that these open strings will have nonzero length, and in the D3-brane's 
worldvolume theory the point charges will have nonzero mass. We can only separate the 
D-branes in directions orthogonal to both stacks, which in this case are the Xs and Xg 
directions. In the Xs-Xg plane, both stacks of D-branes look like points. Without loss of 
generality we can separate these points in the Xg direction. Doing so does not break any 
more supersymmetries. 

3.2 D3-Brane Worldvolume Theory 

At low energy, what is the field theory living on the D3-branes? We still have the fields of 
M = 4 SYM, but we will now also have flavor fields, preserving only M = 2 supersymme- 
try and transforming in the fundamental representation of the gauge group. The matter 
multiplet of A/" = 2 supersymmetry is called the hypermultiplet, which is comprised of two 

complex scalars and two Weyl spinors of opposite chirality, or equivalently, one Dirac spinor. 
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One scalar and one Weyl spinor transform in the A'^c of SU{Nc), and the other scalar and 
spinor transform in the conjugate representation, Nc- We will ignore this subtlety and refer 
to the fields as if they are all in the fundamental representation. In particular, we will refer 
to the hypermultiplet fermions as "quarks," in analogy with QCD. A big difference from 
QCD, however, is the presence of scalars in the fundamental representation, which, being the 
superpartners of the quarks, we will call "squarks." All of the fields in the hypermultiplet 
will of course have the same mass, m. With J\f = 2 supcrsymmetry, the Lagrangian mass 
and the physical mass are actually identical, i.e. the Lagrangian mass is not renormalizcd. 

The Lagrangian is that of = 4 SYM theory plus terms for the hypermultiplet fields. 
Again, we will not reproduce the Lagrangian here. It is written explicitly in ref. [108]. It 
will of course include the kinetic terms for the flavor fields, with gauge-covariant derivatives 
providing couplings to the gauge fields, as well as mass terms for the hypermultiplet fields, 
and quartic couplings for the hypermultiplet scalars. That is not everything, but to explain 
the remainder we need some notation. Let $ and ^ represent a complex scalar and Dirac 
fermion (respectively) of the A/" = 4 multiplet, and let (p and represent a complex scalar 
and the Dirac fermion of the hypermultiplet. Speaking very schematically, the Lagrangian 
contains Yukawa couplings of the form '^^cp and 'tp^ip., quartic couplings of the form 0^^$^$^, 
and a cubic coupling of the form m^^$^. 

What are the global symmetries of A/" = 4 SYM coupled to M = 2 hypermultiplets? 
We begin with the fermionic generators. We have already suggested that, of the sixteen 
supercharges of A/" = 4, only the eight of A/" = 2 supersymmetry remain. As for the sixteen 
superconformal generators, if m = 0, then only eight are broken, while if m is nonzero, then 
all sixteen are broken. 

What about the R-symmctry? The flavor fields break the SO{&)r part of the R- 
symmetry to 50(4)^ x SO{2)r. The 50(4)^ includes the expected SU{2)r and an SU{2) 
that rotates two of the complex scalars in the J\f = 4 multiplet as a doublet, namely the two 
scalars that do not participate in the m(j)^^(j) coupling. In the D-brane picture, the breaking 
of SO{&)r to SO{A)r X SO{2)r is easy to see from the array above. Before we introduce 
the D7-branes, the six directions orthogonal to the D3-branes have an 50(6) rotation sym- 
metry. The D7-branes clearly break this to the 5*0(4) rotations in the (A4, A5, Ag, A7) 
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directions and S0{2) rotations in the Xg-Xg plane. 

This SO{2)r ^ U{1)r symmetry acts as a chiral symmetry: the left- and right-handed 
Weyl fermions in the hyper multiplet carry opposite charges under the ^7(1)^. Of course, 
the J\f = 4 fields are also charged under this U{1)r: one complex scalar has charge +2 (by 
convention), two Weyl fermions have charge and two have charge —1. The scalars in 
the hypermultiplet are neutral. The crucial point is that this U{1)r is analogous to the 
axial symmetry of QCD, and in particular is anomalous when Nc is finite, for the same 
reasons as in QCD. (To see the anomaly in the D-brane description requires introducing the 
dilaton-axion, which we will not discuss.) Also as in QCD, a finite quark mass explicitly 
breaks the axial symmetry. Again, this is easily seen in the D-brane picture. At zero mass, 
the D3-branes and D7-branes are at the same point, right on top of each other in the Xg- 
Xg plane, so that 50(2) rotations about that point are a good symmetry (again, ignoring 
the axion-dilaton). Once we separate the D3-branes and the D7-branes, however, 50(2) 
rotations about the D3-branes are explicitly broken by the presence of the D7-branes. 

The flavor fields will not spoil the Poincare invariancc of the theory, but what about 
the conformal symmetries? A finite mass m will obviously break scale invariance, so the 
only place where we could hope to have conformal symmetry is the massless case, m = 0. 
Even with m = 0, if Nf is on the order of Nc, then conformal invariance is lost: the theory 
dynamically generates a scale. More specifically, quantum effects due to the flavor flelds 
cause the coupling to run, so that the coupling of the theory grows stronger in the UV (like 
QED, unlike QCD), suggesting that the theory requires some UV completion, that is, some 
new physics at short distances. Even at zero mass, then, the conformal group is broken to 
the Poincare group (in the quantum theory, when Nf is on the order of Nc)- 

If all A'^^ hypermultiplets have the same mass m, that is, if the mass matrix in flavor space 
is proportional to the identity, then the theory also has a global, bosonic U{Nf) symmetry, 
just like the vector symmetry of QCD (with mass-degenerate quarks). Notice that, in 
analogy with QCD, we can identify the diagonal C/(l) as baryon number. More generally, 
we could introduce a different mass for each flavor, in which case the flavor symmetry would 
be broken to U{l)'^f . For simplicity, though, we will only study the case in which all flavors 
have the same mass m. In the D-brane picture, this is the statement that in the Xg-Xg 
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plane, we will keep all Nf D7-branes right on top of each other (a single point) rather than 
separating them (many points) and placing them at different distances from the D3-branes. 

Looking ahead a bit, we know that we will eventually take Nc — oo. When we do so, 
we will always keep Nf finite, so that Nc':$> Nf, and we will never go beyond order Nf Nc 
in anything we compute (recall that generically the leading order is N^). We call this the 
"probe limit." Prom a field theory point of view, the probe limit is a simplification because 
the quantum effects of the flavor fields arc highly suppressed, and if we only work to order 
NfNc, their effects are not apparent at all. In the language of perturbation theory, this 
would be the statement that no diagrams with quark (or squark) loops contribute at the 
order to which we work. Two nice things result from this. 

The first nice thing is that, if m = 0, conformal invar iance is restored! More precisely, in 
the probe limit we neglect precisely those quantum effects, due to the flavor fields, that cause 
the coupling to run. In this limit, then, no dynamically generated scale will be apparent, 
and the theory will appear to be conformally invariant. Of course, a nonzero mass will 
explicitly break conformal invariancc, but nevertheless, the probe limit is nice because the 
only scales in any problem will be the ones we introduce by hand. 

The second nice thing that results from the probe limit is that we also neglect the 
quantum effects that produce the U{l)ji anomaly, so that, when m = 0, the U{l)ji will 
appear to be a genuine symmetry of the theory. Again, however, a nonzero mass will 
explicitly break the C/(l)j?. 

As a final note, we have been describing the low-energy physics on the D3-branes' world- 
volume, but we could just as well ask about the low-energy physics of the D7-branes' world- 
volume. The difference is just one of perspective. The low-energy theory on the D7-branes 
will be a (7-|-l)-dimensional Yang-Mills theory. What was a global symmetry on the D3- 
branes, U{Nf), is now the gauge invariance of the D7-brane worldvolume theory. The theory 
will be supersymmetric and contain adjoint fields, from open strings with both ends on the 
D7-branes, as well as fields in the fundamental representation of U{Nf), from open strings 
stretched between the D3-branes and D7-branes. These flavor fields will be restricted to 
propagate along the (3-|-l)-dimensional intersection with the D3-branes, and will have a 
U{Nc) global symmetry. 
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3.3 D7-Branes in AdS^ x 
3.3.1 The Probe Limit 

We now want to repeat Maldacena's arguments, but with the extra ingredient of D7-branes. 
That is, we want to know the low-energy closed string description of the system, i.e. the 
supergravity description. In supergravity, the D3-branes will again act as a source for the 
metric and self-dual five- form, F^^\ The D7-branes will also source the metric, and in 
addition the axion-dilaton, which is a (pseudo-)scalar field. 

The D7-branes actually complicate the picture enough that the supergravity solution 
for the D3/D7 intersection is not known for all values of Nc and Nf. Many problems arise 
when Nf is on the order of Nc, for example, the curvature diverges close to the D-branes, 
while the dilaton diverges far away from them^. 

Here is where the probe limit takes the spotlight. If we take Nc 3> Nf, so that the 
D3-branes vastly outnumber the D7-branes, then we expect the effect of the D7-branes to 
be negligible. In particular, if we have Nc oo D3-brancs and, say, just a single D7- 
brane, then we expect the D3-brancs to warp the geometry much more than the D7-branes. 
We may then study a single D7-brane "cleanly probing" the geometry produced by the 
D3-branes: the D7-branes will be embedded in the geometry but will not deform it. 

These statements can of course be made precise by examining the supergravity equations 
of motion, although we will not do so explicitly. The key point is that the D3-branes source 
the metric and five-form with a strength ggNc, while the D7-branes source the metric and 
axion-dilaton with a strength QsNf. As before, if we take gg ^ Q and Nc ^ oo such that 
gsNc is fixed, then if Nf remains finite wc will have gs^f 0, so we may safely ignore the 
axion-dilaton and the D7-branes' effect on the metric. 

The background is the one we saw before, in section 2.4, where the geometry near the 
D3-branes is AdS^ x S^, and the five-form solution is still eq. (2.4.17). The six coordinates 
transverse to the D3-branes became the radial coordinate of AdSs and the directions of the 
S^. The D7-brane is extended in four of those directions, and hence will be extended along 

^More details about the supergravity solution of the D3/D7 intersection appear in chapter 4 ref. [109], 

and references therein. 
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the radial direction and along an inside S^. We thus expect the D7-brane to be extended 
along AdS^ x in the near-D3-brane geometry. 

Lastly, taking the limit in which QgNc ~ L^/^i 3> 1, we arrive at an expanded version 
of the AdS/CFT correspondence, one that includes flavor fields: J\f = 4 SYM, at large- AT^ 
and large 't Hooft coupling, is dual to supergravity in AdS^ x S^, and a number Nf <^ Nc 
of AT = 2 hypermultiplet fields are dual to a number Nf of D7-branes extended along 
AdS5 X S^. 

3.3.2 D7-brane Embeddings 

We now want to make all of this precise, using equations. In particular, we will derive 
explicit, analytic solutions for the embedding of the D7-brane in the background geometry 
[48, 110]. Given these solutions, we will explain how to extract field theory quantities, such 
as the mass m of the hypermultiplets. 

The action that we will need for the probe D7-branes is just the DBI action. We will 
only need the Abelian version, 

Sd7 = -NfTD7 J d^C^-det{gj^J + {27ra')F^,) (3.3.3) 

Wc have Nf D7-branes, so the action is just Nf copies of the action for a single D7-branc, 
hence the overall factor of Nf. We can translate the D7-brane tension to field theory 
quantities as 

Tn7 = {27r)-'97'a'-' = ^,. (3-3.4) 

The D7-brane action is thus order \NfNc. 

To describe probe D7-brane embeddings, we will use two different coordinate systems, 
those of eq. (2.4.18) and the Feffer man- Graham coordinates of eq. (2.2.6). The latter 
are especially convenient for holographic renormalization. In this section we will study the 
worldvolume scalars of the D7-brane, not the gauge fields, which we leave for later sections. 
We will also focus only on zero temperature embeddings, leaving the finite-temperature case 
for future chapters. 

We begin by rewriting the metric of eq. (2.4.18) slightly, in a fashion that makes the 
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symmetries of the D7-brane manifest, 

dsls = Zire) {-dt^ + dx^) + Z{rQ)-^ {dr^ + ds^s + dy^ + . (3.3.5) 

The warp factor is ^(re) = r|. We have divided the six directions orthogonal to the D3- 
branes into the four in which the D7-brane is extended and the two in which neither kind of 
D-brane is extended. In the four directions parallel to the D7-branes, we have introduced 
spherical coordinates, with radial coordinate r and the metric of denoted ds'^-^. In the 
remaining two directions, which in the asymptotic geometry form the Xg-Xg plane, we have 
introduced polar coordinates, with radial coordinate y and the metric of denoted ds|i. 
Notice, then, that rg = + y^. 

The D7-brane is extended along all of the Minkowski directions, as well as r and the 
directions, and can fluctuate in directions orthogonal to itself, namely y and the direction. 
We will not consider motion in the direction. We are very interested in motion in y, 
however: we know that a nonzero separation in the asymptotic Xg-Xg plane corresponds to 
a nonzero hypermultiplet mass in the dual field theory. 

We thus consider only motion in y. As a D7-brane worldvolume scalar, y can in principle 
depend on any of the worldvolume coordinates: y{^'^)- To simplify our search for solutions, 
we will demand that the D7-brane embeddings preserve certain symmetries. For example, 
we want to preserve the Lorentz invariance of the dual field theory, so we will not allow y to 
depend on the Minkowski coordinates. For simplicity, we will also assume the embedding 
preserves the SO{4) isometry of the S^, hence y will not depend on the coordinates. The 
only remaining possibility, then, is to write an ansatz y{r). 

As an aside, wc mention that all of these assumptions may be relaxed, and indeed, more 
general embeddings contain a lot of interesting physics, about D-branes and about the field 
theory. 

With the ansatz y(r), we can use eq. (2.3.9) to write the induced metric of the D7-brane, 

dsl^ = Zire) {-dt^ + dx^) + ^(re)-^ (dr^ (1 + y'{rf) + ds%) . (3.3.6) 

The new ingredient is the y'{r) term (primes denote ^). Notice that this metric is only 
AdS^ X if y{r) = 0, the meaning of which we will discuss shortly. Taking minus the 
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determinant of the D7-brane metric and simplifying, the DBI action becomes 

Sd7 = -Nf Td7 Vm3.i Vs3 J dr ^/l + y'{rY (3.3.7) 

We have integrated over the "spectator" directions, namely the Minkowski and directions, 
producing a factor of T^s.i , the (divergent) volume of Minkowski space, as well as V53 = 2it'^, 
the volume of a unit-radius S^. In what follows, we will always divide by V^s,!, and work 
with action densities, unless stated otherwise. We thus define the notation S = S'/jy/Vjj.-i,! 
for the action density, which we will henceforth refer to as the action. We will also define a 
more compact notation for the remaining constants in front. 



N = NfTD^Vsz = J^NfNc, (3.3.8) 



not to be confused with the M of supersymmetry (the difference should always be clear from 
the context). We have also written above the translation of J\f to field theory quantities. 

The Lagrangian only depends on y'{r), so the system has a constant of motion, which 
is given by the equation of motion for y{r), 

where the quantity in parentheses is "conserved." Clearly, setting y'{r) = is a solution, 
that is, a constant y{r) is a solution. We could also consider a more general solution, by 
setting the quantity in parentheses equal to a constant, solving for y'{r), and integrating to 
find y{r). These solutions actually have larger action than the y'{r) = solutions, as we 
will show in Chapter 4, so for now we focus on the y'{r) = solutions. 

We have found that y{r) = c, for any constant c, describes the embedding of a D7-brane. 
We would like to interpret this solution, on both the supergravity side and the field theory 
side of the correspondence. 

To understand what y{r) = c means in the supergravity description, we return to the 
induced D7-brane metric, eq. (3.3.6), setting y'{r) = and y(r) = c. Notice that for the 
D7-brane, this means r| = + c^, so that the metric becomes 

dsm = + c^) {-dt^ + dx'^) + {dr^ + ds^) . (3.3.10) 
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Figure 3.1: Cartoon of D7-brane embedded in AdS space. 

As we decrease r, from r = oo down to r = 0, we find that rg decreases from rg = oo 
to re = c, so for nonzero c the D7-brane only reaches some minimum value of vq. As we 
decrease r, from the metric factor in front of ds^a, we can also see that the volume of the 
S'^ dcscrcascs. Indeed, at r = 0, the volume of the S"^ is zero! Recalling that rg is the 
radial coordinate of AdS^, we have found that, for finite c, the D7-brane extends from the 
boundary at rg = oo down into the bulk of AdS^, with the shrinking as it goes, until at 
rg = c the has collapsed completely^ and the D7-brane simply terminates, extending no 
further in rg. 

The embedding of the D7-brane is depicted (crudely) in figure 3.1, which is a side-on 
view of AdS as depicted in figure 2.1. The radial direction rg runs vertically, while one of 
the Minkowski spatial directions runs horizontally. The boundary at rg = oo is at the top 
while the Poincarc horizon rg = is at the bottom. The D7-branc extends throughout the 
region of dashed lines, but ends at some rg = c, and does not extend to smaller values of rg. 

What is the field theory interpretation of y{r) = c? These solutions actually describe 
massive hypermultiplet fields. To see this, recall that if the D3-branes and D7-branes are 
separated in the Xg-Xg plane (in the asymptotic geometry) then the stretched open strings 
between them represent massive flavor fields in the D3-brane theory. Recall also that, in 
the asymptotic geometry, y is precisely the separation in the Xg-Xg plane. Given y(r), we 
find the asymptotic separation by taking by lim^_>oo y{r), which is of course c. The mass of 



■^The 5^ is a topologically trivial cycle of the S^, which makes this possible. 
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the flavor fields is just the string tension times this separation: m = 

We have thus discovered what massive flavor fields look like holographically: D7-branes 
that "end" somewhere in AdS^. Such a description has a nice interpretation in the field 
theory. Far enough in the IR, we expect any massive fields to be absent from the dynamics. 
Indeed, the D7-brane is absent for re < c. Conversely, far in the UV we expect a finite mass 
to be unimportant. As r — *■ oo, we can drop the (? part of r| = + c^, so that the induced 
D7-brane metric becomes pure AdS's x S*'^. The AdS^ factor indicates that in the UV the 
conformal symmetry is restored: the finite mass is not apparent. 

If m precisely equals zero, then re = r and the D7-brane extends all the way to = 0, so 
that the D7-brane geometry is exactly AdS^ x S^. In the field theory, conformal symmetry 
is completely restored, and not just approximate in the UV. 

We end this section by re- writing the above action and solutions, making two changes: 
first, we will use Fefferman-Graham coordinates, and second, we will use a different param- 
eterization of the worldvolume scalar. We rewrite the AdS^ x metric as 

dsjjs = ^ {dz'^ - dt'^ + dx^) + de^ + sin^ edsli + cos'^edsls (3.3.11) 

These coordinates are related to the previous ones as: z = l/r^, r = cos 6, and y = tq sin 9. 
The worldvolume scalar is now 6. Using the same symmetry arguments as before, we 
consider an ansatz in which 6 depends only on the radial coordinate, z. In figure 3.2 we 
have depicted the as a sphere, and the that the D7-brane wraps as a circle, with 
the angle indicated. 9 runs from zero, where the D7-brane wraps the maximum-volume 
equatorial 5^, to 7r/2, where the has collapsed to zero volume (at the "North pole"). 
The induced D7-brane metric is now 

dsjj^ = + G'i^)^ + ^ {-dt^ + dx^) + cos^ 9{z) ds% (3.3.12) 

where the prime on 9'{z) denotes d/dz. The D7-brane action becomes 

S = -M j dz^ cos^9{z)^l + z'^ 9'{zY (3.3.13) 

The equation of motion for 9{z) is not particularly illuminating, so we will not reproduce it. 
Fortunately, however, we already know that y{r) = c is a solution, so we can just translate 
that solution into these coordinates. Using y = re sin ^, we have 9{z) = arcsin(cz). 
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Figure 3.2: Cartoon of S^, with D7-brane wrapping S^, and the definition of the angle 6. 

We can translate what wc know about these solutions into our new notation. The 
boundary is now at z = 0, where the arcsin is zero, hence 9{z) = 0. Looking at the D7- 
brane metric, we can see that the has maximum volume (cos 6{z) = 1) there. As we move 
into the bulk, increasing z, we see that 6{z) increases, until we reach 2; = 1/c, at which point 
9{z) reaches its maximum, 6{z) = 7r/2, where the volume of the is zero (cos^(2;) = 0). 
The asymptotic form of the solution is 9{z) = cz+0{z^), so in these coordinates we identify 
the mass from the coefficient of the term linear in z. 

3.4 Holographic Renormalization 

3.4-1 Regulator and Counterterms 

Recalling the precise statement of the AdS/CFT correspondence, eq. (2.4.20), we know 
that eventually we will need to compute the on-shell D7-brane action and take functional 

derivatives with respect to the leading asymptotic values of the D7-brane worldvolumc fields. 
Let us try to do so naively. Inserting the solution y{r) = c into the action eq. (3.3.7), we 



where we have now explicitly written the endpoints of the r integration. The integral is 
plainly divergent! In this section we will explain how to regulate and renormalize the D7- 



find 




(3.4.14) 
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brane action, using the method of holographic renormalization [87-92], as first done in 
ref. [49]. 

We will work in Feffer man- Graham coordinates throughout this section. We will thus 
use eq. (3.3.13) as the form of the D7-brane action. The main reason for doing so is that 
the scalar 6{z) explicitly depends on the AdS^ radial coordinate, and hence looks like a 
scalar field in AdS^, as opposed to y, which depends on r rather than vq (or z = l/re). 

Holographic renormalization proceeds in two steps (the same steps as in any textbook 
treatment of renormalization): introduce a regulator, and then introduce countcrtcrms to 
cancel divergences as the regulator is removed. The procedure is best illustrated by example, 
so we will describe it as we go. 

The action diverges because of the integration over the radial direction: AdS^ has infinite 
volume, heuristically because "the boundary is infinitely far away." We thus introduce a 
regulator simply by integrating, not all the way to the boundary at z = 0, but only to some 
z = e, with the intention of sending e — in the end. We thus write a regulated action, 

/oo 1 
dz ^ cos^ 9{z) + e'{zy (3.4.15) 

We next need to determine the divergences. To do so, we need the asymptotic form of 6{z) 
(the form near the boundary). We can use 9{zys equation of motion to show that it must 
have the asymptotic form 

9{z) = z ( 0(0) + 0(2) + 0{z^) ) . (3.4.16) 

In the 6{z) = arcsin(c2:) example, we have 0(0) = c and 0(2) = g c^, but in fact the above 
asymptotic form will be the same in any background geometry that is asymptotically AdS^ , 
as we will discuss later (around cq. (3.4.27)). We will therefore leave 0(o) and 0(2) as generic 
asymptotic coefficients, only plugging in 0(o) = c and 0(2) = g in explicit examples. 

The equation of motion actually determines all the higher-order coefficients in terms of 
these two, that is, all the higher order coefficients are just some functions of 0(o) and 0(2) 
that we can compute by demanding that the equation of motion vanish order-by-order in 
powers of z. The equation of motion does not determine 0(o) or 0(2), however. We know 
that the first is (related to) the mass rn of the hypermultiplets, which we should be free to 



56 



vary. What determines ^(2)? The equation of motion is (in this case) just some ordinary 
differential equation. To solve it, we must specify boundary conditions. One boundary 
condition we specify is the leading asymptotic value, that is, the value of O^oy We then 
imagine "integrating in" the solution into the bulk. Somewhere in the bulk (the IR), we 
require another boundary condition, such as smoothness of the solution 0{z). The entire 
solution is thus specified, and in particular, the value of ^(2) is fixed. In short, ^(2) is fixed 
by ^?(o)5 but '^^ '^ili i^ot see that in the asymptotic expansion of the equation of motion, since 
the way it is fixed requires information about the boundary conditions in the bulk of AdS^. 
These statements are in fact true generally for fields in AdS^: in the asymptotic expansion, 
the leading and sub-leading coefficients determine all higher coefficients, and the sub-leading 
term is fixed by the leading term via some boundary condition in the bulk [95, 111]. 

In a little while, we will provide a field theory interpretation of these statements. Before 
we can do so, we need to know how to translate 6(2) into field theory quantities. We will do 
that in the next section. 

Wc can now plug the asymptotic solution into the regulated action, Taylor expand in z, 
and integrate to determine the e — > divergences, 

Sre, = -^J/^{^ - ^(0)^ + 0(.)) = -M (—.-Idfo)^ + (3.4.17) 

Notice that we have not specified the upper endpoint of integration, since that may vary 
from solution to solution. For example, in the 9{z) = arcsin(c2;) solution, we should only 
integrate to z = 1/c. 

We now want to introduce some kind of counterterms to cancel the e ^ divergences. 
These counterterms will be localized on the hypersurface z = e. To understand why, we 
invoke the UV/IR relation: we know that these IR divergences in supergravity translate into 
UV divergences in the field theory. In the field theory, the cancellation of UV divergences 
does not involve IR physics. On the supergravity side, then, we should be able to cancel 
the divergences just with counterterms^ "near the boundary," at 2; = e. 

In the field theory we know that the regulator and counterterms will preserve all global 

^On the gravity side, the counterterms are often also necessary for variational problems to be well-posed. 
See for example ref. [112]. 
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symmetries, except for anomalous symmetries. On the supcrgravity side, we thus expect the 
counterterms to respect all the gauge invariances of the theory, and in particular, coordinate 
transformations, or at least coordinate transformations on the z = e hypersurface. 

We define the induced metric on the regulator hypersurface as ^ij. Written explicitly 
for AdS^ in Fefferman-Graham coordinates, ^ij is 

dsl = -fijdx'dx^ = ^ {-dt'^ + dx'^) (3.4.18) 

We define 7 = det 7^^ , which for AdS^ is simply 7 = — 1 /e^ . We can now write the coun- 
terterms we need to cancel the e — > divergences of the D7-brane action [49]: 

Li = +^^f^/^ L2 = ~^f^/^e{ef (3.4.19) 

These have an e ^ expansion 

Li = +-^1^^ L2 = -^flefo) ^ - ^(^(0) 0(2) + 0{e') (3.4.20) 

which, by construction, cancel the divergences in eq. (3.4.17). 

Notice that the leading divergence, and the first counterterm, Li, are independent of 
0{z). The leading divergence is just the divergence from the infinite volume of AdS^, and 
will be present for any solution 0{z). Also notice that the second counterterm, L2, includes 
a finite contribution of the form ^(g) ^(2)) which will be important shortly. 

We construct the renormalized action by summing the regulated action and countert- 
erms, and then removing the regulator, that is, taking e ^ 0. We will call the sum of the 
regulated action and counterterms the "subtracted action," 

Ssub = Sreg + ^ -^i (3.4.21) 

i 

because we have subtracted the terms that diverge as e — > 0, and the "renormalized action" 
as the e — > limit of the subtracted action, 

Sren = lim Sguh- (3.4.22) 

e— >0 

Actually, one ambiguity remains: we can introduce a finite counterterm: 

Lf = CN^e{ef = CMet. + 0(e2) (3.4.23) 
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where C is, at the moment, an undetermined coefficient. How can we fix C? 

We can determine C as follows. We know that the solution 0{z) = arcsin(cz) is supersym- 
metric^. Supersymmetry demands that the action evaluated on the solution, i.e. the energy 
of the D7-brane, be zero. We will adjust C to make this so. Plugging 6{z) = arcsin(c2;) into 
the action, we find 

SregL, = -M j'" dz ^ {l - ) 

Notice that the upper endpoint of the z integration is 1/c. Recalling that ^(q) = c, we see 
that the divergences have the expected form shown in eq. (3.4.17). Adding Li and L2 and 

taking e ^ 0, we find that the renormalized action is 

Sren = -M^c'' (3.4.25) 

which is clearly nonzero. To get zero, we fix the coefficient C of the finite counterterm: 

Lf = +^AfV^eie)' (3.4.26) 

Starting now, the renormalized action will always include Lf as well as Li and L2. 

The method we have outlined is actually valid for a wide class of asymptotically AdS 
geometries. To be specific, consider a metric of the form 

dz'^ 1 

ds^ = ^ + -^gii(x\z)dx'dx^. (3.4.27) 
From Einstein's equation, we can show that gij{x^, z) has an asymptotic expansion 

g.,{x\z) = glf{x^) + gif{x^)z' + 0{z'), (3.4.28) 

For the moment, let us assume that gf^^ is Ricci-flat. In such straightforward 
exercise reveals that in fact g^^^ = 0, and the first sub-leading term enters at order z'^ [89]. 

Any such metric has the same near-boundary behavior as pure AdS^, that is, the spatial 
slices diverge as when z ^ hence anything that depended only on the asymptotic 

^Recall that this solution describes hypermultiplets in the field theory with an A/' = 2 supersymmetry- 
preserving mass m. 
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behavior in z will be unchanged. In particular, the asymptotic form of ^(-z), eq. (3.4.16), 
will be unchanged, as will the divergences of the on-shell D7-brane action. Some of the 
intermediate equations may look different, for example, factors of \J — det g'(^) will appear, 
and 7jj will have a different form, but the counterterms are written covariantly, so their 
form is unchanged. As long as g^^^ is Ricci-flat, Li, L2 and Lj will suffice. 

As an example, the AdS-Schwarzschild metric in eq. (2.2.6) has the form of eq. (3.4.27) 
with Ricci-flat so the counterterms Li, Li and will be sufficient. This agrees nicely 
with our field theory expectations. A general statement in thermal field theory is that 
introducing a finite temperature does not introduce new UV divergences. In other words, 
in finite-temperature field theory, the zero-temperature counterterms are sufficient. 

If g^^ is not Ricci-flat, additional counterterms are required. To be covariant, these 
counterterms must be built from the Ricci tensor (with indices properly contracted) and 
Ricci scalar of 7ij. Introducing a g^^ with nonzero curvature corresponds to studying a 
field theory on a curved space. Many interesting physical questions arise when we study 
A/" = 4 SYM on a curved manifold, simply because the manifold introduces a new scale into 
the theory: the curvature. We will only be interested in field theories in flat space, though, 
so we will not need the additional counterterms. They are written explicitly in ref. [49]. 

L\ and Li are actually not unique to the D7-brane. L\ is present merely to renormalize 
the infinite volume of ^^55, and hence will be present in any action that involves integrating 
over AdSs, or at least over the region near the boundary. Li is present for a free scalar in 
AdS^^ with a particular mass. Both of these are not surprising: the D7-brane fills AdS^ and 
possesses a worldvolume scalar. Indeed, suppose we Taylor expand the action in powers of 
Q{z)- The result will have the schematic form 

s ~ ^dz^g[\^)^ {g^'di.ed.e + M^e^) + 0(^(2)^)^ (3.4.29) 

with g/j^i, the metric of AdS^ and g its determinant. Clearly, the first term is just the volume 
of AdS^ while the quadratic term is that of a free scalar with mass M^. The remainder 
looks like an infinite series of interaction terms. 

The mass-squared of the worldvolume scalar turns out to be = — 3 (in units of the 

AdS^ radius), which looks worrisome at first, being negative. Fortunately, in AdS space. 
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scalar fields can have a slightly negative and still be stable: due to the curvature of the 
space, the total energy remains finite as long as > —4. This is called the Breitenlohner- 
Preedman bound [113, 114]. The worldvolume scalar's mass is clearly in the safe region. 

Only Lf is really special to the D7-brane. In fact, we will see in the next section that 
L f must be present to extract the correct answer for the one-point function of the operator 
dual to 0{z). 

3.4-2 The Dual Operator and Its One-Point Function 

In this section we identify the field theory operator dual to 6{z) and use holographic renor- 
malization to compute its expectation value (one-point function) from the renormalized 
on-shell D7-brane action. 

What field theory operator is dual to the D7-branc worldvolume scalar? Wc know 
that the leading asymptotic value of the worldvolume scalar is, up to some constants, the 
hypermultiplet mass m. Recalling the precise definition of the AdS/CFT correspondence, 
eq. (2.4.20), we know that the leading value of the bulk field acts as a source for the dual 
operator. The dual operator must therefore be given by taking d/dm of the field theory 
Lagrangian, and must therefore include at least the mass operator of the fermions in the 
hypermultiplet and m times the mass operator for the scalar fields in the hypermultiplet. 
Notice the operator will thus explicitly depend on m. We will denote the operator as Om- 

A precise derivation of Om was performed in rcf. [35]. Due to supcrsymmetry, Om 
actually includes more than just the mass operators of the hypermultiplet fields. It also 
includes an interaction term with an adjoint field, from the m(j^<^(j) coupling mentioned 
in section 3.2. We will just write the operator schematically. Let represent the two 
complex scalars in the hypermultiplet (a = 1,2) and represent the Dirac fermion in the 
hypermultiplet, and let $ represent one of the complex scalars in the A/" = 4 multiplet. Om 
has the form 

Om = i^^ + {^^i^a + V2ct>\^ct>c) + h.c. (3.4.30) 

a 

For most practical purposes, however, we can think of Om as being just the hypermultiplet 
mass operator. 
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Wc should be able to compute correlation functions of Om from the on-shcll D7-brane 
action. We will illustrate the procedure by computing the one-point function for Om- 
According to the precise statement of the correspondence, the expectation value of Om is 
given by taking the derivative of the on-shell D7-brane action with respect to the leading 
asymptotic value of 0{z): 

(C„) = ^ = (W)|^^ (3.4.31) 

After introducing the regulator, {Om) can actually be re-written as 

(0„.)=(W)lhn^^^. (3,4.32) 

Essentially, all we have done here is introduce some factors of e so that we can write (Om) as 
a variational derivative with respect to 9{e), which is easier to compute than the derivative 
with respect to 0(o) • The contribution from the regulated action is 



+Ar^ cos^e{e) , ^^^'^^^ . (3.4.33) 



59{€) 59'{z) 

where we have defined the Lagrangian via Sreg = J^d^ L and in the first equality we have 
used the equation of motion, 

5L d 5L , „ ,x 

(3.4.34) 



seiz) dz 8e'{z) 

and then performed the z integration. Using the asymptotic form of 9{z), we find 

? 7^ ^ = -^1? + - ^Ofo,) + 0(.-) (3.4.35) 

The contribution from the counterterms is simple. Li is independent of 9{z) and hence 
contributes nothing. For the remaining counterterms we find 

= -MV^9{e), ^ = +lMV^9{ef (3.4.36) 



S9{e) " ' ' 59{e) 3 

so that these contribute 



= -M^ -N9,2)+0{ef (3.4.37a) 



-Af9t. + 0{e^) (3.4.37b) 
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Summing everything and taking e — > 0, we find 

{O^) = N{2T,a') (29^2) - ^^fo)) (3-4.38) 

Given a solution 9{z), we can compute the associated value of {Om) by extracting the 
asymptotic coefficients 0(o) ^'^d 0(2), and combining them as indicated in eq. (3.4.38). We 
now also understand the meaning of 9(2) i^i the field theory: its value determines (Cm)- 

In fact, this is a generic feature in gauge-gravity duality: given some solution of some 
field, in the asymptotic expansion the coefficient of the leading term acts as a source for 
the dual operator while the coefficient of the sub-leading term determines the expectation 
value of the dual operator [95, 111]. The leading term, being the lower power in z, will be 
non-normalizable, meaning that if we insert just that term into the action and attempt to 
integrate all the way to the boundary, we find a divergence. The sub-leading term will be 
normalizable, meaning that in the same process the result will be finite. 

We can now return to a point we mentioned above. Prom the supcrgravity point of view, 
we have a field with some equation of motion, in this case 9{z), whose equation of motion 
is just an ordinary differential equation in z. We specify a boundary condition at the AdS 
boundary, in our case the value of 9(f)), and integrate the solution into the bulk, where we 
impose some other condition, for example regularity of the solution. The rest of the solution 
is then determined. In particular, our 9(2) will be determined. 

We are now ready to interpret these statements in the field theory. Choosing the value 
of 9(Q) means choosing m, a parameter in the Lagrangian. Prom a field theory point of 
view, we then expect the dynamics of the theory, including IR physics, to determine (Cm). 
That is what we sec happening in the supcrgravity description when we see the boundary 
condition in the bulk (the IR) fixing the value of 0(2). 

We end with an example: what value of (Cm) does the solution 9{z) = arcsin(cz) 
produce? Recalling that ^(q) = c and 0(^) = \(?, we find precisely {Om) = 0. Indeed, 
supersymmetry demands that (Om) = [50]. A fermion bilinear, such as the mass operator 
of the hypermultiplet fermions, may be written as the supersymmetric variation of another 
operator, which must vanish. Notice that the finite counterterm Lf, which was required by 
supersymmetry, was essential to get (Om) = 0- 
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If we break supersymmetry, however, we may find a nonzero value for (Om)- We can 
break supersymmetry by introducing a finite temperature, density, or electric and magnetic 
fields. The effect of doing so is the subject of the next two chapters, although most of our 
attention will be on the on-shell value of the D7-brane action, rather than the value of (Om)- 
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Chapter 4 

THERMODYNAMICS OF FLAVOR FIELDS 

We are now ready to study thermodynamics of flavor fields using the holographic de- 
scription. We will add one ingredient at a time: temperature, chemical potential, electric 
fields and magnetic fields. The emphasis will be on the physics and the qualitative form of 
the phase diagram. We will only present explicit results in two cases, both at finite chemical 
potential: finite temperature and zero mass (T > 0, m = 0), and zero temperature and finite 
mass (T = 0,m > 0). In each case, we will present exact solutions for D7-brane worldvol- 
ume fields, and using these and the on-shell D7-brane action, draw conclusions about the 
thermodynamics of the field theory. In the former case, we will find no phase transitions, 
while in the latter case we will find a second-order transition. We will describe how these 
results fit into the "big picture" of the full phase diagram as we go. 

4.1 Finite Temperature 

4.1.1 Thermal Phase Transition 

In A/" = 4 SYM, in the 't Hooft limit and with large 't Hooft coupling, a state in ther- 
modynamic equilibrium at temperature T is dual holographically to supergravity on the 
(4-|-l)-dimensional AdS-Schwarzschild geometry times S^. The temperature and entropy of 
the field theory are identified with the temperature and entropy of the AdS-Schwarzschild 
black hole. The AdS-Schwarzschild metric is written explicitly in eqs. (2.2.4) and (2.2.6). 
We will work in Euclidean signature, with compact time direction, until stated otherwise. 

We introduce flavor fields as before, with a holographic description as probe D7-branes. 
The on-shell action of the D7-branes now corresponds (up to a factor of T) to their contri- 
bution to the free energy, which will be order XNfNc because their action is order XNfNc. 

In the AdS-Schwarzschild background, two types of D7-brane embedding are possible 
[50-53]. As before, we will have D7-brane embeddings in which the D7-brane ends at 
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some radial position, where the collapses to a zero volume, outside of the black hole 
horizon. Starting now, we will follow convention and call these "Minkowski embeddings." 
The presence of the horizon admits a second class of embeddings, however, the so-called 
"black hole embeddings." These are D7-branes for which the shrinks but never fully 
collapses as the D7-brane extends into the bulk: the D7-brane exends all the way to, and 
intersects, the black hole horizon. Intuitively, we can think of the D7-brane as extending 
into, and ending somewhere inside, the black hole. The black hole interior is not part 
of the geometry, however, so we should not take this too seriously. An observer in AdS- 
Schwarzschild will simply see the D7-brane intersect the horizon. 

The two classes of embeddings have distinct topology. For a Minkowski embedding, the 

C collapses to zero volume, while for a black hole embedding it does not. Instead, for 
a black hole embedding, at the horizon the gu component of the metric vanishes, indicating 
that the Euclidean time circle collapses to zero size there. 

Invoking the UV-IR relation and recalling that the radial direction acts as an energy scale 
in the field theory, the position of the horizon corresponds to the scale T, while the position 
where the D7-brane ends corresponds to the scale ni. Roughly speaking, the Minkowski 
embeddings correspond to large m, compared to T, while the black hole embeddings corre- 
spond to small m. For example, the zero-mass solution, y(r) = 0, will obviously intersect 
the horizon: the sits at the equator of the for all values of r. 

In the field theory, we now have two scales, T and m, so we have the possibility of a 
phase transition. Suppose we fix m to be some large value and then change the temperature. 
We must measure T in units of m, so we imagine increasing T/m and asking whether a 
phase transition occurs. 

The supergravity analysis reveals that, indeed, a phase transition does occur [50-56]. In 
the holographic picture, we are starting with a Minkowski D7-brane and increasing T, so 
the horizon begins moving up towards the D7-brane. Eventually, the horizon will reach the 
D7-brane's endpoint. Following convention, we will call this the "critical embedding" of the 
D7-brane. As we continue increasing T, the D7-brane will intersect the horizon, becoming a 
black hole embedding. The topology of the D7-brane changes in this process, so we expect 
some kind of discontinuous transition. 
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By computing the value of the on-sheU D7-brane action for these embeddings, a tran- 
sition has been discovered: the minimum of the D7-brane's action does jump from a 
Minkowski embedding to a black hole embedding. In fact, this occurs before the hori- 
zon reaches the Minkowski D7-brane's endpoint. In other words, the D7-brane "jumps 
into" the black hole before reaching the critical embedding. 

To show this explicitly, we show in figure 4.1 the results of the numerical supergravity 
calculation. To perform the calculation, wc use Fefferman-Graham coordinates and the 
worldvolume scalar field 0{z), with the counterterms for the action written in section 3.4. 
In the plots, the vertical axis is the value of the renormalized on-shell D7-brane action, Sren, 
divided by the normalization factor Af. On the horizontal axis is the leading asymptotic 
coefficient, 0(o), in units of the AdS^ radius. The black curves come from black hole em- 
beddings while the gray curves come from Minkowski embeddings. The critical embedding 
is thus where the black and gray curves meet. 

Recall that the hypcrmultiplet mass is given by m = and conversely 0(o) = 27r\^m, 
again in units of the AdS^ radius. Our calculations are done with T = 1 in units of the 
AdS^ radius. We may thus write 0(o) = 27r\/Am/T. 

In figure 4.1 (a.) we show the large-scale structure. Notice in particular that as 0(o) 
grows, the action approaches zero. This conforms to our field theory expectation that as 
m ^ oo, the fiavor fields' contribution to the free energy should vanish, as they decouple 
from the dynamics. Notice that, on these scales, the action appears to be single- valued. 

In figure 4.1 (b.) we show a close-up near the critical embedding, which reveals that 
the action is not single- valued. The ground state will be the state with the smallest value 
of the on-shcll action (in field theory language: the lowest free energy). Suppose we fix 
m and increase T, which appears in this plot as moving from the right to the left, since 
0(0) ~ m/T. As we enter figure 4.1 (b.) from the right, the smallest values are on the gray 
curve, i.e. the Minkowski embeddings. We can see that two other branches of solutions 
also appear, corresponding to black hole embeddings, with larger action. Once we reach 
the critical solution at 0(o) ~ 2.90, however, the minimum values of the action are now on 
the black curve, i.e. the black hole embeddings. As we continue moving to the left, the 
minimum action solutions remain the black hole embeddings. The system never accesses 
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the critical embedding, at 9(0) ~ 2.88. The action/free energy is thus continuous, but has a 
discontinuous first derivative, a "kink," at the phase transition point, 9(o) ^ 2.90. 

Translating to field theory language, we conclude that the flavor fields experience a first- 
order phase transition at T 2.90 / 2'K\f\m. We then expect any observable associated 
with the flavor flelds to exhibit a discontinuous change at that point. A simple observable 
to compute on the supergravity side is (O^), which is the flrst derivative of the free energy 
with respect to m, and so of course displays a discontinuity at the transition point [50-56]. 

A more detailed analysis of near-critical D7-brane solutions revealed that, in fact, the 
solutions exhibit a self-similar structure. Such self-similarity of the supergravity solutions 
naturally suggests that the phase transition in the field theory will be first order [53,56]. 

Meson "Melting" 

The simplest way to characterize the phase transition is in terms of the meson spectrum in 
the field theory. We have not yet described the meson spectrum in the field theory, or how 
to compute it from supergravity, so let us briefly review the "qualitative universal" question, 
"How do we compute a meson spectrum holographically?" We begin at zero temperature. 

A/" = 4 SYM plus M = 1 hypermultiplets, in the 't Hooft limit, with large 't Hooft 
coupling and iV/ ^ iVc, is not a confining theory^, and will not have mesons in the same 



^The potential between quarks (and/or squarks) is known in some limits, see for example ref. [110]. 
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sense as QCD, i.e. a quark and anti-quark connected by a flux tube. The theory may have 
quark-anti-quark bound states, however. At zero temperature, we know from dimensional 
analysis that the mass of any such meson must be proportional to m. 

Holographically, a meson is represented by an open string with both ends attached to 
the D7-brane: one endpoint represents the quark and the other endpoint represents the 
anti-quark. In flat space (not AdS^ x S^), we know that such strings will have massless 
excitations, in which the string collapses to zero length, and massive excitations, with 
masses on the order of £j^. The former give rise to the massless world volume fields of the 
D7-brane, which appear in the DBI action. Now embedding the D7-brane and open strings 
in AdS^ X S^, in the limit where the string length is much smaller than the AdS^ radius, 
we expect the masses of all states to acquire a correction on the order of the inverse AdS^ 
radius. The massless excitations thus acquire masses of order the inverse AdS^ radius. The 
spectrum of excited string states will be largely unchanged, but will be far heavier, by at 
least a factor of £j^, than the lightest states. In the field theory, we want to know the 
lightest degrees of freedom, as these dominate the low-energy physics, so we will focus only 
on the lightest string excitations. 

The lightest states excite the D7-brane worldvolume fields, which we know how to de- 
scribe using the DBI action. Our question thus becomes: what is the spectrum of excitations 
of the D7-brane worldvolume fields? In practical terms, we first pick a particular D7-brane 
embedding and then compute the action for fiuctuations about that solution. As we only 
want the spectrum, we only need the action for such fiuctuations to quadratic order. We 
will then need to solve the linearized equation of motion of the fiuctuations, which in general 
is a complicated partial difi"crcntial equation. 

Remarkably, this calculation has been done analytically for the probe D7-branc in AdS^ x 
S^, in ref. [110]. In other words, at zero temperature, exact solutions for the fiuctuations 
of the D7-brane fields have been computed, and hence the exact meson spectrum is known. 
We will skip the details. The three salient features are: the spectrum is gapped, discrete, 
and the meson masses scale as m/y/X. The 't Hooft coupling, A, is very big, so these mesons 
are actually much lighter than the fundamental quarks and squarks. Put another way, these 
mesons are very deeply bound. 
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That the spectrum is discrete is easy to understand just from basic facts about differential 
equations. If m is nonzero, then the D7-brane ends at some value of the radial coordinate. A 
fluctuation will then have either Neumann or Dirichlet boundary conditions at the endpoint. 
If we imagine the fluctuation as a wave traveling down the D7-brane, then at the endpoint 
the wave will be reflected, giving rise to a discrete set of eigenfrequencies. 

A special case is zero mass. With no endpoint to "bounce off of," the fluctuation will 
simply continue traveling all the way to the Poincare horizon. The spectrum of excitations, 
and hence the spectrum of mesons in the field theory, will be gapless and continuous. 

Now we turn to finite temperature, where we have two classes of embeddings. For 
Minkowski embeddings, the spectrum resembles the zero-temperature spectrum: it is gapped, 
discrete, and the masses are proportional to m/y/X [115]. Again, we can understand this 
just from the boundary conditions on fluctuations at the D7-brane's endpoint. 

The story for black hole embeddings, however, resembles the story from zero temperature 
and zero mass. Again, picture a wave traveling down the D7-brane. Such a wave can be 
absorbed, at least partially, into the black hole. In other words, we must impose a boundary 
condition that any such wave is purely outgoing (from the bulk, into the black hole). The 
spectrum will thus be not only gapless and continuous, but also contain eigenfrequencies 
with nonzero imaginary part [115]. 

Translating to the field theory language: at high enough temperature, the flavor fields 
undergo a first-order transition, and in the high-temperature phase the mesons acquire some 
width, which in fact increases as the temperature rises, suggesting that the mesons "melt" 
into the plasma [115]. (A broadening of quasi-particle peaks is also visible in spectral 
functions [116].) The first-order thermal phase transition thus leaves its most dramatic 
imprint on the meson spectrum. To emphasize this, starting now we will refer to the 
transition as a "meson melting" transition [115]. 

In the following sections, we will introduce additional thermodynamic variables, namely 
a density and electric or magnetic fields. We will characterize phases of the field theory by 
their meson spectrum. If the phase is described holographically by a Minkowski embedding, 
we will refer to the phase as a phase with stable mesons. If the phase is described by a 
black hole embedding, we will refer to IS clS cl phase with melted mesons. 



70 



4.2 Finite Density 

4-2.1 Only Black Hole Embeddings are Allowed 

As mentioned in section 3.2, the Nf flavors of hypermultiplet possess a global U{1) symmetry 
that we may identify as baryon number, and which we will denote as U{1)b- We will thus 
have a conserved current, J^, for baryon number. Again, let ip represent the Dirac fermion 
in the hypermultiplet and 4>a the two complex scalars. J'* then has the schematic form 

2 

J'' = iV^V + i'f'aD^'f'- - {D>'4>a)^<t>a) + h.C. (4.2.1) 

a=l 

with D'^ the gauge-covariant derivative. We want to study the field theory with a finite 
baryon number density. In other words, we want to study the theory in a state with (J*) ^ 0. 
Here we arrive at another question of the "universal qualitative" variety: what does a finite 
baryon number density look like holographically? 

To answer this question, we need only recall one of the basic facts about the AdS/CFT 
correspondence from section 2.4: a global symmetry in the field theory will be dual to a 
gauge invariance of the gravity theory. In this case, the global symmetry is the U{\)b 
associated with the flavor fields, which must be dual to a U{1) gauge invariance associated 
with the D7-branes, which must be the D7-brane worldvolume gauge field. More explicitly, 
the field theory operator J'* will be dual to the D7-brane's gauge field A^, hence J* will be 
dual to At, so if we want a state with nonzero (J*), then we need a nontrivial At. 

At can in principle depend on all of the D7-brane worldvolume coordinates. What 
ansatz for At should we write? As we did for the worldvolume scalar, we make symmetry 
arguments. We of course want to preserve the rotation and time translation invariance of 
the field theory, so At should not depend on the field theory directions. For simplicity, we 
do not want At to depend on the directions. The remaining option is for At to depend 
on the AdS radial direction. We thus want At{r) (here we are using the radial coordinate 
of eq. 3.3.5). We will always work in a gauge where the component of A^ in the radial 
direction vanishes: A^. = 0. We then have an electric field, F^-t = drAt, on the D7-brane, 
pointing in the radial direction. 

Suppose we are working at finite temperature, so that we have both Minkowski and 
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black hole embeddings. We can see that a problem will arise if we "turn on" At{r) for a 
Minkowski embedding. We have electric field lines pointing in the radial direction. What 
happens at the D7-brane's endpoint? The field lines cannot simply terminate: they must 
end on something. To turn the question around, we must ask what source is producing 
these field lines. 

The most natural source is a finite density of strings^. Recall again that a single string 
endpoint holographically represents a quark in the field theory. We want to study a finite 
density of quarks (and squarks), which means that in the D7-brane picture we want to study 
a finite density of strings, each with one endpoint on the D7-brane and the other endpoint 
piercing the horizon. The string endpoints on the D7-brane of course act as a source for 
the gauge field on the D7-brane, producing At{r). 

A straightforward analysis reveals, however, that the force such a density of strings 
exerts on the D7-branc is greater than the tension of the D7-branc [35]. We thus expect 
the strings to pull the D7-brane into the horizon, so that the electric field lines then end 
safely on the horizon. Actually showing that this occurs is difficult, however. Doing so 
requires computing a time-dependent solution of the D7-brane's evolution. What can be 
shown explicitly is that, with nonzero At{r), black hole embeddings are sufficient to cover 
the entire range of masses, from m = to m = oo [35]. Our intuitive, physical picture thus 
seems to be consistent, so we conclude that when we introduce At{r), we need only include 
black hole embeddings in our analysis. 

An interesting question, though, is what large-mass solutions look like. A large-mass 
Minkowski embedding is a D7-brane that ends near the boundary. When we turn on At(r), 
however, we are claiming that the D7-brane must somehow reach all the way to the horizon. 
How can this be? An analysis of the D7-brane action reveals that in such a limit, the D7- 
brane develops a "spike," in which the almost collapses to zero volume somewhere near 
the boundary, but then remains finite- volume all the way to the horizon [35]. In fact, the 
action for the spike, as derived from the D7-brane action, is identical in form to the action 



^Another possible source is a baryon, which in AdS/CFT appears as a D5-brane wrapping the factor 
of AdSr, X . Solutions for the D5/D7 bound state are not yet known [110], however, so we will not 
consider this possibility. 
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of a bundle of strings (that is, a number (J*) of strings) [35]. 

Crucially, black hole embeddings force upon us a particular boundary condition for At{r). 
At the horizon, the Euclidean time circle collapses to zero size. For ^4^ to be well-defined 
as a one-form there^. At must vanish: At{rh) = 0. Notice that this boundary condition, 
forced upon us by the geometry of the space, restricts gauge transformations, as only the 
subset of gauge transformations that preserve At{rh) = are allowed. 

4-2.2 D7-brane Action and Gauge Field Solution 

We are now ready to write the D7-brane action, including the worldvolume electric field 
Fj-^t = drQAt{rQ) = A'i.{re), for which we will solve. We will use the bulk metric 

ds'^ = - h{re) dt^ + rl dx^ + dO"^ + sin^ 9 ds% + cos^ 6 ds% (4.2.2) 

where /^(re) = Tq 1 and is the black hole horizon, related to the temperature by 

Th = ttT. Notice also that we have chosen Lorentzian signature. 

We will describe the D7-brane embedding by 6' (re). Introducing the worldvolume gauge 
field, the D7-brane action becomes 

S = -M j dre 4 cos^ 9 ^1 + /^(re) ^'^ - (27ra')2 (4.2.3) 

The action only depends on the derivative of At, so we have a constant of motion, which 
for the moment we will call D, 

D = Nrl cos' 9 ('^^^')'^t (4.2.4) 

+ h{re) ^'2 _ (27ra')2 A'^ 

We can now solve for ^^(re), 

D 



A'tire) 



1 + h{re) 9'^ 



(4.2.5) 



Inserting this solution back into the action, we find 
S = -M^ J dre 4 cos^ 6 



1 + hire) 9'^ 



\^'rlcos^9 + j^ 



(4.2.6) 



^The precise statement is: the event horizon of the geometry is a KiUing horizon which contains a 
bifurcation surface where the Kilhng vector dt vanishes. For the gauge field to be regular as a one-form, 
we must have At{rh) — [35]. 
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We may obtain the equation of motion for ^(rg) in two ways. We may either derive it from 
eq. (4.2.3) and then plug in the gauge field eq. (4.2.5), or we may eliminate the gauge 
field at the level of the action via a Legendre transform and then derive 0(r6)'s equation of 
motion. The Legendre-transformed action, which we denote S, is 

S = S - j dreFr.t ^ = " / '^^e + Kr^) cos^ e + j^^- (4.2.7) 

We would like to translate into field theory quantities. Recalling that z = l/re, we 
know from the analysis in section 3.4 that ^(re) will have an asymptotic expansion with 
the leading, non-normalizable term going as l/rg and the subleading, normalizable term 
going as l/rf. The coefficient of the leading term will be, up to a factor of (27ra'), the 
hypermultiplet mass ?Ti, and the coefficient of the sub-leading term will determine (C^^). 
In particular, these statements do not change in the AdS-Schwarzschild geometry, or in the 
presence of At{rQ), as can be confirmed explicitly from ^(r6)'s equation of motion. 

We also know that we should interpret the on-shell action, eq. (4.2.6), up to a Wick 
rotation to Euclidean signature and a factor of temperature, as the Gibbs free energy (grand 
canonical ensemble) and the Legendre transform, eq. (4.2.7) as the Helmholtz free energy 
(canonical ensemble) [103,104]. 

How do we extract (J*) and the chemical potential from the supergravity description? 
At{rQ) is dual to J*. The coefficient of the leading, non-normalizable term in At{rQ)'s 
asymptotic expansion should give the source for J*, which is just the baryon number chem- 
ical potential, fi. The coefficient of the sub-leading, normalizable term should determine 
the expectation value (J*). Let us see how this works. Prom eq. (4.2.5) we find Aj(r6)'s 
asymptotic expansion, 

A'tire) = ^ \ + o(\) (4.2.8) 

Integrating in re, we find 

1 D 1 „ / 1 



The leading behavior of At{rQ) is a constant, which we identify as the chemical potential. 
Notice that, normally, such a constant would not be physical: we could make a gauge 
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transformation to remove it. In our case, however, the boundary condition At{rh) = 
forbids us from making such gauge transformations, so the constant is physical. In fact, we 
can write /x as the integral of a gauge-invariant quantity, the field strength: 

= At(oo) = / dreA'^ire) = / dreFtr, (4.2.10) 

Jrh J oo 

where the second equality relies on the fact At{rh) = 0. We can thus interpret /i as the work 
done to push a string endpoint against the radial electric field, Ftr^, from the boundary (the 
UV) to the horizon (the IR). This is a nice holographic version of a field theory statement: 
the chemical potential is the energy needed to add a charge to the system. 

The coefficient of the sub-leading term should determine the expectation value of the 
dual operator, (J*). The quick way to determine the exact relation is to invoke holographic 
renormalization. (A longer, but more physical, argument appears in ref. [35].) We introduce 
a regulator by cutting off rg at some large value A. We can then write the expectation value 
for any current component J'*, 

{J^) = lim ^ — ^ (4.2.11) 

where 7 is the determinant of the induced metric on the re = A slice. Notice that \/— 7 ~ A^. 
For J*, we want the variation with respect to At{rQ), with all other fields held fixed. Defining 
the Lagrangian via S = J dr^L, we have 

= / drG—-—drMt = D dredr,5At = D{SAt{A)-5At{rh)), (4.2.12) 

Enforcing 6At{rh) = we find = D and hence (J*) = D. 

4-2.3 Finite Density Phase Diagram 

Using the above holographic setup, or more precisely the Euclidean-signature version of 
the above setup, much of the phase diagram of A/" = 4 SYM coupled to massive M = 2 
hypermultiplets at finite baryon density has been mapped out, in both the canonical and 
grand canonical ensembles [35,57,58,64]. In this section we will just present a qualitative 
sketch of the phase diagram, in its current form. 
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{J^)/m^ vs. T/m 



ji/m vs. T/m 




(a.) 



(6.) 



Figure 4.2: Phase diagrams in (a.) the canonical ensemble and (b.) the grand canonical 
ensemble, as presently understood. 



We begin in the canonical ensemble, in which the density is fixed, i.e. is our thermo- 
dynamic control parameter, and hence an input, and the chemical potential is an output. 
By continuity we expect the transition observed at zero density (discussed in section 4.1.1) 
to persist to finite density, at least for small density. In other words, in our phase dia- 
gram, which is now two-dimensional, with axes (J*)/m^ and T/m, we expect a line of phase 
transitions to emerge from the transition point on the {J*)/m^ = axis. 

The analysis of ref. [35] revealed that a line of transitions does emerge from the ( J*)/m^ = 
axis, but eventually ends in a critical point. This is depicted schematically in figure 4.2 
(a.). In the figure, the vertical axis is (J*)/m'^ and the horizontal axis is T/m. The solid 
line represent the line of first-order phase transitions that emerge from the transition point 
on the horizontal axis. The black dot represents the critical point. 

Immediately to the left of the line of transitions is a region of thermodynamic instability, 
discovered in ref. [35] and indicated in the plot with dashed lines (distinct from the dotted 
line, which we will discuss shortly). In this region, the condition > 0, for stability of 
the system against fluctuations in the density, is violated. We must then ask what the true 
ground state is in this region. The grand canonical ensemble will offer some insight into 
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this question. 

An important lesson emerges from the hne of first-order transitions at finite density. 
Recall that at zero density, the first-order transition appeared in the bulk as a topology 
change: the D7-brane transition was from a Minkowski to a black hole embedding. Here, 
however, as we only have black hole embeddings, the transition is between a black hole em- 
bedding and another black hole embedding, so the topology does not change. Nevertheless, 
the transition is discontinuous. 

What about the grand canonical ensemble, where /x is the control parameter and (J*) 
is an output of the dynamics? Here we learn an important lesson. To describe a finite 
density in the field theory, we must use black hole embeddings of the D7-brane. The field 
theory includes states, however, with finite jj. but zero density. In a theory with a mass 
gap, chemical potentials of magnitude less than the mass gap do not lead to any finite 
density: if the lightest charge carriers have mass m, then a chemical potential fj, < m 
cannot produce a nonzero density. In the canonical ensemble they never appear, since zero 
density corresponds to fi = m. In the grand canonical ensemble we should of course be able 
to dial fj, all the way down to zero, but the physics for /x < m is still expected to be trivial 
in this range, with (J*) = and the free energy being independent of fi. 

What D7-brane embeddings describe states with ^ > but (J*) = 0? The answer is: 
Minkowski embeddings [57,58,64]. How can this be? The D7-brane action only depends on 
dertivatives of the gauge potential, in this case, derivatives of ^((re). A valid solution of 
AtireYs equation of motion, in the absence of sources, is a constant, At{rQ) = /x. Our nice 
physical intuition of as the work needed to bring a string cndpoint from the boundary to 
the horizon fails here, but ^ is nevertheless the leading asymptotic value of At{rQ), and hence 
the chemical potential. Most importantly, D = for such a solution. We can thus have a 
Minkowski embedding with At{rQ) = fj,, describing a state in the field theory with finite /x 
but zero density. Notice also that the value of the on-shell action will be the same as the 
At{re) = solution, indicating that in the field theory the free energy is independent of /x, 
as expected. In the grand canonical ensemble, we must include such ^((re) = /x Minkowski 
embeddings if we wish to include all physically allowed values of /x. 

Including such embeddings, the phase diagram in the grand canonical ensemble was 
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computed in refs. [57,58,64], and is depicted schematically in figure 4.2 (b.). The vertical axis 
is /x/m and the horizontal axis is T/m. The solid line indicates a line of first-order transitions 
emerging from the known transition on the horizontal axis. The solid line intersects the 
vertical axis precisely at /u/m = 1. In supergravity language, the transition across the 
line is from Minkowski embeddings with constant ^((re), below the line, to black hole 
embeddings with nontrivial ^^(re), above the line. In field theory language, the transition 
is from a state with fi > but (J*) = 0, with a meson spectrum that is gapped and 
discrete, to a state with /i > and (J*) > 0, with a meson spectrum that is gapless and 
continuous. Notice, then, that the transition involves charge condensation, as (J*) jumps 
from zero to nonzero values across the transition line. Furthermore, the discontinuities 
in thermodynamic quantities across the phase transition line actually shrink as the system 
approaches the vertical axis [57] , raising the possibility that the line of first-order transitions 
may end in a critical point when T = and ji = m. In the next section, we will show that 
this is indeed the case. 

The two ensembles should be equivalent, that is, should describe the same physics. To 
explain what this means in practical terms, consider the canonical ensemble. Here we choose 
values of T and (J*), and the Helmholtz free energy then gives us the value of jj, = iJ,{T, {J*}). 
If, in the grand canonical ensemble, we choose the same T and the same value of /i, the 
Gibbs free energy will return the value of (J*) that we started with. We can thus think 
of the two ensembles geometrically, as two equivalent parameterizations of a surface, the 
collection of points (T, /x(T, ( J*)), ( J*)) in the canonical ensemble or (T, /x, ( J*)(r, /x)) in 
the grand canonical ensemble. This surface is the manifold of equilibrium states in the 
parameter space. We should ask whether the phase diagrams in the two ensembles are 
consistent with one another. 

If we map the line of first-order transitions in the grand canonical ensemble to the phase 
diagram in the canonical ensemble, we find the dotted line in figure 4.2 (a.) [57]. Below the 
line, Minkowski embeddings with constant At{rG) are thermodynamically preferred. These 
have (J*) = 0, so the grand canonical ensemble appears to have a gap. Notice also that the 
entire unstable region found in the canonical ensemble is within the region below the line. 

The question is thus how the system accesses the states below the dotted line in figure 
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4.2 (a.). The suggestion of ref. [57] was that the ground state in this region may actuaUy be 
an inhomogeneous mixture of the two phases, in analogy with the hquid-to-gas transition 
of ordinary water. Such a thing is difficult to describe in supergravity, however. It would 
be something like an inhomogeneous mixture of Minkowski and black hole embeddings. 
The important lesson is that the simple homogenous states that are easy to explore using 
supergravity are almost certainly insufficient to describe the entire phase diagram. 

We emphasize that most of the work that we have summarized was numerical: the phase 
diagrams were derived from numerical solutions for D7-brane embeddings. We will now 
present analytic solutions in two limits: zero mass and zero temperature. These solutions 
will nicely confirm the general picture outlined above. In particular, the zero-temperature 
analytic solution will reveal that the transition at /x = m is second order. 

4-2.4 Exact Solution at Finite Density: T > 0, m = 

Reviewing the phase diagrams in figures 4.2 (a.) and (b.), we can see that the m = 
limit actually looks pretty boring: it is the extreme upper right-hand corner of each plot, 
where wc expect no phase transitions. We will compute an exact solution for the D7-brane 
embedding scalar and gauge field in this limit, and explicitly confirm the absence of any 
interesting phase structure. 

We hasten to add, however, that the current phase diagram is probably incomplete. 
The key observation is that a chemical potential acts as a negative mass-squared for scalars. 
Suppose for the moment that we studied our theory at zero temperature, zero mass, and 
zero 't Hooft coupling, but with finite chemical potential. The theory will have no ground 
state: the potential for the scalars in the hypermultiplct would be unbounded from below. 
Introducing a finite mass or temperature can stabilize the potential, as can the quartic 
self-interactions and Yukawa couplings. In such cases we generically expect Bose-Einstein 
condensation: the hypermultiplet scalars will develop a nonzero expectation value, break 
the U{1)b symmetry, etc. How to see this in the supergravity description remains an open, 
and very important, question, one that we will not address here. We will work only with 
the D7-brane configurations described above. 
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Our goal is thus to integrate the solution for A[(rQ) in eq. (4.2.5) and the on-shcU action, 
eq. (4.2.6). To do so in general requires an explicit solution for 6 (re). We know the solution 
in the massless case: 9{re) = 0, which also describes (Cm) = in the field theory. 

With 9{rQ) = 0, the integral for At{rQ) becomes 



L»2 



,2 \ -1/2 



Mre) = I ( -6 + ) ■ (4.2.13) 

This integral can be performed analytically using incomplete Beta functions. We have 
collected the details in an Appendix to this chapter. The result is 

, , , 11/ D ( ( 4 I l\ ( rl 1 l\\ 



(4.2.14) 

Notice that At{rh) = 0, as it should be. The value of /i is given by taking re — >■ oo, for 
which the first term in parentheses goes to B « 8.41. 

The integral for the on-shell action is: 

^ = -^" j ^6 4 + . (4.2.15) 

The re integral is of course divergent. We regulate the divergence by cutting off the re inte- 
gration at some large value A. Rather than use counterterms, we will subtract a background 
solution. The simplest solution is a D7-brane with ^(re) = At{re) = 0, for which we denote 
the regulated on-shell action as Sq: 

So = -^f r dr^rl = - + \^frl (4.2.16) 

We then define the renormalized action Sren as 



lim {S-So), (4.2.17) 

A— »oo 



with the result 




Sren = -\Mrt - '-Af-^^ b{ ^ ^ | . (4.2.18) 

The result has no interesting features, i.e. is completely monotonic, and its derivatives 
display no discontinuities, which is consistent with the numerical results computed in refs. 

[35,57,58,64] and displayed in figure 4.2. 
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4-2.5 Exact Solution at Finite Density: T = 0, m > 

The main feature of the phase diagram of figure 4.2 (b.) is the hne of first-order transitions 
between a phase with a gapped, discrete spectrum of stable mesons and zero baryon density 
and a pliasc with a gapless, continuous spectrum of unstable mesons and finite baryon 
density. As the line of transitions approaches the T/m = (vertical) axis, the first order 
transition grows weaker: the discontinuities in the derivatives of thermodyanmic potentials 
shrink. A natural question is whether the transition eventually becomes continuous. In 
other words, does the line of first-order transitions end in a critical point when T = and 
H = ml In this section we will compute the thermodynamics precisely along the T/m = 
axis, and discover a second-order transition when fx = m. 

In supcrgravity language, our goal is to solve for the D7-branc fields, and use them to 
compute the D7-brane's on-shcU action. Even when T = 0, finding a solution for ^(rg) and 
then performing the integrals to find At{rQ) and the on-shell action is very difficult. 

The key to finding exact solutions is to change coordinates. We switch to the AdS^ x 
metric 



ds^ = Z{re) {-dt^ + dx^) + ^(re)"^ (dr^ + ds^s + dy^ + ds%) , (4.2.19) 



with Z{rQ) = Vq = r'^ + y^. The D7-brane scalar is now y{r) and the gauge field will be 
At{r). With this ansatz for the fields, the D7-brane action becomes 



Remarkably, the action depends only on derivtiaves of both y{r) and At{r). We thus have two 




(4.2.20) 



constants of motion. Defining the Lagrangian by 5 = J drL, we denote the two constants 



of motion as c and d: 





(2W)24 



d. 



Vl + y'2 _ (27ra')2Af 



Vl + y'2 _ (27ra')24^ 



(4.2.21) 



The ratio implies 




(27ra')^ c- 




(4.2.22) 
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so that if we can perform the integration for one field, wc automatically find the answer 
for the other. We can solve algebraically for y'{r) and Af{r) in terms of the integration 
constants c and d, 

y = / n =' ^* = ; ' =, (4-2.23) 

which can be integrated using incomplete Beta functions. The result depends on the sign 
of ^2^1^')'^ ~ When c = d = 0, we obtain the solution with y'{r) = and A[{r) = 0, so 
y{r) and At{r) are constants. If At{r) = 0, we recover the embeddings of section 3.3.2, with 
y{r) constant. We will investigate the other possibilities in what follows. 
The D7-brane action evaluated on these solutions is 

fA 



/ 



■M / drr 



3 



(4.2.24) 



The lower endpoint of integration will depend on the sign of (^2iTa')'^ ~ '^^^ integral 
diverges if we integrate to r = oo, of course, so we have regulated the integral with a cutoff 

at r = A. The divergence is already present for c = d = 0: 

/•A 1 

So = -Af drr^ = —M k"^ (4.2.25) 
JO 4 

so we will s-gciiii d.6fiiie th.6 r6iiorinctlized. oii-sh.6ll &.ctioii, Sf^n^ cis 

Sren = Hm {S - Sq) (4.2.26) 

A— >oo 

In the field theory, the thermodynamic potential Q of the grand canonical ensemble is 
given by^ O = —Sren- We can thus determine whether solutions with nonzero d and c are 
thermodynamically favored relative to the c = d = case as follows: Sren > means the 
configuration is favored and Sren < means the configuration is disfavored. 

We Legendre transform to obtain the free energy density in the canonical ensemble, 
F = O, + ii{J^). Recall that at zero temperature, free energy and energy are the same, so 
at zero temperature F is also the energy density. 



*In section 2.5, when T > and = 0, we made the identification F = TSsugra, and mentioned that 
if /X > we should identify Q = TSsugra- Here we have T = and /x > 0, however, for which the 
identification above is the correct one. 
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Figure 4.3: Cartoons of D7-brane embeddings, for various values of (2^')^ " ^ 



The conserved charges c and d determine (Om) and ( J ) as foUows, 

where in each case when we vary one field we hold the other fixed. We can write 

6S = Jdr (^.^^^drSAtir) + JL_dr6y{r)j = dMi(oo) - c6y{oo) (4.2.28) 

where we demand that 6At(r) and 6y{r) are always zero at the lower endpoint of the r 
integration. If we vary At{r) while holding y{r) fixed {Sy{r) = 0), we find (J*) = d, and 
similarly we find (Cm) = (27ra')c. Starting now, we will stick to the notation c and d and 
refer to these as the condensate and density. 

We now turn to the various possibilities for the values of c and d, and in particular 
the sign of ^^2^1')'^ ~ '^^^ possibilities are depicted, in a very schematic way, in figure 
4.3. In the figure, the vertical axis is one of the directions transverse to the D3-brancs (in 
Maldacena's construction) but parallel to the D7-branes, while the horizontal axis is an 
overall transverse direction. Figure 4.3 (a.) depicts the c = d = case, in which At{r) = 
and y(r) is a constant. Figure 4.3 (b.) depicts the (^2na')'^ — < case, and figure 4.3 (c.) 
depicts the (2na')'^ — c > case, which we will now describe in turn. 

We begin with f^2na'y^ — < 0, and furthermore, we restrict to d = 0. Such solutions 
were first studied in ref. [117]. These are, in fact, precisely the solutions we mentioned in 
section 3.3.2, below eq. (3.3.9), when we first discussed D7-brane embeddings. 
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For such d = embeddings, y'{r) diverges as ^^^^^ at a critieal radius ro = AA~^/^c^/^. 
This behavior indicates that the D7-brane has a turn-around point where the D7-brane 
smoothly matches onto a second branch. This solution thus describes a D7/anti-D7-brane 
pair connected by a smooth neck, as depicted schematically in figure 4.3 (b.). 

In the asymptotic region, we may specify the distances of the D7-brane and the anti- 
D7-brane to the D3-branes separately. In other words, we may associate two masses to 
such a configuration. Wc denote these as m and m. A D7-branc solution exists for any 
m and {Om)- The mass (or asymptotic separation) of the anti-D7-brane is then fixed in 
terms of these two input parameters. Prom the field theory point of view, we would like 
to interpret the two masses m and fn as input parameters, so that the above construction 
dynamically determines the value of the condensate (Om)- Prom the supergravity picture, 
we can see that the condensate only depends on the difference of the masses, m — fn. As 
we will show (and as first noted in ref. [117]), all such solutions have Sren < 0, and hence 
will be meta-stable, at best. The true ground state of the system has vanishing (Om)- 

Now we re-introduce nonzero d. To connect the two halves of the D7-brane smoothly we 
need the gauge field and its derivative to be continuous at the neck. In particular, the field 
strength on the two branches should be equal in magnitude but opposite in direction. The 
baryon number density is then equal and opposite at the two ends, so this configuration 
does not describe any net baryon number! The Lagrangian of the field theory has a global 
U{1) X U{1) symmetry, and, at least in the grand canonical ensemble, we can imagine 
turning on different chemical potentials for the two i7(l)'s at the level of the Lagrangian, 
even though the symmetry is spontaneously broken to the diagonal U{1)b in the state 
represented by the connected D7-brane configuration. If we wanted to have a branc/anti- 
brane system with finite baryon density we would have to include an explicit source at the 
neck. Such additional sources will back-react and alter the configuration, so we will not 
discuss this option further. 

In any case, these brane/anti-brane systems correspond to a different field theory La- 
grangian than that for systems with a single D7-brane reaching the asymptotic region. The 
dual field theory for the brane/ anti-brane is A/" = 4 SYM coupled to two hypermultiplets that 
preserve opposite N = 2 supersymmetries, so that the full theory is non-supersymmetric. 
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ro=^f-'/'{^'-77^2] (4-2.29) 



When comparing all possible bulk configurations that correspond to hypermultiplets pre- 
serving M = 2 supersymmetry, the brane/anti-brane configurations never contribute. 

For nonzero d, keeping in mind that for the brane/anti-brane system this does not 
correspond to a net baryon density, the D7-brane reaches the turn-around point 

{2na' 

The integrals for y{r), Mir) and the regulated on-shell action are performed in the Ap- 
pendix. We find 

and At{r) is simply J2tu7y\ times y(r), as shown in eq. (4.2.22). Notice y(ro) = At{rQ) = 0. 
The mass and chemical potential are given by the asymptotic values at r — >■ oo, with 
limr_*oo-B |, =0. The renormalized on-shell action for the D7-brane is 

S„„ = -i^-V3(^-^)^'^(-|i) (4.2.31, 

The on-shell action for the anti-D7-brane is identical, so the total action is twice Sren- 
Notice that Sren < 0, so these connected brane/anti-brane configurations always have a 
higher free energy than the c = d = case, and hence are thermodynamically disfavored. 

At this point we should also consider c and d nonzero, but such that (^2na'y^ — = 0. 
This configuration shares a similar feature with the brane/anti-brane configurations: it does 
not satisfy the required UV boundary conditions for a theory with a single D7-brane added. 
When ^2mx')'^ ~ ^ =0) wc have y'(r) ~ ^j, so the D7-brane misses the D3-branes at the 
origin and returns to the asymptotic region at r ^ oo. We will therefore not study the 
c = case further. 

Now we come to the interesting case: {^f^i^ji — > 0. We will show that these are the 
zero-temperature limit of black hole embeddings. Notice first that these solutions have no 
turn-around point, and describe D7-branes alone. 

For a ^2^0.')^ — > solution, suppose we take c = 0. The solution with constant y(r) 
is still allowed. Geometrically, this is the zero-temperature Minkowski embedding, but now 



85 



with nonzero A[(r). As mentioned before, the constant charge density along the D7-brane 
requires a source at r = 0, so this embedding is not physical unless we add extra strings 
connecting the D3-branes and D7-branes, or some other source. 

We can avoid this issue, however, by demanding that the D7-brane touch the horizon at 
r = 0. Since the horizon is located at = r| = + this immediately implies a boundary 
condition on y, namely y{r = 0) = 0. The embedding satisfying this boundary condition is 
the zero temperature analogue of black hole embeddings. We will stick to this name, even 
though at zero temperature "horizon crossing" would be more appropriate. This black hole 
embedding is displayed schematically in figure 4.3 (c). 

The integral for y'(r) is done in the Appendix, with the result 

- 4^-"' ( - - i 

and again, At{r) is (^S^ryi^ times y{r). Notice that y(p) = At{0) = 0. Identifying m and /x 
from the asymptotic values y{oo) and At{oo), where as r — oo the incomplete Beta function 
becomes S (g, we find 

c = jX{27ra'f (//^ - m^) m (4.2.33a) 
= -fM(2Tra'f ii? - rr?) (4.2.33b) 

where we have defined the constant 



7=(isQ,i)) « 0.363. (4.2.34) 



Notice that since these black hole embeddings were valid only when (^2^^')'^ — c > 0, we can 
immediately see that they only exist for > m, as expected for a theory with a mass gap, 
in a state with nonzero density. 

The renormalized on-shell action in this case is 

2/3 



V-'"^-'" - ^^)''' = (m^ - rr?? (4.2.35) 



Notice that Sren > 0; so these embeddings are thermodynamically favored relative to the 
c = d = case. Notice that we observe no instability in O. Stability requires |^ > 0, which 
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is clearly satisfied in eq. (4.2.33). This is consistent with the numerical results of ref. [57] 
in the zero-temperature limit. 

The thermodynamic potential in the canonical ensemble is 

F^n^ ,iJ') ^ 1^- V3,- V3 _ (3^ + c=) (4.2.30, 

so that F > 0. The black hole embedding still turns out to be the configuration with lowest 

F, however. 

To see the phase transition between black hole and Minkowski embeddings, let us return 
to c = and constant y{r). These are solutions everywhere away from r = 0. These 
D7-branes do not obey the boundary condition y(0) = 0. In other words, when r = the 
D7-branes are still a distance y away from the D3-branes. 

The integral for At{r) will be unchanged, however. In particular, we still integrate from 
r = to r = 00, and At(0) = 0. Prom eq. (4.2.33) we read off the relation between ^ and d, 

d = 7A/'(27^a')V^ (4.2.37) 

The integral for the on-shell action is also unchanged, so we simply set c = in the first 
equality of eq. (4.2.35), 

Sren = ^7-'/'AA-'/'(2W)-^/3^^/3 ^ 1^^(2W)V (4-2.38) 

and hence J7 = —Sren ~ Naively this J7, for any nonzero d (or /j, > m), is smaller 

(more negative) than the one for the black hole embedding, eq. (4.2.35). Notice, however, 
that unlike the nice interpretation we found for the relations among c, d, n and m for the 
black hole embeddings, this result for fi, a pure quartic in ^u, is completely independent 
of the mass and therefore appears to be unphysical. For example, in the decoupling limit 
m — > DO, should go to zero, while this result clearly does not. 

The problem of course is that the Minkowski embedding, as it stands, is not consistent: 
we completely ignored the contribution to J7 due to whatever object sources the D7-brane 
gauge field. The simplest source is a finite density d of fundamental strings stretching from 
the D3-branes to the D7-brane. These D7-branes with strings attached are not solutions 
to the full equations of motion following from the combined DBI and Nambu-Goto actions. 
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however. Naively, we might think that the string back-reaction on the D7-brane can be 
neglected in the large-A^c and large-A limit. For a single string this is certainly true. The 
prefactor of the Nambu-Goto action scales as \/A, while the prefactor J\f of the D7-brane 
action scales as NfNf.X. We want a finite density d of strings, however. If we keep all 
geometric distances of order one in units of the AdS curvature radius, the mass will scale as 
\/A, so we want the chemical potential also to scale as VX- Prom eq. (4.2.33), we can see 
that the density d will then scale as NfNc\^. With this, d times the Nambu-Goto action 
has the same NfNcX scaling as the D7-brane action, and so the back-reaction is order one. 

To include the back-reaction we would have to re-solve the equations of motion including 
the extra source term. We already found the most general solution to the equations of motion 
with sources localized at r = 0, so the back-reacted solution including any such sources 
must be within this class. The only well-behaved solution in this class which asymptotically 
becomes a single brane is the black hole embedding, so it must be the back-reacted solution. 

To show that this is physically reasonable we can show that the brane-plus-strings con- 
figuration, which could serve as a consistent initial data for a full time-dependent physical 
solution, has higher energy than the black hole embedding. Wc add a term dm to the action 
representing the finite density d of strings with length y = (27ra')m. The full free energy 
including the D7-branes and strings in this case is 

which is positive at /x = m and is disfavored relative to the black hole embedding at all 
values of n > m. We can also allow some mixture, where the D7-brane embedding satisfies 
the boundary condition y{r = 0) = yo for some m > > 0, with a density d of strings 
extending from y = to y = yo- In figure 4.3, we imagine sliding the black hole embedding 
to the right a distance yo, such that the total asymptotic separation remains y, and then 
introducing the strings. The free energy is then 

and we can easily check that in the relevant regime /i > m > > the derivative of this 
expression with respect to yo is strictly positive. The D7-brane with strings attached can 
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continuously lower its energy until it turns into the black hole embedding with all strings 
dissolved in the D7-brane. 

While none of the finite-density Minkowski embeddings solve the equations of motion, 
the trivial Minkwoski embeddings with constant At{r) do. As only Frt enters the action 
and not At{r) itself, these embeddings have F = Q, = d = 0, and is a free parameter. For 
fi < m these are the only allowed configurations, so they dominate the ensemble. The fact 
that they are indistinguishable from the vacuum is completely natural from a field theory 
point of view, since the theory has the mass gap m. When /j. > m these solutions still exist, 
but in that regime the black hole embedding, with negative Q,, dominates the ensemble. 

At /X = m we hence have a second order phase transition between the trivial Minkowski 
embedding and the black hole embeddings. Both Q and its first derivative with respect to /j, 
are zero at that point, but the second derivative jumps from zero, for the trivial Minkowski 
embedding, to —2^J\f {27ra')'^m'^ , for the black hole embedding. The critical exponents are 
given by their mean-field values, that is, for /x = m -|- e we have ft ^ e^, c ^ e, d ^ e. 

We have thus shown that, indeed, the line of first-order transitions in the phase diagram 
of figure 4.2 ends in a critical point on the T = axis, explicitly confirming the numerical 
results of ref . [57] . 

4.3 Finite Electric and Magnetic Fields 

In the next chapter we will study the transport properties of fiavor fields using the D7-brane 
description. We will introduce external electric and magnetic fields and then compute the 
resulting charge currents, from which we can extract a conductivity tensor. Before studying 
transport, however, we must study thermodynamics: we must know what the equilibrium 
state is before we perturb around it. In this section we very briefiy review what is currently 
known about the thermodynamics of the fiavor fields in the presence of constant electric and 
magnetic fields [59-63]. We stress that all work to date regarding the thermodynamics in 
the presence of external fields has been done with zero baryon density: the phase diagram 
in the full parameter space (m, T, (J*), E, and B) remains largely incomplete. 
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4-3.1 Electric Field 

How can we introduce an electric field in the field theory? The flavor fields have the baryon 
number symmetry U{1)b- We can introduce an electric field that couples only to excitations 
that carry this U{1)b charge. Notice that we will not gauge the ^7(1)^: wc will not have 
a dynamical "photon." The electric field we introduce will be an external parameter, one 
that is only "felt" by fields carrying U{1)b charge. 

The primary motivation for introducing an electric field in this fashion is because the 
supergravity description is relatively easy. In the field theory, we want a nonzero field 
strength F*^ = E, describing a constant electric field E in the x direction. Under the right 
circumstances, we expect a current in the x direction, that is, an expectation value (J^). 
The former we expect to arise as the leading, non-normalizable mode of some D7-brane 
field, while the latter we expect to arise from a sub-leading, normalizable mode. The U{1)b 
in question is described by the U{1) gauge field on the D7-brane worldvolume. We can thus 
introduce E and ( J^) by assuming an ansatz for the D7-brane gauge field that includes [65] 

A^{t,re) = -Et + Ure) (4.3.41) 

with fxire) an rg-dependent function for which we must solve. The first term produces 
the F*^ in the gauge theory, while we should be able to extract (J^) as a coefficient of a 
sub-leading term in fxireYs asymptotic expansion. We may thus study the field theory at 
finite temperature, baryon density, and finite E by adding Ax{t,rG) to our above ansatz for 
the fields. Currently, however, the thermodynamics with a finite electric field has only been 
studied at zero density. 

Wc can actually anticipate what will occur with a finite electric field, using our field 
theory intuition [63]. Suppose we are at zero density, and with a finite temperature below 
the meson melting transition. The low-energy degrees of freedom are then stable mesons. 
We expect that an electric field will tend to pull a meson apart: if the meson is a quark 
and anti-quark, the electric field will pull the quark one way and the anti-quark the other. 
In other words, the binding energy of the quark is reduced in the presence of an electric 
field. We should then not need to heat the system up as much to see mesons melt. Our 
expectation is thus that the meson melting temperature should decrease. 
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Once the mesons have been pulled apart and/or melted, we also expect a current, (J^), 
since we will then have free charge carriers (quarks and squarks) in the presence of an 
electric field. The meson melting transition in this case may be called an "insulator-metal" 
transition, at least in the sense that charge carriers are liberated. 

Overall, the supergravity analysis is consistent with the above expectations, though 
with some minor subtleties. Prom the supergravity point of view, the electric field on the 
D7-brane worldvolume tends to "push the D7-brane into the horizon," facilitating the first- 
order transition from Minkowski to black hole embeddings. Some D7-brane embeddings were 
found (numerically) to be singular, however, meaning these embeddings, as they stand, are 
probably unphysical [62,63]. A small part of the phase diagram thus remains in question. 
The status of these embeddings is still unresolved. 

The result for the phase diagram, as derived from D7-brane solutions [62,63], is depicted 
schematically in figure 4.4 (a.). In the figure the vertical axis is E/T"^ while the horizontal 
axis is T/ra. The solid line represents the phase transition from Minkowski to black hole 
embeddings: Minkowski on the left, black hole on the right. In field theory language, the 
region to the left of the solid line is a phase with stable mesons and no current (J^), while 
to the right of the line is a phase with unstable mesons and a net current. The solid line 
emerges from the known transition point on the horizontal axis, where EjT"^ = 0. Clearly 
the transition temperature decreases as E/T'^ increases. The region marked with dashed 
lines corresponds to singular D7-brane embeddings. The ground state in this region is 
currently unknown. 

4.3.2 Magnetic Field 

We introduce a magnetic field in the field theory is essentially the same way as we introduced 
the electric field: as a non-dynamical external field that couples to anything carrying U(1)b 
charge. We then want a nonzero F^^ = B. If we also introduce E as above, we expect a 
current (J^) but also now, via the Hall effect, (J^). To our ansatz for D7-brane field we 
should thus add [66] 

Ay{x,rQ) = B X + fy{rQ) (4.3.42) 
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The leading term will produce F^^ while fy{r%) should produce (J^). In this section, 
however, we consider only = 0, so that the currents are zero. 

Two interesting things occur with nonzero B. The first occurs at zero temperature and 
zero hypermuhiplct mass: chiral symmetry is spontaneously broken [59,60]. We mean by 
this that, even with zero temperature and zero mass, nonzero B produces a nonzero (Cm)- 
As {Ojn) is charged under the U{1) that acts on the quarks as a chiral symmetry (dual to 
rotations in the X^-Xq plane), we call this chiral symmetry breaking. A Goldstone boson 
should then appear in the meson spectrum, which a supergravity analysis confirms [59,60]. 
The Goldstone boson obeys the Gell-Mann-Oakes-Renner relation, as it should [59,60]. 
Additionally, the spectrum of massive mesons exhibits Zeeman splitting [59,60]. 

The second interesting thing occurs at finite temperature: we find a critical value of 
the magnetic field above which mesons are always stable [61,62]. In heuristic terms, the 
magnetic field acts to "hold mesons together." Suppose we fix m and T such that the 
system is in a phase with stable mesons. As B increases, no transition will occur. If instead 
we choose m and T such that the system is in the phase with unstable mesons, then as B 
increases we will eventually see a transition to the phase with stable mesons. 
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Prom a supergravity point of view, the black hole of the AdS-Schwarzschild background 
pulls on the D7-brane, while the worldvolume magnetic field acts to pull the D7-brane away 
from the black hole, so that the temperature of the transition from Minkowski to black 
hole embeddings increases. For fixed m and T, Above a critical value of B only Minkowski 
embeddings exist. These embeddings can then be shown to describe the entire range of m, 
and indeed, for m = 0, the supergravity solution produces {Om) 0, hence chiral symmetry 
is spontaneously broken. 

The result for the phase diagram is depicted schematically in figure 4.4 (b.). In the 
figure the vertical axis is iJ/T^ and the horizontal axis is Tjm. We have included in the 
plot the value Tjm, = oo, indicated with the dashed line on the right. The solid line 
represents transitions from Minkowski embeddings, on the left, to black hole embeddings, 
on the right, or in field theory language, transitions from states with stable mesons to states 
with unstable mesons. The line emerges from the known transition point on the horizontal 
axis, where B /T^ = 0. 

Clearly, if we start on the B/T'^ = axis, to the left of the transition point, and then 
increase B/T"^, moving up in the figure, we encounter no transition. If we start to the 
right of the transition point, we can see that eventually we encounter a transition to the 
Minkowski/stable meson phase. Indeed, even when T/m = oo, which we could achieve with 
T finite and m = 0, if we increase past a critical value, indicated in the figure by the 

dotted line, a transition occurs to the Minkowski/stable meson phase. 

Appendix: Beta functions and Incomplete Beta functions 

The integrals we need for this chaper all take the form of Beta functions or incomplete Beta 
functions. Here we collect some of the important properties of these functions, and perform 
the integrals necessary for sections 4.2.4 and 4.2.5. We use the notation of section 4.2.5, 
and make note of which integrals are relevant for section 4.2.4. 
The Beta function is defined as 
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and the incomplete Beta function as 

B {x; a, b) = r dt{l - t)^-^"-^ = f'^^ du{l + (2) 
Jo Jo 

These satisfy the recursion relation 

B{x;a,b) = B{a,b) - B{1- x;b,a) (3) 

and are related to the hypergeometric function as 

B{x; a, b) = a~^x"-F{a, 1 - 6; a + 1; x). (4) 

The divergence of S in section 4.2.5 appears in the expansion, for a = —2/3 about 
x = 0, 

B(x;-|6)=-^x-2/3 + 0(a;V3) (5) 

where for S we will have x ~ (see for example eq. (8) below). 

The integrals we need can easily be brought into the form of eq. (2). For notational 
simplicity, we will replace (^2na'y^ with just (P. We begin with the (P — < case. The 
integral needed for y{r) and At(r) is 

7<(r) = fdr {Af^r'' - (c^ - d^))"'/' , (6) 

Jro 

with ro = J\f~^^^ (c^ — d'^Y^^. We change variables to t = r^/r^, 

I<(r) = ^N-^l\^-d^)-^/^ j\lt{l~t)-^'h-^l^ (7) 

The same change of variables works for the regulated on-shell action. 

Now wc turn to the d? — (? > case of section 4.2.5. The integrals here arc also relevant 
for section 4.2.4. The integral needed for y{r) and At{r) (and of the same form as the 
integral for At{rQ) in section 4.2.4) is 

(r) = f dr {J\fr^ + d^ - c^) . (9) 
Jo 



94 



We change variables to u = ^ri^, 

/^(r) = lAf-"='(<i=-<;^)"'" r<J«(l + «)-"\-=/« (10) 

= V./3(,2_,.)-i/3j( ATV 1 ly 

6 ^ ^ VA^^'^ + t^ -c^ 6 37 

The integral for the regulated on-shell action (also relevant for the on-shell action in section 
4.2.4) is, with A' = ^ 

o Jo 
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Chapter 5 

TRANSPORT OF FLAVOR FIELDS 

Having seen some of the equihbrium thermodynamics of flavor fields, we are now ready 
to study the transport of flavor flelds, that is, we are ready to study perturbations about 
equihbrium. Such perturbations are characterized by transport coefficients, such as the 
shear and bulk viscosity, diffusion constants, and conductivities. In this chapter we will use 
the supergravity description of flavor fields to compute their conductivity tensor. We first 
introduce electric and magnetic fields, and then use the D7-brane description to compute 
the resulting currents. From Ohm's law, we can then simply read ofi" the conductivity. As 
an added bonus, we can compute the drag force experienced by the flavor fields as they are 
pushed through the M = 4 SYM plasma, or more accurately the drag coefficient, in the 
limit of large hypermultiplet mass m, where the flavor fields behave as quasi-particles. 

Most of the material in this chapter comes from refs. [65,66]. 

5.1 Preliminaries 

To start, we recall some basic definitions. The conductivity tensor aij measures the response 
of a conducting medium to externally applied fields. It is defined by 

(Ji) = (^ijEj 

where Ej arc externally applied electric fields and {Ji) are the currents induced in the 
medium. An external magnetic field B produces off-diagonal elements in aij-. the induced 
current is perpendicular to both E and B. This is the Hall effect. For a rotationally- 
invariant system with E in the x direction and B perpendicular to the xy plane, axx = (^yy 
and Uxy = —CTyx- The component axx is called the Ohmic conductivity and cr^y the Hall 
conductivity. 

We will now review two results from classical electromagnetism that we will reproduce 
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from our supergravity calculation in appropriate limits. These will provide some nice checks 
of our supergravity answer for aij. 

Imagine filling the vacuum with a charge density (J*). In the lab frame we may introduce 
a magnetic field B. In a frame moving with velocity —v relative to the lab frame we will 
find a current J = ( J*)v and an electric field 

E = -vx B = --^Jx B. (5.1.1) 

If we take B = (0, 0, B) we find the conductivity 

(Txx = 0, a^y = {J')/B. (5.1.2) 

Notice that this argument does not require that the charge density be comprised of quasi- 
particle charge carriers. Indeed, this argument relies only on Lorentz invariance. 

Now imagine a density (J*) of massive quasi-particles propagating non-relativistically 
through an isotropic, homogeneous, neutral medium. In the rest frame of the medium we 
introduce an electric field E in the x direction in addition to the magnetic field. The force 
on a quasi-particle is then 

^ = E + vxB-iidP, (5.1.3) 

where our quasi-particle has charge +1 and is the drag coefficient (hence the subscript, 
to distinguish it from the chemical potential). We replace the momentum with the velocity 
using p = Mv for quasi-particle mass M. We then replace the velocity with the induced 
current using v = (J)/(J*). Imposing the steady-state condition ^ = and solving for {J) 
yields 

<Jo _ _ ao{B/fidM) ^^^^^ 



(B//i,M)2 + l' {B/^idMY + l 

where ctq = {,P)/iJ,dM is the conductivity when B = 0, and is also known as the Drude form 
of the conductivity, although Drude would have called (J*) the number density of electrons 
and M the electron mass. 

5.2 The Physics of Flavor Transport 

We want to study the transport properties of a number Nf <^ Nc of Af = 2 hypermultiplet 

fields propagating in an A/" = 4 SYM plasma, at large- iVg and large 't Hooft coupling A. We 
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will introduce a finite density (J*) of flavor fields, and electric and magnetic fields E and B 
that couple to U{1)b charge. We will thus think of the density (J*) as a number density of 
charge carriers. We will define the baryon charge of the flavor fields to be +1. 

The physical picture of what will happen in the presence of nonzero E is fairly intuitive. 
The A/" = 4 SYM fields will not "feel" the electric field since they do not carry U{1)b 
charge. The Af = 2 hypermultiplet fields, however, will be pushed by E. They will still be 
interacting with the = 4 SYM fields, of course, with coupling strength A. The = 4 
SYM plasma will thus provide a drag force, as the charge carriers must push through the 
plasma. We may hope that when the force due to the electric field balances the drag force 
we will find a steady-state current ( J^). With nonzero B, we also expect a current ( J^) from 
the Hall effect. Our goal is to compute these currents, using the supergravity description, 
and extract from them the components of the conductivity tensor. 

Ultimately, the simple physical picture described above will turn out to be correct, but 
the reasons why this is so are actually quite subtle. To explain why, let us consider for a 
moment the same system, but with Nf on the order of Nc- In this case, we do not expect 
a finite answer for aij. Why not? The two crucial features of the system are: 

1. ) The system has a finite charge density in a constant electric field, and 

2. ) The system is translationally invariant. 

The first point tells us that the electric field will be doing net work on the system. The 
system will thus be gaining energy at a constant rate as time passes. The second point 
tells us that momentum will be conserved in the system. As the electric field does work 
on the charge carriers, they will gain energy and momentum. If momentum is conserved, 
the charge carriers will transfer momentum to the = 4 SYM plasma. Put more simply, 
the charge carriers will begin to drag the plasma along with them. Our constant electric 
field will continue to do work on the system, dumping in more and more energy, so that 
over time the entire system will move faster and faster, without bound. The end result is 
definitely not a steady state with constant currents (J^) and (J^). Instead, we end up with 
a plasma moving infinitely fast! In a word, the problem is dissipation, or rather, the lack 
thereof: our system seems to have no mechanism by which the charge carriers can dissipate 
momentum, which is necessary to produce a finite conductivity. 
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To address this problem, we could use an electric field with harmonic time dependence, 
E{t) ~ e"^*, in which case the electric field will do no net work on the system. Indeed, 
other gauge-gravity systems have been studied with harmonic fields, giving rise to finite, 
frequency-dependent conductivities [118-121]. We will work with constant E, however. 

We could also try to break translation invar iance, as occurs in real materials. A crude 
picture of a metal, for example, is a gas of conduction electrons propagating through a 
lattice of ions, where the lattice clearly breaks translation invariance. The system then does 
not conserve momentum, and the conduction electrons can experience dissipation. Impurity 
potentials that break translation invariance have been introduced in gauge-gravity systems, 
giving rise to dissipation and hence finite conductivity, even in a constant external electric 
field [122]. We will not break translation invariance, however. 

How on earth can we find a finite conductivity, then? The answer comes from the probe 
limit, in which we can basically fake dissipation. In the probe limit, which in particular 
means Nf Nc, the stress-energy tensor of the system will cleanly separate into two terms, 
one of order N^, coming from the adjoint fields, and one of order NjNc, coming from the 
fiavor fields. We depict this schematically as 

T^, ~ 0{N^)f,, + 0{NfN^)^,. (5.2.5) 

The flavor fields may then transfer momentum to the A/" = 4 SYM fields at a constant rate 
for a time that is parametrically large in A^c, that is, a time on the order of Nc- At earlier 
times, the velocity of the A/" = 4 SYM plasma will be negligible. Only at very late times will 
the flavor flelds have transferred enough momentum to the plasma for the velocity of the 
plasma to be non-negligible. For any finite time, then, the A/" = 4 SYM plasma basically 
acts as an energy-momentum reservoir for the flavor fields. In simple terms, all we are 
saying is that because the system has so few flavor flelds, we can neglect their effect on the 
plasma, at least for a very long time. 

In the end, then, we will flnd a flnite conductivity, even with flnite charge density, 
constant electric field, and translation invariance, because we work in the probe limit. We 
should always bear in mind, however, that our result will not be valid for all time. 

If we take m ^ oo, making the flavor flelds very heavy, then we expect the charge carriers 
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to behave as quasi-particles: as m — >■ oo, the Compton wavelength of the flavor excitations 
will become arbitraily shorter than their mean free path, so that they will behave essentially 
as isolated particles. In other words, the charge carriers will be very heavy probes lumbering 
through the plasma, their dynamics essentially classical. 

In the m — DO limit, then, the quasi-particles will obey a classical force law, 

^ = -HdP + E + vxB (5.2.6) 

where fi^ is a drag coefficient and p is the momentum of the quasi-particles. Here we see 
explicitly the drag force due to the J\f = 4 SYM plasma and the Lorentz force due to the 
external fields. 

We must actually be more precise with what we mean by "m — ^ oo." To explain why, 
let us consider a single quark, rather than a density (J*) of quarks (and squarks). A single 
quark is represented in the supergravity description as a single string stretched from the 
endpoint of a Minkowski-embedded D7-brane to the AdS-Schwarzschild horizon. The rest 
mass of the quark is given by the length of the string times the string tension. The rest 
mass thus differs from the Lagrangian mass m. Defining Am = ^ \/AT, the rest mass M^est 
is, for m 3> Am [36] , 

Mrest = m- Am + O (Am^/m) (5.2.7) 

For arbitrary m, the rest mass also differs from the kinetic mass, Mkm, which will be the 
mass that appears in the force law above. More specifically, M^j^ is the mass we must use 
if we write p in terms of v. (Notice we used the kinetic mass M in section 5.1, around eq. 
(5.1.4).) If m 2> Am, then the kinetic mass and thermal rest mass only differ by [36] 

Mkin = Mrest + 0{^m^/m) (5.2.8) 

Our point is that the large- mass limit m — oo actually means not just m ^ T, but 
m 3> Am. In this limit, we will drop all corrections and equate Mfei„ = Mrest = 'm- 

Analysis of single-string solutions [36] revealed that in this limit the correct relation 
between velocity and momentum is relativistic, that is, 
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with 7 the usual rclativistic factor 7 = -j^^. In the second equahty we have invoked the 
m — > 00 Umit and replaced Mjtj„ with m. For a single string moving under the influence of 
the D7-brane electric field, in the m S> Am limit, the result for the drag coeffcient is [36,37] 

|ldm = ^^^\T'^ (5.2.10) 

We will study a finite density (J*) of strings, with (J*) on the order of NfNc, in the 
presence of external E and B fields, and calculate /i^m in the m ^ 00 limit. Using the 
relativistic relation eq. (5.2.9), we will find exactly the same result as the single-string 
calculation. The reason why we find the same answer comes from the probe limit. 

In the probe limit, the plasma contains order adjoint degrees of freedom and order 
NfNc <■ N^ flavor de grees of freedom. The flavor excitations are thus dilute in the large- 
Nc limit. In a perturbative analysis, the flavor excitations will be more likely to scatter 
off of adjoint degrees of freedom than off of other flavor excitations. Scatterings with 
adjoint degrees of freedom will thus be the flavor excitations' primary mechanism for the 
microscopic energy loss that results in the macroscopic drag force. Introducing a density (J*) 
of order NfNc will not change this to leading order in large- AT^. Increasing the stength of the 
coupling muddies the picture of isolated scatterings but does not affect the argument, which 
relies only on laigc-Nc counting. Taking m — > 00, and in particular m S> /U (the chemical 
potential, not the drag coefficient), serves only to dilute the charge carriers further. We 
therefore expect to recover the zero-density result at leading order in the Nf <^ Nc limit. 

The E and B independence follows from this, simply because the zero-density result 
was already, curiously, independent of the quasi-particle momentum, or equivalently 
of m and v [36,37]. As v is determined by E and B, and is the only place where E and B 
could appear in the answer, we expect the answer to be independent of E and B. 

The m — 00 limit is also where we will reproduce the classical result eq. (5.1.4), which 
we derived assuming the charge carriers were quasi-particles. 

We now turn to the details of the supergravity calculation. 
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5.3 The D7-brane Solution 

In this section, we will present the solution for the D7-brane fields describing a field theory 
state with nonzero m, T, (J*), E and B. 

We will use an AdS-Schwarzschild metric in Fefferman- Graham coordinates, as in eq. 
(2.2.6), 

Recall that the boundary is at z = and the black hole horizon is at z = zh with z]j^ = -^T. 
We will not actually need an explicit form for the metric in most of what follows, since most 
of our arguments are based on general properties of AdS-Schwarzschild, in particular that 
gtt vanishes at the horizon and that the metric has a second-order pole at the boundary. 
We will use an S*^ metric 

dsls = d9^ + sm^ edsli + cos^ 9 ds^s. (5.3.12) 

The embedding of the D7-brane will be described by 0{z). We will not explicitly solve for 
0{z), but we will consider limits. We will be considering only black hole embeddings (as we 
want finite density in the field theory), so m = is represented by 6{z) = 0. When m ^ oo 
we expect the D7-branc to form a "spike," as explained in section 4.2.1. In this case, we 
may approximate 9{z) « 7r/2 or cos^(2;) « 0, at least for z along the spike. In particular, 
we will make this approximation for z near the horizon. 

In the field theory, we want a density (J*), external fields E and B, and the resulting 
induced currents (J^) and (J^). As explained in sections 4.3.1 and 4.3.2, we thus introduce 
worldvolume gauge field components At(z) and 

A^{z,t) = -Et + U{z), Ay{z,x)=Bx + fy{z) (5.3.13) 

so that at the boundary we have electric and magnetic fields F*^ = E and F^'^ = B. As 
part of our gauge choice we take = 0. 

As our gauge fields only depend on {z,t,x,y), the D7-brane action is simply a (3+1)- 
dimensional Born-Infeld action, with some "extra" factors in front from the and the 
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extra spatial direction, which may be written as 



S = -M j d^xcos^ Ogl'f\j-g - {2-Ka'Y]^gF^ - {2T^a'Y^ (F A Pf (5.3.14) 

In contrast to our previous definition of S, here we define S as Sd7 divided by a single factor 
of the volume of M. We have also defined d'^x = dz dt dx dy, and defined g = g^z gtt g'ix ^ 
the determinant of the induced metric in the {z,t,x,y) subspace, with gzz = ^/z^ + 9'{zf. 
Writing = F'^'^F^^j, where Greek indices run over {z,t,x,y), and F'^^ = \e^''''^^ F^f^ for 
totally antisymmetric e^^'^l^ with = +1^ we have explicitly 

^5^^ = glxA't + gugxxA!^ + gugxx^'^ + gzzgxxAl + gzzgttAl (5.3. 15a) 

i (F A Ff = Qf'^^F^.) ' = AlA? + AlA'y' + 2AyA[A,A'y. (5.3.15b) 

where dots, A, denote derivatives with respect to t, primes, A', denote derivatives with 
respect to z, and bars. A, denote derivatives with respect to x. 

The action only depends on the derivatives of At{z), fx{z) and fy{z), so we will have 
three conserved charges. In the Appendix to this chapter we identify these as (J*), (J^) 
and {jy), 

M{27ra'fgii' cos3 9 -dlxA't - {27r a' f{AlA[ + Ay A, A ^ ^ ^^^^ ^^^^^^^ 
^-g - {2T,a' f\gF^ - {2T,a'Y\ (F A Ff 

M{27ra'fgl/f cos^ 6 , \fJn\g.r.,A: ^ ^ (5.3.16b) 

^-g - {2'Ka'f\gF^ - {2'Ka'Y\ (F A Ff 

,oi/o Q \Qtt\QxxAi, - (27ra'f(AlAi, + AyA:i:A'A 
Af{2Tra'fgl/f cos^ 6 ^ ^ ^ ) K x y x t) ^ ^jy^ (5.3.16c) 



-g - {2TTa' f\gF^ - {2TTa'Y\ (F A F)' 
With a little algebra we solve for the gauge fields from eq. (5.3.16), 

4(.) = -^^i4^i^ (5.3.17) 

where we have introduced the coefficients 

i = Igttlglx-C^^c^'fF'^F^ 

= \gu\glx + i27ra' f{\gu\B^-gccxE^) (5.3.18a) 
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X = \gtt\gLN\2T^a'f<^os^e-{27:a'f{J,){Jn 

= |ff«|ffLAA'(2W)4cos6^+ (2W)2 {\gu\{J'f-g.. {{r? + (J^)2)X5.3.18b) 

a = -(2W)2f^'^(J^) 

= {2Tra')Wgtt\{J')B + g^^{jy)E) (5.3.18c) 

Notice that ^ is simply —det{gab + {2Tra')Fab) in the {t,x,y) subspace, and that cos0(2:) 
appears only in %. We also have 

^U^) = M^^, 4(-) = ,/^i:9i±S (5.3.19) 

In the original action we may now replace the gauge fields with the conserved charges. The 
resulting effective action has only the single dynamical field 0{z), 

S = -J^ii-Kolf J d^xcos^ 0gl^,/^;^\-—i== (5.3.20) 

We may obtain the equation of motion for 9{z) in two ways. Wc may derive it from the 
original action eq. (5.3.14) and then plug in our gauge field solutions eqs. (5.3.17) and 
(5.3.19), or we may Legendre transform to eliminate the gauge fields at the level of the 
action. The Legendre-transformed action S is 

S . ,5.3.21) 

= - (27^^)2 / d'^x9li'^\9tt\~^^^gxxV^X - 

where = A!f{z), = A'^{z) and -0^ = A'y{z) reproduce eqs. (5.3.17) and (5.3.19). 

Specifying the boundary conditions will then determine the D7-brane solution com- 
pletely. At the horizon, At{z) must of course obey At{zH) = 0. We are free to choose 
the leading asymptotic values of the fields near the boundary z — > 0. We first choose the 
asymptotic value Oq of 9{z). The gauge fields asymptotically approach the boundary as 

A,(!,) = -Et + c^ + \jj!^^^z'' + 0(z*) (5.3.22b) 
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Ay{z) = Bx + Cy + ^jS^L_z^ + 0{z') (5.3.22c) 

where /Lt, Cx and Cy are constants of integration. The leading asymptotic value fJL is the U (1) b 
chemical potential. For A^. and Ay wc impose the boundary condition c^- = Cj, = because 
we do not want to source the corresponding operators in the field theory Lagrangian. In 
what follows, we will not actually make use of this boundary condition, however, since we 
will only work with the derivatives (or field strengths) A!^{z) and A!y{z). 

5.4 The Conductivity Tensor 

We focus now on the quantity ^/^X~ appearing in the effective action eq. (5.3.20). 
We will find that demanding reality of the effective action allows us to solve for ( J^) and 
(J^), and hence the conductivity, in terms of E, B and (J*). The action must be real: an 
imaginary action signals an instability of the solution. 

In eq. (5.3.18a) we see that, as a function of z, ^ has a zero: ^ < at the horizon where 
\gtt\ = 0, whereas ^ > near the boundary z ^ 0. We denote the zero of ^ as z*, for which 
we can solve (here we must use the explicit form of the metric in eq. (5.3.11)), 
z^ 



^ = e^-h^ + ^(e2 - 62)2 + 2(e2 + b'^) + l (5.4.23) 



z^ 



-J ((e2 - 62) + ^(e2 _ 62)2 + 2{e^ + 62) + _i 



where we have defined the dimensionless quantities 

1 E \ B 
e = -(27ra')Ezl = — b = -(2Tra')Bzl = ^— (5.4.24) 

and converted to field theory quantities. Knowing that ^ is the {t,x,y) part of —det{gab + 
{27Ta')Fab), we will interpret as an effective horizon on the D7-brane worldvolume, or in 
other words, as the location in z where the determinant of the D7-brane action vanishes. 
Notice also that when £^ = we find = zh, so in this case the effective horizon coincides 
with the actual horizon of the background geometry. As wc increase E, then, wc should think 
of Zif as starting at the horizon and then moving up the D7-brane, toward the boundary. 
We will also need Qxxi^*) = 7r^r^jr(e, 6) where 

^(e, 6) = i (l + e2 - 62 + ^(e2 - 62)2 + 2(e2 + 62) + l) (5.4.25) 
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For later use notice that J-'{0, b) = 1. 

In fact all three functions, ^, x and a must share the same zero z*. Prom eq. (5.3.18b) 
we see that at the horizon x < while at the boundary % > 0, so % also has a zero. In 
particular > at the horizon and at the boundary. If ^ and % have distinct zeroes, then 
in the region between those zeroes one would change sign while the other would not, hence 
in that region < and the effective action would be imaginary. The only consistent 
possibility is for ^ and x to share the zero at z*. We must also have < ~^ as z ^ z*, 
so that a ^ at as well. 

We thus set all of eqs. (5.3.18) to zero at and solve for (J^) and (J^), 

9xx "T {^T^o: ) n 

(5.4.26a) 

with all functions of z evaluated at z*. Converting to field theory quantities, we find 



pb 



(5.4.27b) 



(5.4.28) 



where we have defined p similarly to e and b, 

but while e and b are dimensionless, p has dimension one. 

We can interpret our result as follows. Two types of charge carriers contribute to the 
conductivity. The first are the charge carriers we have introduced explicitly in p. Taking 
p = leaves a nonzero axx, however, so we must have another source of charge carriers. We 
will guess that these come from pair production in the plasma: although the net charge p is 
zero, the flavor fields are still part of the plasma, just with equal numbers of particles and 
anti-particles. When we introduce the electric field E, the particles go one way and the anti- 
particles go the other way, producing a net current. Notice that the hypermultiplet mass 
m is hiding in cos^(z*), whose behavior is consistent with our interpretation: cos6{z*) — 
as m ^ oo, and cos 9{z^) = 1 for m = 0. 
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Notice also that our answer depends on the electric field E. We are therefore capturing 
effects beyond linear response: in linear response theory, the conductivity would be constant 
(as the current would be linear in the electric field: (J^) = axxE, etc.). We can attribute 
the nonlinearities of our answer to the fact that we started with a nonlinear action, the 
Born-Infeld action. As explained in section 2.3.2, however, we can linearize the Born-Infeld 
action, that is, we can Taylor expand in the gauge field and recover the Maxwell action as 
the leading term. We can thus recover the linear-response answer by linearizing in E, which 
in practical terms means simply setting E = in our result. 

We will now check our answer in two limits, to reproduce the forms in section 5.1. 

To recover cq. (5.1.2), wc must linearize in the electric field, setting e = 0, and hence 
J^{0,b) = 1, in eq. (5.4.27). We also restore Lorentz invariance by taking T ^ 0. We find 
(^xx = and axy = {J^)/B, as expected. 

To recover eq. (5.1.4), we return to finite T and again linearize in the electric field. We 
additionally take the m — > oo limit cos9{z^) ^ 0. The conductivity becomes 



- J^- (5-4.29) 



62 + 1' 52 + 1- 

As shown in section 5.5, in the m ^ oo limit we identify ^y/XT"^ = Hafn, where /U^ is the 
drag coefficient. We thus have p = and h = so the conductivity indeed has the 

form expected for quasi-particles propagating through an isotropic, homogeneous medium, 
eq. (5.1.4). 

We end with a question: where in the phase diagram of the field theory is our result 
valid? The short answer is in supergravity language. Our argument relied on the fact that 
the D7-brane touches the AdS-Schwarzschild horizon, such that the quantities under the 
square root in the D7-brane action change sign at some z* . Our solution should thus be valid 
whenever the ground state of the system is a black hole embedding of the D7-brane. Notice 
in particular that our answer will be valid at zero density, but only in the high-temperature, 
melted-meson phase, described by black hole embeddings. 
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5.5 The Drag Force 

In this section we compute the drag force on the charge carriers, in the large-mass limit. In 
particular, we will compute the product fidm. 

In the m ^ oo limit where cos ^ ~ 0, we expect the flavor excitations to be well-described 
as a collection of quasi-particles, which will have an equation of motion 

^=-l^dP + E + vxB, (5.5.30) 



with V is the quasi-particle velocity and fi^ is the drag coefficient. Our first goal is to 

dp 
dt 

i\p\ = ^Je^ + v^B^ + 2E-{vxB) (5.5.31) 



compute the magnitude of the drag force, fidlPl- In the steady-state, ^ = 0, which implies 



m 

As m — DO, we expect pair creation to be suppressed, so only the charge carriers in (J*) 
should contribute to (J), hence (J) = {J^)v. We immediately read off v'^ = \gtt\/gxx by 
setting X to zero at and dropping the cos^(2;*) term. Setting ^ = at gives us 

Setting a = at gives us the component of v in the y direction, 

(J«) \ga\B .,B 

We then have 2E ■ {v x B) = 2EBvy = -2B^v^. The drag force is then 

fJ-dlpl = ^~^9xxiz*)v (5.5.34) 

We can now compute /Xdm. We employ the relativistic relation \p\ = jmv with 7 = ■-j^=, 
and find 

l^d'm= - — j^/gxxiz*)"^ - \9tt{z*)\9xxiz*) = — jz]^^ = -Vxt"^ (5.5.35) 

This result for ^4 m is identical to the drag force computed from a single-string solution in 
ref. [36], as advertised. In other words, the value of /x^m for a single quark or squark, when 
(J*) = and = 0, is identical to eq. (5.5.35), as expected. 
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Appendix: Holographic Renormalization of Worldvolume Gauge Fields 

In this section we study the holographic renormahzation of the D7-brane's worldvolume 
fields, in particular the gauge fields introduced in section 5.3. Our main goal is to identify 
the constants of motion in eqs. (5.3.16) as the expectation values (J^). 

We find from its equation of motion that 9{z) has the usual asymptotic expansion 

e{z) = 9oz + 92Z^ + 0{z^). (1) 

where the leading coefficient 9o is the source for the dual operator, i.e. gives the hypermul- 
tiplet mass via = (27ra')m. 

Plugging eq. (1) into the regulated action we find, with nonzero E and B, the divergences 

dzL = -A/'y dz [z-^ - 9lz-^ + ^{2'Kolf{B'^ - E^)z-^ + 0{z)\ (2) 

We again need the counterterms 

Li = L2 = Lf = M^V^9{er (3) 

with the induced metric at z = e and 7 its determinant. Notice that ^—7 = + 
O(e^). We suppress J dtdxdy unless stated otherwise. The last divergence requires a new 
counterterm 

Lf = ~^^{2^^a'f^/^F'^Fij log € = -^Mi27ra'fiB^ - E^)loge + Oie^ log e) (4) 

The generating functional of the field theory is then the e — limit of S = Sreg + Yli ^i- 
We want the expectation values (J*), (J^) and (J^). In holo-rg, (J^) is 

(J'^) = limi-L-^ (5) 

For (J*), we need 

fZH XT XT rZH XT 

where we have used the fact that is z-independent on-shell. Enforcing 5Ai{zH) = 

we find = 3^ and hence (J*) = J^^, as advertised. 
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For (J^), we reinstate / dt because Ax is time-dependent, 

We employ precisely the same argument as before for the first term. For the second term 
we observe that is t-independent on-shell and hence 



dtJ^^dtSA, = / dtdtSAx = (8) 



/-«,..-^-^ .... 

where we demand that the fluctuation be weh-behaved (vanishing) at t = ±00. The coun- 
terterm Lf gives a vanishing contribution to (J^) for the same reason, 

SLf = -^M{2na')^^Y'l'' J dtj^{FikFji)dt6Axloge (9) 

= +^M{27ra'f j dtAx{e)dt5Ax log e + 0(e^ log e) 

= ©(e^loge) 

We then have = and hence (J^) = ^1^. 

{jy) is very similar. Ay depends on x so we reinstate J dx. We have 

"'"^ [wa^'-'^' + 

The same argument as above applies for the first term, and for the second term we observe 
that j§^j^ is x-independent on-shell. Demanding that the fluctuation be well-behaved at 
X = ±00 gives J" dxdxSAy = and no contribution from Lp. We thus have (J^) = j§^- 
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Chapter 6 
CONCLUSION 

We have covered much ground, so let us summarize the sahent points to take away from 
ah of this. We studied the holographic dual of a strongly-coupled non-Abelian gauge theory, 
focusing on the thermodynamics and transport properties of the fields in the fundamental 
representation of the gauge group. Specifically, wc studied the theory at zero temperature 
and finite baryon number density, where, using the holographic description, we discovered a 
second-order phase transition (in the grand canonical ensemble) when the chemical potential 
was equal to the mass of the flavor fields. The transition is between a zero-density phase 
with a gapped, discrete spectrum of stable mesons and a finite-density phase with a gapless, 
continuous spectrum of unstable mesons. At finite temperature, we introduced external 
electric and magnetic fields and, again using the holographic description, computed the 
resulting currents of fiavor fields, from which wc extracted the conductivity tensor. As an 
added bonus, we computed the drag force acting on the fiavor fields. 

Much work remains to be done. The most obvious task is to complete the phase diagram 
of the theory in the full parameter space of m, T, (J*), E and B, in both the canonical and 
grand canonical ensembles. The D7-brane solutions discussed above are almost certainly in- 
sufficient. As mentioned in section 4.2.3, these D7-brane solutions all describe homogeneous 
phases in the gauge theory, whereas we have good reasons to expect that mixed phases may 
play a pivotal role in certain regions of the phase diagram. As mentioned in section 4.2.4, 
additional homogeneous phases may also be present: at low enough temperature and high 
enough chemical potential, we expect the hypermultiplet scalars to undergo Bose-Einstein 
condensation, for instance. The D7-brane solutions above do not describe such a conden- 
sate. Some regions of the phase diagram are currently not understood at all, for example 
the regions described by singular D7-brane emebeddings, as mentioned in section 4.3.1. Ad- 
ditionally, the complete meson spectrum (and/or spectral functions) should be computed 
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in the full parameter space, to provide a complete picture of the physics in any phase. And 
another obvious direction to go is to continue studying the transport properties of the flavor 
fields. For example, what current is generated by a thermal gradient (the Nernst effect)? 

We argued in the introduction that gauge-gravity duality may teach us something about 
strongly-coupled non-Abelian gauge theories. More recently, however, a new and exciting 
application of gauge-gravity duality has emerged. Many condensed matter systems may 
be described by strongly-coupled, scale-invariant field theories, for example cold atoms at 
unitarity or high- Tc superconductors near their quantum critical points. Gauge-gravity 
duality has already provided toy models for such systems, and promises much more [118- 
124]. Only recently (as of the time of this writing), gauge-gravity duality has even described 
superconductivity [125-127]! 

One key feature noticably lacking in these models, however, is a Fermi surface. One 
big obstacle to realizing a Fermi surface holographically is supersymmetry: gauge-gravity 
systems usually describe fermionic and bosonic charge carriers. Breaking supersymmetry 
means we are not guaranteed stability, however. To find a holographic description of a finite 
density of fermionic charge carriers, then, we will probably have to break supersymmetry 
and then add some other ingredient to ensure stability. 

We hasten to point out, however, that looking for a Fermi surface is a fool's errand: 
generically, the Fermi surface is unstable. For fermions at low temperature and high den- 
sity, where we expect a Fermi surface, the appropriate degrees of freedom are excitations 
about the Fermi surface. Standard renormalization arguments then reveal that, for exci- 
tations with equal but opposite momenta, under renormalization group flow toward the 
Fermi surface, any repulsive interaction becomes irrelevant while any attractive interaction 
becomes relevant. As a result, fermions on opposite sides of the Fermi surface bind into 
pairs, called Cooper pairs, and the resulting scalar condensate breaks the associated U{1) 
symmetry (or gauge invariance) . This is of course the BCS theory of superconductivity, ap- 
plicable to electrons in a low-temperature metal or quarks in massless, three-flavor QCD at 
asymptotically high density. Good explanations of BCS theory appear in refs. [6,128-130]. 
Some canonical papers about BCS theory in high-density QCD are refs. [131-135]. 

The question we should ask is thus not, "What is a holographic description of a Fermi 
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surface?" but rather the more subtle question, "What is the holographic description of the 
ground state of a low- temperature, finite-density system of strongly-interacting fermions?" 
Any insights that holography can provide into such systems would be valuable for both 
condensed matter physics and QCD. 

We hasten to point out still another caveat: the large- A?c limit changes the physics signifi- 
cantly. For massless QCD at low temperature and high density, and in the 't Hooft limit, the 
ground state of the system may actually be spatially inhomogcncous, forming a so-called 
"chiral density wave," analogous to similar plane-wave ground states in low-dimensional 
electron systems [136, 137]. In practical terms, this would require, for example, studying 
D7-brane embeddings that depend on Minkowski directions. These are typically very diffi- 
cult to construct. The rewards for doing so, however, may provide sufficient motivation for 
us to "shut up and calculate." 

Gauge-gravity duality is clearly a versatile tool. We hope only that, in this dissertation, 
we have convinced the reader that gauge-gravity duality is useful for studying strongly- 
coupled gauge theories, as a tool on par with and complementary to other methods, such as 
lattice simulation. We believe that holographic methods are worth learning, as they should 
prove fruitful for a long time to come. 
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